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
41       real*8 x1,x2  variable
61         x1=xr(0*p+i)
63         call feval(x1,x2,val)
78         x1=xq(0*q+i)
80         call dfeval(x1,x2,val)
88         call dfeval(x2,x1,val)