fortran 做 函数图像 r=a*(1+cos(b))