PROGRAM EXA_11
      C     USE DISLIN         for Fortran 90!
            CHARACTER*60 CTIT1,CTIT2
            EXTERNAL ZFUN
            CTIT1='Surface Plot (SURFUN)'
            CTIT2='F(X,Y) = 2*SIN(X)*SIN(Y)' 
            CALL SETPAG('DA4P')
            CALL DISINI
            CALL PAGERA
            CALL COMPLX
            CALL AXSPOS(200,2600)
            CALL AXSLEN(1800,1800)
            CALL NAME('X-axis','X')
            CALL NAME('Y-axis','Y')
            CALL NAME('Z-axis','Z')
            CALL TITLIN(CTIT1,2)
            CALL TITLIN(CTIT2,4)
            CALL VIEW3D(-5.,-5.,4.,'ABS')
            CALL GRAF3D(0.,360.,0.,90.,0.,360.,0.,90.,-3.,3.,-3.,1.)
            CALL HEIGHT(50)
            CALL TITLE
            CALL SURFUN(ZFUN,1,10.,1,10.)
            CALL DISFIN
            END     
            FUNCTION ZFUN(X,Y)
            FPI=3.14159/180.
            ZFUN=2*SIN(X*FPI)*SIN(Y*FPI)
            END
 Figure B.12: Surface Plot