real theta(512), singlebeam_theta(512), phi(1024), singlebeam_phi(1024)








        call pgbegin(0,'?',2,1)
        call pgslw(5)
        call pgpap(0.0,0.5)
        call pgadvance
        call pgsch(2.5)
        call pgscf(2)
        call pgvport(0.25,0.9,0.25,0.9)
        call pgwindow(0.,90.,0.,1.15)
        call pgsci(3)
        call pgbox('BCTNS',0.,0.,'BCTNSV',0.,0.)
        call pglabel('\(2185)(deg)','', '')
        call pgsch(1.0)
        call pgsci(3)
        call pgline(512,theta,singlebeam_theta)

        call pgadvance
        call pgsch(2.5)
        call pgscf(2)
        call pgvport(0.25,0.9,0.25,0.9)
        call pgwindow(0.,180.,0.,1.15)
        call pgsci(3)
        call pgbox('BCTNS',0.,0.,'BCTNSV',0.,0.)
        call pglabel('\(2147)(deg)','', '')
        call pgsci(3)
!        call pgpoint(512,phi,singlebeam,12)
        call pgline(1024,phi,singlebeam_phi)
        call pgend

        end
这类的程序 不太会 那位高人懂这个的帮讲解一下呗 想绘制f(x,y)=cos((2*3.14159*0.832*(-0.093*sin(x)*cos(y)+(-7.709+7.709)*sin(x)*sin(y))+23.124)/2)的图象 不知中间空白处该怎么编