How to find skewness, kurtosis?

I have two vectors, charge 'q' and corresponding phase angle 'Ph', each of size 89905x1. 'Ph'- min=0 degree, max=360 degree & 'q'- min=-4e-12, max=4.5e-12. I have ploted these data, in which the x-axis is my phase angle vector which is divided into small phase windows of width,say 1.5 degree and y axis are: 1. number of charges present in that phase window. 2. sum of charge values observed in phase window. 3. average value of charge value observed in phase window. 4. maximum value of charge observed in phase window. How can I find skeness, kurtosis of this distributions ? How can I get a distribution fit of my plots? How can I find skeness, kurtosis of this distribution fits ?