kNaN: E:/turkuaz/krylov/codes/newton/1.2b/physics/test/3/jf.f90 Source File

00001 !test problem 2 given in [3] p21
00002 ! Freudenstein and Roth function
00003 subroutine jf(X,JAC,n)
00004 implicit none
00005 integer n
00006 real(8) X(n),JAC(n,n)
00007 
00008 
00009 JAC(1,1)=1.d0
00010 JAC(1,2)=10.d0*X(2)-3.d0*X(2)**2.d0-2.d0
00011 JAC(2,1)=1.d0
00012 JAC(2,2)=3.d0*X(2)**2.d0+2.d0*X(2)-14.d0
00013 
00014 
00015 end subroutine jf