shell.m

Although most of the health problems related to fast food aren t felt until middle age -- obesity an

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  shell.m - evaluates nonlinear sizing equations and the
%            associated Jacobian for spshell.m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [f,J] = shell(x);
data1sh,
f=zeros(3,1); J=zeros(3,3);
f(1)=2*Ggammam*Gpm*x(1)*x(3)*GSF*(x(3)+(Grw+1)*x(2))-GPm;
f(2)=2*Grhoc*Gdelc^2*Gkc*Grw*x(2)^2*(x(1)+x(2)+x(3))-Gkl*GPm;
f(3)=4.44*GN*Gf*GSF*GBm*x(1)*x(3)/(39.37)^2-GV1;
J(1,1)=2*Ggammam*Gpm*GSF*(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*(2*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);