I Have used Sage command
sage: S = SBox(14,4,13,1,2,15,11,8,3,10,6,12,5,9,0,7) sage: S.difference_distribution_table() the output is [16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0] [ 0 0 0 2 0 0 0 2 0 2 4 0 4 2 0 0] [ 0 0 0 2 0 6 2 2 0 2 0 0 0 0 2 0] [ 0 0 2 0 2 0 0 0 0 4 2 0 2 0 0 4] [ 0 0 0 2 0 0 6 0 0 2 0 4 2 0 0 0] [ 0 4 0 0 0 2 2 0 0 0 4 0 2 0 0 2] [ 0 0 0 4 0 4 0 0 0 0 0 0 2 2 2 2] [ 0 0 2 2 2 0 2 0 0 2 2 0 0 0 0 4] [ 0 0 0 0 0 0 2 2 0 0 0 4 0 4 2 2] [ 0 2 0 0 2 0 0 4 2 0 2 2 2 0 0 0] [ 0 2 2 0 0 0 0 0 6 0 0 2 0 0 4 0] [ 0 0 8 0 0 2 0 2 0 0 0 0 0 2 0 2] [ 0 2 0 0 2 2 2 0 0 0 0 2 0 6 0 0] [ 0 4 0 0 0 0 0 4 2 0 2 0 2 0 2 0] [ 0 0 2 4 2 0 0 0 6 0 0 0 0 0 2 0] [ 0 2 0 0 6 0 0 0 0 4 0 2 0 0 2 0] This output is produced from one first line of S box. How to find the difference table of the whole S box