%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  core.m - evaluates nonlinear sizing equations and the
%           associated Jacobian for spcore.m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [f,J] = core(x);
data1cor,
f=zeros(3,1); J=zeros(3,3); 
f(1)=2*Ggammam*Gpm*x(1)*x(3)*GSF*(2*x(3)+(Grw+1)*x(2))-GPm;
f(2)=2*Grhoc*Gdelc^2*Gkc*Grw*x(2)^2*(x(3)+x(1)+x(2))-Gkl*GPm;
f(3)=4.44*GN*Gf*GSF*GBm*x(1)*x(3)/(39.37)^2-2*GV1;
J(1,1)=2*Ggammam*Gpm*GSF*(2*x(3)^2+(Grw+1)*x(2)*x(3));
J(1,2)=2*Ggammam*Gpm*GSF*((Grw+1)*x(1)*x(3));
J(1,3)=2*Ggammam*Gpm*GSF*(4*x(1)*x(3)+(Grw+1)*x(1)*x(2));
J(2,1)=2*Grhoc*Gdelc^2*Gkc*Grw*x(2)^2;
J(2,2)=2*Grhoc*Gdelc^2*Gkc*Grw*(2*x(2)*x(3)+2*x(1)*x(2)+3*x(2)^2);
J(2,3)=2*Grhoc*Gdelc^2*Gkc*Grw*x(2)^2;
J(3,1)=4.44*GN*Gf*GSF*GBm/(39.37)^2*x(3); 
J(3,2)=0; 
J(3,3)=4.44*GN*Gf*GSF*GBm/(39.37)^2*x(1);