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
