program scrsho_main
       real fx
       external fx
c
       call scrsho(fx)
c
       end



       function fx(a)
c
c      fx(a) = a**2
c
       real fx, a
c      fx = a**3
       fx = 7*sin(a) + a
       end