E:/turkuaz/krylov/codes/newton/1.2b/physics/test/1/jf.f90

Go to the documentation of this file.

!example problem 6.5.1 given in [1] p149
subroutine jf(X,JAC,n)
implicit none
integer n
real(8) X(n),JAC(n,n)


JAC(1,1)=2.d0*X(1)
JAC(1,2)=2.d0*X(2)
JAC(2,1)=dexp(X(1)-1.d0)
JAC(2,2)=3.d0*X(2)**2.d0


end subroutine jf

---