Computing skewness and other distributional Characteristics when pdf is known

Hi everybody
I want to compute skewnesss and other distributional characteristics of a gamma distribution when its shape and scale parameters are known (pdf is known, not through its variates/observations). How can I do this through R, or other software?

Thanks for your reply.