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
62         x2=xr(1*p+i)
63         call feval(x1,x2,val)
79         x2=xq(1*q+i)
80         call dfeval(x1,x2,val)
88         call dfeval(x2,x1,val)