7       subroutine feval(x1,x2,val)  argument
8       real*8 x1,x2,val  local
10       val=x1*x1+x2*x2+x1*x2+1
14       subroutine dfeval(x1,x2,val)  argument
15       real*8 x1,x2,val  local
17       val=2*x1+x2
40       real*8 val,diff  variable
63         call feval(x1,x2,val)
64         iinput(i)=val
80         call dfeval(x1,x2,val)
81         diff=val-ooutput(0*q+i+ooffset)
85      &     ' != ',val
88         call dfeval(x2,x1,val)
89         diff=val-ooutput(1*q+i+ooffset)
93      &     ' != ',val