function J = jacobian(y)
F = y(1);
B = y(2);
Y = y(3);
J = y(4);
C = y(5);
K = y(6);
A = y(7);
P = y(8);

J = zeros(8,8);

%G1 = F';
%J(1,1) = G1_F
J(1,1)= 10*(1-2*F) - 2*B - 2*Y;
%J(1,2) = G1_B
J(1,2)=  - 2*F;
%J(1,3) = G1_Y
J(1,3) = - 2*F;


return

