subroutine fdjac1(fcn, n, x, y)
integer :: n, x, y
call fcn(x, y)
end
