subroutine jfcn(t, y, jf)
use work
! Calculeaza membrii II ai ec. dif. de ordinul I
implicit real(8) (a-h,o-z)
real(8) y(n), ff(n), jf(n,n)


!ff(1) =y(2)
!ff(2) =lambda*(1-y(1)**2)*y(2) -y(1) + A*dcos(omega*t)

! fcn(1) = c(0)*y(1) -h*ff(1)  +c(1)*xa(0,:) +c(2)*xa(-1,:) +c(3)*xa(-2,:) +c(4)*xa(-3,:) +c(5)*xa(-4,:)
! fcn(2) = c(0)*y(2) -h*ff(2)  +c(1)*xa(0,:) +c(2)*xa(-1,:) +c(3)*xa(-2,:) +c(4)*xa(-3,:) +c(5)*xa(-4,:)

jf(1,1) =c(k,0)	; jf(1,2) =-h;
jf(2,1) =-h*(-lambda*2*y(1)*y(2) -1)	; jf(2,2) =c(k,0) -h*(lambda*(1-y(1)**2))	;

END subroutine jfcn

