老师要求用fortran编程实现GPD分布函数的参数,求代码