function [residual, g1, g2, g3] = CA_BMZ12_dynamic(y, x, params, steady_state, it_)
%
% Status : Computes dynamic model for Dynare
%
% Inputs :
%   y         [#dynamic variables by 1] double    vector of endogenous variables in the order stored
%                                                 in M_.lead_lag_incidence; see the Manual
%   x         [nperiods by M_.exo_nbr] double     matrix of exogenous variables (in declaration order)
%                                                 for all simulation periods
%   steady_state  [M_.endo_nbr by 1] double       vector of steady state values
%   params    [M_.param_nbr by 1] double          vector of parameter values in declaration order
%   it_       scalar double                       time period for exogenous variables for which to evaluate the model
%
% Outputs:
%   residual  [M_.endo_nbr by 1] double    vector of residuals of the dynamic model equations in order of 
%                                          declaration of the equations.
%                                          Dynare may prepend auxiliary equations, see M_.aux_vars
%   g1        [M_.endo_nbr by #dynamic variables] double    Jacobian matrix of the dynamic model equations;
%                                                           rows: equations in order of declaration
%                                                           columns: variables in order stored in M_.lead_lag_incidence followed by the ones in M_.exo_names
%   g2        [M_.endo_nbr by (#dynamic variables)^2] double   Hessian matrix of the dynamic model equations;
%                                                              rows: equations in order of declaration
%                                                              columns: variables in order stored in M_.lead_lag_incidence followed by the ones in M_.exo_names
%   g3        [M_.endo_nbr by (#dynamic variables)^3] double   Third order derivative matrix of the dynamic model equations;
%                                                              rows: equations in order of declaration
%                                                              columns: variables in order stored in M_.lead_lag_incidence followed by the ones in M_.exo_names
%
%
% Warning : this file is generated automatically by Dynare
%           from model file (.mod)

%
% Model equations
%

residual = zeros(90, 1);
T11 = y(47)/params(59);
T19 = y(42)/params(54);
T40 = y(50)/params(72);
T47 = y(50)/y(63);
T203 = (1-params(35))*y(50)*y(39)/y(2);
T254 = y(49)^params(51);
T268 = params(41)/2;
T270 = T268*(y(36)/y(2)-params(37))^2;
T271 = y(2)*T270;
T273 = params(43)/(2*T254);
T274 = y(36)^2;
T281 = (y(51)*y(35))^params(35);
T282 = y(2)^(1-params(35));
T288 = params(39)*y(35)^params(38);
T302 = params(42)*y(42)^params(40);
T357 = y(76)/params(35);
T360 = y(46)*params(35)/((1-params(35))*y(76));
T361 = T360^(1-params(35));
T368 = y(39)*y(50)*(1-params(35))/y(2);
T373 = y(51)^params(35);
T379 = y(45)*y(37)/y(41);
T381 = T379^params(36);
T408 = params(40)/(params(40)-1);
T417 = (1-params(42)*y(42)^(params(40)-1))/(1-params(42));
T425 = params(34)*params(42)*y(122)/y(38);
T426 = y(123)^params(40);
T427 = T425*T426;
T432 = y(123)^(params(40)-1);
T433 = T425*T432;
T481 = T268*(y(55)/y(10)-params(37))^2;
T489 = (y(51)*y(54))^params(35);
T490 = y(10)^(1-params(35));
lhs =y(94);
rhs =log(T11)*400;
residual(1)= lhs-rhs;
lhs =y(96);
rhs =400*log(T19);
residual(2)= lhs-rhs;
lhs =y(95);
rhs =0.25*(y(96)+y(15)+y(18)+y(19));
residual(3)= lhs-rhs;
lhs =y(98);
rhs =log(T40)*100;
residual(4)= lhs-rhs;
lhs =y(97);
rhs =100*log(T47);
residual(5)= lhs-rhs;
lhs =y(93);
rhs =0.5*(0.01*y(97)^2+(y(95)-y(96))^2);
residual(6)= lhs-rhs;
lhs =y(94);
rhs =params(32)*x(it_, 5)+params(31)*y(139)+params(30)*y(138)+params(29)*y(137)+params(28)*y(130)+params(27)*y(22)+params(26)*y(21)+params(25)*y(20)+params(24)*y(17)+y(98)*params(23)+params(22)*y(136)+params(21)*y(135)+params(20)*y(134)+params(19)*y(129)+params(18)*y(25)+params(17)*y(24)+params(16)*y(23)+params(15)*y(16)+y(97)*params(14)+params(13)*y(133)+params(12)*y(132)+params(11)*y(131)+params(10)*y(128)+y(15)*params(6)+y(96)*params(5)+params(1)*y(14)+params(2)*y(26)+params(3)*y(27)+params(4)*y(28)+params(7)*y(18)+params(8)*y(19)+params(9)*y(29);
residual(7)= lhs-rhs;
lhs =y(38);
rhs =y(43)/y(33);
residual(8)= lhs-rhs;
lhs =y(38);
rhs =y(47)*params(34)*y(122)/y(123);
residual(9)= lhs-rhs;
lhs =y(46);
rhs =(T203*y(40)+(1-params(37))*y(45))/y(6);
residual(10)= lhs-rhs;
lhs =y(30);
rhs =y(47)/y(123);
residual(11)= lhs-rhs;
lhs =y(85);
rhs =y(123);
residual(12)= lhs-rhs;
lhs =y(124);
rhs =y(30)*y(48);
residual(13)= lhs-rhs;
lhs =y(90);
rhs =y(30)*y(48);
residual(14)= lhs-rhs;
lhs =y(41);
rhs =params(96)*(y(2)*y(46)*y(6)-y(7)*y(1)*(y(2)*y(6)-y(4)));
residual(15)= lhs-rhs;
lhs =y(88);
rhs =y(41)*(1-params(96));
residual(16)= lhs-rhs;
lhs =y(45)*y(49);
rhs =1+params(41)*(y(36)/y(2)-params(37))+y(36)*params(43)/T254;
residual(17)= lhs-rhs;
lhs =y(37);
rhs =y(2)*(1-params(37))+y(49)*y(36);
residual(18)= lhs-rhs;
lhs =y(34);
rhs =y(36)*params(43)/(2*T254);
residual(19)= lhs-rhs;
lhs =y(50);
rhs =y(33)+y(36)+T271+T273*T274;
residual(20)= lhs-rhs;
lhs =y(50)*y(39);
rhs =T281*T282;
residual(21)= lhs-rhs;
lhs =T288;
rhs =y(40)*y(38)*params(35)*y(50)*y(39)/y(35);
residual(22)= lhs-rhs;
lhs =y(39);
rhs =(1-params(42))*y(44)^(-params(40))+T302*y(3);
residual(23)= lhs-rhs;
lhs =y(65);
rhs =log(y(50))-log(params(72));
residual(24)= lhs-rhs;
lhs =y(66);
rhs =log(y(36))-log(params(71));
residual(25)= lhs-rhs;
lhs =y(91);
rhs =log(y(33))-params(39)*y(35)^(1+params(38))/(1+params(38));
residual(26)= lhs-rhs;
lhs =y(92);
rhs =y(91)+params(34)*y(127);
residual(27)= lhs-rhs;
lhs =y(73);
rhs =y(39)*y(50)*params(35)/y(35);
residual(28)= lhs-rhs;
lhs =y(75);
rhs =y(39)*y(50)*(1-params(35))/y(37);
residual(29)= lhs-rhs;
lhs =y(74);
rhs =y(37)/y(35);
residual(30)= lhs-rhs;
lhs =y(76);
rhs =T288/y(38);
residual(31)= lhs-rhs;
lhs =y(77);
rhs =y(76)/y(73);
residual(32)= lhs-rhs;
lhs =y(78);
rhs =y(46)/y(75);
residual(33)= lhs-rhs;
lhs =y(79);
rhs =y(77)+y(78);
residual(34)= lhs-rhs;
lhs =y(80);
rhs =T357*T361;
residual(35)= lhs-rhs;
lhs =y(81);
rhs =y(76)/y(73);
residual(36)= lhs-rhs;
lhs =y(82);
rhs =(y(46)*y(6)-(1-params(37))*y(45))/T368;
residual(37)= lhs-rhs;
lhs =y(83)^params(35);
rhs =y(50)*y(39)/(T282*T373);
residual(38)= lhs-rhs;
lhs =y(48);
rhs =y(68)*T381*y(69);
residual(39)= lhs-rhs;
lhs =y(89);
rhs =T381;
residual(40)= lhs-rhs;
lhs =y(70);
rhs =y(45)*y(37)-y(41);
residual(41)= lhs-rhs;
lhs =y(69);
rhs =(y(71)/params(108))^params(105);
residual(42)= lhs-rhs;
lhs =y(71);
rhs =y(42)*y(70)/y(13);
residual(43)= lhs-rhs;
lhs =y(72);
rhs =y(70)/y(13);
residual(44)= lhs-rhs;
lhs =y(86);
rhs =T379;
residual(45)= lhs-rhs;
lhs =y(87);
rhs =T271;
residual(46)= lhs-rhs;
lhs =y(44);
rhs =T408*y(31)/y(32);
residual(47)= lhs-rhs;
lhs =y(44);
rhs =T417^(1/(1-params(40)));
residual(48)= lhs-rhs;
lhs =y(31);
rhs =y(50)*y(40)+T427*y(120);
residual(49)= lhs-rhs;
lhs =y(32);
rhs =y(50)+T433*y(121);
residual(50)= lhs-rhs;
lhs =y(57);
rhs =y(43)/y(53);
residual(51)= lhs-rhs;
lhs =y(57);
rhs =params(34)*y(62)*y(125)/params(54);
residual(52)= lhs-rhs;
lhs =y(61);
rhs =((1-params(35))*y(63)/y(10)*y(59)+(1-params(37))*y(60))/y(11);
residual(53)= lhs-rhs;
lhs =y(52);
rhs =y(62)/params(54);
residual(54)= lhs-rhs;
lhs =y(126);
rhs =y(52)*params(61);
residual(55)= lhs-rhs;
lhs =y(49)*y(60);
rhs =1+params(41)*(y(55)/y(10)-params(37));
residual(56)= lhs-rhs;
lhs =y(56);
rhs =(1-params(37))*y(10)+y(49)*y(55);
residual(57)= lhs-rhs;
lhs =y(63);
rhs =y(53)+y(55)+y(10)*T481;
residual(58)= lhs-rhs;
lhs =y(63)*y(58);
rhs =T489*T490;
residual(59)= lhs-rhs;
lhs =params(39)*y(54)^params(38);
rhs =y(59)*params(35)*y(57)*y(63)*y(58)/y(54);
residual(60)= lhs-rhs;
lhs =y(58);
rhs =1;
residual(61)= lhs-rhs;
lhs =y(64);
rhs =y(50)-y(63);
residual(62)= lhs-rhs;
lhs =y(67);
rhs =log(y(64))-log(params(95));
residual(63)= lhs-rhs;
lhs =y(84);
rhs =y(45)-y(60);
residual(64)= lhs-rhs;
lhs =y(59);
rhs =(params(40)-1)/params(40)*params(93);
residual(65)= lhs-rhs;
lhs =log(y(51));
rhs =params(47)*log(y(9))+x(it_, 1);
residual(66)= lhs-rhs;
lhs =log(y(49));
rhs =params(49)*log(y(8))+x(it_, 2);
residual(67)= lhs-rhs;
lhs =log(y(43));
rhs =params(52)*log(y(5))+x(it_, 3);
residual(68)= lhs-rhs;
lhs =log(y(68));
rhs =params(103)*log(y(12))-x(it_, 4);
residual(69)= lhs-rhs;
lhs =y(99);
rhs =y(128);
residual(70)= lhs-rhs;
lhs =y(100);
rhs =y(131);
residual(71)= lhs-rhs;
lhs =y(101);
rhs =y(132);
residual(72)= lhs-rhs;
lhs =y(102);
rhs =y(129);
residual(73)= lhs-rhs;
lhs =y(103);
rhs =y(134);
residual(74)= lhs-rhs;
lhs =y(104);
rhs =y(135);
residual(75)= lhs-rhs;
lhs =y(105);
rhs =y(130);
residual(76)= lhs-rhs;
lhs =y(106);
rhs =y(137);
residual(77)= lhs-rhs;
lhs =y(107);
rhs =y(138);
residual(78)= lhs-rhs;
lhs =y(108);
rhs =y(15);
residual(79)= lhs-rhs;
lhs =y(109);
rhs =y(18);
residual(80)= lhs-rhs;
lhs =y(110);
rhs =y(17);
residual(81)= lhs-rhs;
lhs =y(111);
rhs =y(20);
residual(82)= lhs-rhs;
lhs =y(112);
rhs =y(21);
residual(83)= lhs-rhs;
lhs =y(113);
rhs =y(16);
residual(84)= lhs-rhs;
lhs =y(114);
rhs =y(23);
residual(85)= lhs-rhs;
lhs =y(115);
rhs =y(24);
residual(86)= lhs-rhs;
lhs =y(116);
rhs =y(14);
residual(87)= lhs-rhs;
lhs =y(117);
rhs =y(26);
residual(88)= lhs-rhs;
lhs =y(118);
rhs =y(27);
residual(89)= lhs-rhs;
lhs =y(119);
rhs =y(19);
residual(90)= lhs-rhs;
if nargout >= 2,
  g1 = zeros(90, 144);

  %
  % Jacobian matrix
  %

T3 = (-1);
T713 = getPowerDeriv(y(51)*y(35),params(35),1);
T718 = params(39)*getPowerDeriv(y(35),params(38),1);
T737 = 1/y(2);
T744 = 2*(y(36)/y(2)-params(37));
T747 = y(2)*T268*T737*T744;
T767 = (-y(36))/(y(2)*y(2));
T772 = T268*T744*T767;
T775 = (-(T270+y(2)*T772));
T776 = getPowerDeriv(y(2),1-params(35),1);
T780 = (-(y(39)*y(50)*(1-params(35))))/(y(2)*y(2));
T782 = (-((y(46)*y(6)-(1-params(37))*y(45))*T780));
T797 = y(45)/y(41);
T798 = getPowerDeriv(T379,params(36),1);
T815 = params(34)*params(42)*(-y(122))/(y(38)*y(38));
T826 = params(34)*params(42)*1/y(38);
T847 = y(50)*(1-params(35))/y(2);
T849 = (-((y(46)*y(6)-(1-params(37))*y(45))*T847));
T864 = (-(y(45)*y(37)))/(y(41)*y(41));
T871 = 1/params(54);
T876 = params(42)*getPowerDeriv(y(42),params(40),1);
T883 = (-(params(42)*getPowerDeriv(y(42),params(40)-1,1)))/(1-params(42));
T884 = getPowerDeriv(T417,1/(1-params(40)),1);
T894 = getPowerDeriv(y(123),params(40),1);
T895 = T425*T894;
T898 = getPowerDeriv(y(123),params(40)-1,1);
T899 = T425*T898;
T926 = y(37)/y(41);
T938 = getPowerDeriv(T360,1-params(35),1);
T959 = getPowerDeriv(y(49),params(51),1);
T973 = (-(params(43)*2*T959))/(2*T254*2*T254);
T982 = 1/y(63);
T1003 = (1-params(35))*y(39)/y(2);
T1005 = (-((y(46)*y(6)-(1-params(37))*y(45))*T1003));
T1017 = getPowerDeriv(y(51),params(35),1);
T1023 = getPowerDeriv(y(51)*y(54),params(35),1);
T1043 = 1/y(10);
T1046 = 2*(y(55)/y(10)-params(37));
T1059 = (-y(55))/(y(10)*y(10));
T1063 = T268*T1046*T1059;
T1067 = getPowerDeriv(y(10),1-params(35),1);
T1092 = (-y(50))/(y(63)*y(63));
T1141 = (-((1-params(35))*y(46)*params(35)))/((1-params(35))*y(76)*(1-params(35))*y(76));
T1142 = T938*T1141;
  g1(1,47)=(-(400*1/params(59)/T11));
  g1(1,94)=1;
  g1(2,42)=(-(400*T871/T19));
  g1(2,96)=1;
  g1(3,95)=1;
  g1(3,15)=(-0.25);
  g1(3,96)=(-0.25);
  g1(3,18)=(-0.25);
  g1(3,19)=(-0.25);
  g1(4,50)=(-(100*1/params(72)/T40));
  g1(4,98)=1;
  g1(5,50)=(-(100*T982/T47));
  g1(5,63)=(-(100*T1092/T47));
  g1(5,97)=1;
  g1(6,93)=1;
  g1(6,95)=(-(0.5*2*(y(95)-y(96))));
  g1(6,96)=(-(0.5*(-(2*(y(95)-y(96))))));
  g1(6,97)=(-(0.5*0.01*2*y(97)));
  g1(7,14)=(-params(1));
  g1(7,94)=1;
  g1(7,15)=(-params(6));
  g1(7,96)=(-params(5));
  g1(7,128)=(-params(10));
  g1(7,16)=(-params(15));
  g1(7,97)=(-params(14));
  g1(7,129)=(-params(19));
  g1(7,17)=(-params(24));
  g1(7,98)=(-params(23));
  g1(7,130)=(-params(28));
  g1(7,144)=(-params(32));
  g1(7,131)=(-params(11));
  g1(7,132)=(-params(12));
  g1(7,133)=(-params(13));
  g1(7,134)=(-params(20));
  g1(7,135)=(-params(21));
  g1(7,136)=(-params(22));
  g1(7,137)=(-params(29));
  g1(7,138)=(-params(30));
  g1(7,139)=(-params(31));
  g1(7,18)=(-params(7));
  g1(7,19)=(-params(8));
  g1(7,20)=(-params(25));
  g1(7,21)=(-params(26));
  g1(7,22)=(-params(27));
  g1(7,23)=(-params(16));
  g1(7,24)=(-params(17));
  g1(7,25)=(-params(18));
  g1(7,26)=(-params(2));
  g1(7,27)=(-params(3));
  g1(7,28)=(-params(4));
  g1(7,29)=(-params(9));
  g1(8,33)=(-((-y(43))/(y(33)*y(33))));
  g1(8,38)=1;
  g1(8,43)=(-(1/y(33)));
  g1(9,38)=1;
  g1(9,122)=(-(y(47)*params(34)*1/y(123)));
  g1(9,123)=(-(y(47)*params(34)*(-y(122))/(y(123)*y(123))));
  g1(9,47)=(-(params(34)*y(122)/y(123)));
  g1(10,2)=(-(y(40)*(1-params(35))*(-(y(50)*y(39)))/(y(2)*y(2))/y(6)));
  g1(10,39)=(-(y(40)*(1-params(35))*y(50)/y(2)/y(6)));
  g1(10,40)=(-(T203/y(6)));
  g1(10,6)=(-((-(T203*y(40)+(1-params(37))*y(45)))/(y(6)*y(6))));
  g1(10,45)=(-((1-params(37))/y(6)));
  g1(10,46)=1;
  g1(10,50)=(-(y(40)*(1-params(35))*y(39)/y(2)/y(6)));
  g1(11,30)=1;
  g1(11,123)=(-((-y(47))/(y(123)*y(123))));
  g1(11,47)=(-(1/y(123)));
  g1(12,123)=T3;
  g1(12,85)=1;
  g1(13,30)=(-y(48));
  g1(13,124)=1;
  g1(13,48)=(-y(30));
  g1(14,30)=(-y(48));
  g1(14,48)=(-y(30));
  g1(14,90)=1;
  g1(15,1)=(-(params(96)*(-(y(7)*(y(2)*y(6)-y(4))))));
  g1(15,2)=(-(params(96)*(y(46)*y(6)-y(6)*y(7)*y(1))));
  g1(15,4)=(-(params(96)*y(7)*y(1)));
  g1(15,41)=1;
  g1(15,6)=(-(params(96)*(y(46)*y(2)-y(2)*y(7)*y(1))));
  g1(15,46)=(-(params(96)*y(2)*y(6)));
  g1(15,7)=(-(params(96)*(-(y(1)*(y(2)*y(6)-y(4))))));
  g1(16,41)=(-(1-params(96)));
  g1(16,88)=1;
  g1(17,36)=(-(params(41)*T737+params(43)/T254));
  g1(17,2)=(-(params(41)*T767));
  g1(17,45)=y(49);
  g1(17,49)=y(45)-(-(y(36)*params(43)*T959))/(T254*T254);
  g1(18,36)=(-y(49));
  g1(18,2)=(-(1-params(37)));
  g1(18,37)=1;
  g1(18,49)=(-y(36));
  g1(19,34)=1;
  g1(19,36)=(-T273);
  g1(19,49)=(-((-(y(36)*params(43)*2*T959))/(2*T254*2*T254)));
  g1(20,33)=T3;
  g1(20,36)=(-(1+T747+T273*2*y(36)));
  g1(20,2)=T775;
  g1(20,49)=(-(T274*T973));
  g1(20,50)=1;
  g1(21,35)=(-(T282*y(51)*T713));
  g1(21,2)=(-(T281*T776));
  g1(21,39)=y(50);
  g1(21,50)=y(39);
  g1(21,51)=(-(T282*y(35)*T713));
  g1(22,35)=T718-y(40)*y(38)*params(35)*(-(y(50)*y(39)))/(y(35)*y(35));
  g1(22,38)=(-(y(40)*params(35)*y(50)*y(39)/y(35)));
  g1(22,39)=(-(y(40)*y(38)*params(35)*y(50)/y(35)));
  g1(22,40)=(-(y(38)*params(35)*y(50)*y(39)/y(35)));
  g1(22,50)=(-(y(40)*y(38)*params(35)*y(39)/y(35)));
  g1(23,3)=(-T302);
  g1(23,39)=1;
  g1(23,42)=(-(y(3)*T876));
  g1(23,44)=(-((1-params(42))*getPowerDeriv(y(44),(-params(40)),1)));
  g1(24,50)=(-(1/y(50)));
  g1(24,65)=1;
  g1(25,36)=(-(1/y(36)));
  g1(25,66)=1;
  g1(26,33)=(-(1/y(33)));
  g1(26,35)=params(39)*getPowerDeriv(y(35),1+params(38),1)/(1+params(38));
  g1(26,91)=1;
  g1(27,91)=T3;
  g1(27,92)=1;
  g1(27,127)=(-params(34));
  g1(28,35)=(-((-(y(39)*y(50)*params(35)))/(y(35)*y(35))));
  g1(28,39)=(-(y(50)*params(35)/y(35)));
  g1(28,50)=(-(params(35)*y(39)/y(35)));
  g1(28,73)=1;
  g1(29,37)=(-((-(y(39)*y(50)*(1-params(35))))/(y(37)*y(37))));
  g1(29,39)=(-(y(50)*(1-params(35))/y(37)));
  g1(29,50)=(-((1-params(35))*y(39)/y(37)));
  g1(29,75)=1;
  g1(30,35)=(-((-y(37))/(y(35)*y(35))));
  g1(30,37)=(-(1/y(35)));
  g1(30,74)=1;
  g1(31,35)=(-(T718/y(38)));
  g1(31,38)=(-((-T288)/(y(38)*y(38))));
  g1(31,76)=1;
  g1(32,73)=(-((-y(76))/(y(73)*y(73))));
  g1(32,76)=(-(1/y(73)));
  g1(32,77)=1;
  g1(33,46)=(-(1/y(75)));
  g1(33,75)=(-((-y(46))/(y(75)*y(75))));
  g1(33,78)=1;
  g1(34,77)=T3;
  g1(34,78)=T3;
  g1(34,79)=1;
  g1(35,46)=(-(T357*params(35)/((1-params(35))*y(76))*T938));
  g1(35,76)=(-(T361*1/params(35)+T357*T1142));
  g1(35,80)=1;
  g1(36,73)=(-((-y(76))/(y(73)*y(73))));
  g1(36,76)=(-(1/y(73)));
  g1(36,81)=1;
  g1(37,2)=(-(T782/(T368*T368)));
  g1(37,39)=(-(T849/(T368*T368)));
  g1(37,6)=(-(y(46)/T368));
  g1(37,45)=(-((-(1-params(37)))/T368));
  g1(37,46)=(-(y(6)/T368));
  g1(37,50)=(-(T1005/(T368*T368)));
  g1(37,82)=1;
  g1(38,2)=(-((-(y(50)*y(39)*T373*T776))/(T282*T373*T282*T373)));
  g1(38,39)=(-(y(50)/(T282*T373)));
  g1(38,50)=(-(y(39)/(T282*T373)));
  g1(38,51)=(-((-(y(50)*y(39)*T282*T1017))/(T282*T373*T282*T373)));
  g1(38,83)=getPowerDeriv(y(83),params(35),1);
  g1(39,37)=(-(y(69)*y(68)*T797*T798));
  g1(39,41)=(-(y(69)*y(68)*T798*T864));
  g1(39,45)=(-(y(69)*y(68)*T798*T926));
  g1(39,48)=1;
  g1(39,68)=(-(T381*y(69)));
  g1(39,69)=(-(y(68)*T381));
  g1(40,37)=(-(T797*T798));
  g1(40,41)=(-(T798*T864));
  g1(40,45)=(-(T798*T926));
  g1(40,89)=1;
  g1(41,37)=(-y(45));
  g1(41,41)=1;
  g1(41,45)=(-y(37));
  g1(41,70)=1;
  g1(42,69)=1;
  g1(42,71)=(-(1/params(108)*getPowerDeriv(y(71)/params(108),params(105),1)));
  g1(43,42)=(-(y(70)/y(13)));
  g1(43,13)=(-(y(42)*(-y(70))/(y(13)*y(13))));
  g1(43,70)=(-(y(42)*1/y(13)));
  g1(43,71)=1;
  g1(44,13)=(-((-y(70))/(y(13)*y(13))));
  g1(44,70)=(-(1/y(13)));
  g1(44,72)=1;
  g1(45,37)=(-T797);
  g1(45,41)=(-T864);
  g1(45,45)=(-T926);
  g1(45,86)=1;
  g1(46,36)=(-T747);
  g1(46,2)=T775;
  g1(46,87)=1;
  g1(47,31)=(-(T408*1/y(32)));
  g1(47,32)=(-(T408*(-y(31))/(y(32)*y(32))));
  g1(47,44)=1;
  g1(48,42)=(-(T883*T884));
  g1(48,44)=1;
  g1(49,31)=1;
  g1(49,120)=(-T427);
  g1(49,38)=(-(y(120)*T426*T815));
  g1(49,122)=(-(y(120)*T426*T826));
  g1(49,40)=(-y(50));
  g1(49,123)=(-(y(120)*T895));
  g1(49,50)=(-y(40));
  g1(50,32)=1;
  g1(50,121)=(-T433);
  g1(50,38)=(-(y(121)*T432*T815));
  g1(50,122)=(-(y(121)*T432*T826));
  g1(50,123)=(-(y(121)*T899));
  g1(50,50)=T3;
  g1(51,43)=(-(1/y(53)));
  g1(51,53)=(-((-y(43))/(y(53)*y(53))));
  g1(51,57)=1;
  g1(52,57)=1;
  g1(52,125)=(-(params(34)*y(62)*T871));
  g1(52,62)=(-(params(34)*y(125)/params(54)));
  g1(53,10)=(-(y(59)*(1-params(35))*(-y(63))/(y(10)*y(10))/y(11)));
  g1(53,59)=(-((1-params(35))*y(63)/y(10)/y(11)));
  g1(53,11)=(-((-((1-params(35))*y(63)/y(10)*y(59)+(1-params(37))*y(60)))/(y(11)*y(11))));
  g1(53,60)=(-((1-params(37))/y(11)));
  g1(53,61)=1;
  g1(53,63)=(-(y(59)*(1-params(35))*T1043/y(11)));
  g1(54,52)=1;
  g1(54,62)=(-T871);
  g1(55,52)=(-params(61));
  g1(55,126)=1;
  g1(56,49)=y(60);
  g1(56,55)=(-(params(41)*T1043));
  g1(56,10)=(-(params(41)*T1059));
  g1(56,60)=y(49);
  g1(57,49)=(-y(55));
  g1(57,55)=(-y(49));
  g1(57,10)=(-(1-params(37)));
  g1(57,56)=1;
  g1(58,53)=T3;
  g1(58,55)=(-(1+y(10)*T268*T1043*T1046));
  g1(58,10)=(-(T481+y(10)*T1063));
  g1(58,63)=1;
  g1(59,51)=(-(T490*y(54)*T1023));
  g1(59,54)=(-(T490*y(51)*T1023));
  g1(59,10)=(-(T489*T1067));
  g1(59,58)=y(63);
  g1(59,63)=y(58);
  g1(60,54)=params(39)*getPowerDeriv(y(54),params(38),1)-y(59)*params(35)*y(57)*(-(y(63)*y(58)))/(y(54)*y(54));
  g1(60,57)=(-(y(59)*params(35)*y(63)*y(58)/y(54)));
  g1(60,58)=(-(y(59)*params(35)*y(57)*y(63)/y(54)));
  g1(60,59)=(-(params(35)*y(57)*y(63)*y(58)/y(54)));
  g1(60,63)=(-(y(59)*params(35)*y(57)*y(58)/y(54)));
  g1(61,58)=1;
  g1(62,50)=T3;
  g1(62,63)=1;
  g1(62,64)=1;
  g1(63,64)=(-(1/y(64)));
  g1(63,67)=1;
  g1(64,45)=T3;
  g1(64,60)=1;
  g1(64,84)=1;
  g1(65,59)=1;
  g1(66,9)=(-(params(47)*1/y(9)));
  g1(66,51)=1/y(51);
  g1(66,140)=T3;
  g1(67,8)=(-(params(49)*1/y(8)));
  g1(67,49)=1/y(49);
  g1(67,141)=T3;
  g1(68,5)=(-(params(52)*1/y(5)));
  g1(68,43)=1/y(43);
  g1(68,142)=T3;
  g1(69,12)=(-(params(103)*1/y(12)));
  g1(69,68)=1/y(68);
  g1(69,143)=1;
  g1(70,128)=T3;
  g1(70,99)=1;
  g1(71,131)=T3;
  g1(71,100)=1;
  g1(72,132)=T3;
  g1(72,101)=1;
  g1(73,129)=T3;
  g1(73,102)=1;
  g1(74,134)=T3;
  g1(74,103)=1;
  g1(75,135)=T3;
  g1(75,104)=1;
  g1(76,130)=T3;
  g1(76,105)=1;
  g1(77,137)=T3;
  g1(77,106)=1;
  g1(78,138)=T3;
  g1(78,107)=1;
  g1(79,15)=T3;
  g1(79,108)=1;
  g1(80,18)=T3;
  g1(80,109)=1;
  g1(81,17)=T3;
  g1(81,110)=1;
  g1(82,20)=T3;
  g1(82,111)=1;
  g1(83,21)=T3;
  g1(83,112)=1;
  g1(84,16)=T3;
  g1(84,113)=1;
  g1(85,23)=T3;
  g1(85,114)=1;
  g1(86,24)=T3;
  g1(86,115)=1;
  g1(87,14)=T3;
  g1(87,116)=1;
  g1(88,26)=T3;
  g1(88,117)=1;
  g1(89,27)=T3;
  g1(89,118)=1;
  g1(90,19)=T3;
  g1(90,119)=1;

if nargout >= 3,
  %
  % Hessian matrix
  %

  v2 = zeros(374,3);
T1356 = getPowerDeriv(y(49),params(51),2);
T1392 = (-(T268*T737*T744+y(2)*T268*(T767*2*T737+T744*T3/(y(2)*y(2)))));
T1401 = (-(T772+T772+y(2)*T268*(T767*2*T767+T744*(-((-y(36))*(y(2)+y(2))))/(y(2)*y(2)*y(2)*y(2)))));
T1412 = getPowerDeriv(y(51)*y(35),params(35),2);
T1413 = y(51)*T1412;
T1419 = getPowerDeriv(y(2),1-params(35),2);
T1433 = params(39)*getPowerDeriv(y(35),params(38),2);
T1542 = getPowerDeriv(T360,1-params(35),2);
T1543 = params(35)/((1-params(35))*y(76))*T1542;
T1578 = T368*T780+T368*T780;
T1708 = getPowerDeriv(T379,params(36),2);
T1709 = T797*T1708;
T1729 = T864*T864*T1708+T798*(-((-(y(45)*y(37)))*(y(41)+y(41))))/(y(41)*y(41)*y(41)*y(41));
T1926 = getPowerDeriv(y(51)*y(54),params(35),2);
T1927 = y(54)*T1926;
  v2(1,1)=1;
  v2(1,2)=6671;
  v2(1,3)=(-(400*(-(1/params(59)*1/params(59)))/(T11*T11)));
  v2(2,1)=2;
  v2(2,2)=5946;
  v2(2,3)=(-(400*(-(T871*T871))/(T19*T19)));
  v2(3,1)=4;
  v2(3,2)=7106;
  v2(3,3)=(-(100*(-(1/params(72)*1/params(72)))/(T40*T40)));
  v2(4,1)=5;
  v2(4,2)=7106;
  v2(4,3)=(-(100*(-(T982*T982))/(T47*T47)));
  v2(5,1)=5;
  v2(5,2)=8978;
  v2(5,3)=(-(100*(T47*T3/(y(63)*y(63))-T982*T1092)/(T47*T47)));
  v2(6,1)=5;
  v2(6,2)=7119;
  v2(6,3)=  v2(5,3);
  v2(7,1)=5;
  v2(7,2)=8991;
  v2(7,3)=(-(100*(T47*(-((-y(50))*(y(63)+y(63))))/(y(63)*y(63)*y(63)*y(63))-T1092*T1092)/(T47*T47)));
  v2(8,1)=6;
  v2(8,2)=13631;
  v2(8,3)=T3;
  v2(9,1)=6;
  v2(9,2)=13775;
  v2(9,3)=1;
  v2(10,1)=6;
  v2(10,2)=13632;
  v2(10,3)=  v2(9,3);
  v2(11,1)=6;
  v2(11,2)=13776;
  v2(11,3)=T3;
  v2(12,1)=6;
  v2(12,2)=13921;
  v2(12,3)=(-0.01);
  v2(13,1)=8;
  v2(13,2)=4641;
  v2(13,3)=(-((-((-y(43))*(y(33)+y(33))))/(y(33)*y(33)*y(33)*y(33))));
  v2(14,1)=8;
  v2(14,2)=6081;
  v2(14,3)=(-(T3/(y(33)*y(33))));
  v2(15,1)=8;
  v2(15,2)=4651;
  v2(15,3)=  v2(14,3);
  v2(16,1)=9;
  v2(16,2)=17690;
  v2(16,3)=(-(y(47)*params(34)*T3/(y(123)*y(123))));
  v2(17,1)=9;
  v2(17,2)=17547;
  v2(17,3)=  v2(16,3);
  v2(18,1)=9;
  v2(18,2)=17691;
  v2(18,3)=(-(y(47)*params(34)*(-((-y(122))*(y(123)+y(123))))/(y(123)*y(123)*y(123)*y(123))));
  v2(19,1)=9;
  v2(19,2)=6746;
  v2(19,3)=(-(params(34)*1/y(123)));
  v2(20,1)=9;
  v2(20,2)=17471;
  v2(20,3)=  v2(19,3);
  v2(21,1)=9;
  v2(21,2)=6747;
  v2(21,3)=(-(params(34)*(-y(122))/(y(123)*y(123))));
  v2(22,1)=9;
  v2(22,2)=17615;
  v2(22,3)=  v2(21,3);
  v2(23,1)=10;
  v2(23,2)=146;
  v2(23,3)=(-(y(40)*(1-params(35))*(-((-(y(50)*y(39)))*(y(2)+y(2))))/(y(2)*y(2)*y(2)*y(2))/y(6)));
  v2(24,1)=10;
  v2(24,2)=5474;
  v2(24,3)=(-(y(40)*(1-params(35))*(-y(50))/(y(2)*y(2))/y(6)));
  v2(25,1)=10;
  v2(25,2)=183;
  v2(25,3)=  v2(24,3);
  v2(26,1)=10;
  v2(26,2)=5618;
  v2(26,3)=(-((1-params(35))*(-(y(50)*y(39)))/(y(2)*y(2))/y(6)));
  v2(27,1)=10;
  v2(27,2)=184;
  v2(27,3)=  v2(26,3);
  v2(28,1)=10;
  v2(28,2)=5655;
  v2(28,3)=(-((1-params(35))*y(50)/y(2)/y(6)));
  v2(29,1)=10;
  v2(29,2)=5512;
  v2(29,3)=  v2(28,3);
  v2(30,1)=10;
  v2(30,2)=722;
  v2(30,3)=(-((-(y(40)*(1-params(35))*(-(y(50)*y(39)))/(y(2)*y(2))))/(y(6)*y(6))));
  v2(31,1)=10;
  v2(31,2)=150;
  v2(31,3)=  v2(30,3);
  v2(32,1)=10;
  v2(32,2)=759;
  v2(32,3)=(-((-(y(40)*(1-params(35))*y(50)/y(2)))/(y(6)*y(6))));
  v2(33,1)=10;
  v2(33,2)=5478;
  v2(33,3)=  v2(32,3);
  v2(34,1)=10;
  v2(34,2)=760;
  v2(34,3)=(-((-T203)/(y(6)*y(6))));
  v2(35,1)=10;
  v2(35,2)=5622;
  v2(35,3)=  v2(34,3);
  v2(36,1)=10;
  v2(36,2)=726;
  v2(36,3)=(-((-((-(T203*y(40)+(1-params(37))*y(45)))*(y(6)+y(6))))/(y(6)*y(6)*y(6)*y(6))));
  v2(37,1)=10;
  v2(37,2)=6342;
  v2(37,3)=(-((-(1-params(37)))/(y(6)*y(6))));
  v2(38,1)=10;
  v2(38,2)=765;
  v2(38,3)=  v2(37,3);
  v2(39,1)=10;
  v2(39,2)=7058;
  v2(39,3)=(-(y(40)*(1-params(35))*(-y(39))/(y(2)*y(2))/y(6)));
  v2(40,1)=10;
  v2(40,2)=194;
  v2(40,3)=  v2(39,3);
  v2(41,1)=10;
  v2(41,2)=7095;
  v2(41,3)=(-(y(40)*(1-params(35))*T737/y(6)));
  v2(42,1)=10;
  v2(42,2)=5522;
  v2(42,3)=  v2(41,3);
  v2(43,1)=10;
  v2(43,2)=7096;
  v2(43,3)=(-((1-params(35))*y(39)/y(2)/y(6)));
  v2(44,1)=10;
  v2(44,2)=5666;
  v2(44,3)=  v2(43,3);
  v2(45,1)=10;
  v2(45,2)=7062;
  v2(45,3)=(-((-(y(40)*(1-params(35))*y(39)/y(2)))/(y(6)*y(6))));
  v2(46,1)=10;
  v2(46,2)=770;
  v2(46,3)=  v2(45,3);
  v2(47,1)=11;
  v2(47,2)=17691;
  v2(47,3)=(-((-((-y(47))*(y(123)+y(123))))/(y(123)*y(123)*y(123)*y(123))));
  v2(48,1)=11;
  v2(48,2)=6747;
  v2(48,3)=(-(T3/(y(123)*y(123))));
  v2(49,1)=11;
  v2(49,2)=17615;
  v2(49,3)=  v2(48,3);
  v2(50,1)=13;
  v2(50,2)=6798;
  v2(50,3)=T3;
  v2(51,1)=13;
  v2(51,2)=4224;
  v2(51,3)=  v2(50,3);
  v2(52,1)=14;
  v2(52,2)=6798;
  v2(52,3)=T3;
  v2(53,1)=14;
  v2(53,2)=4224;
  v2(53,3)=  v2(52,3);
  v2(54,1)=15;
  v2(54,2)=145;
  v2(54,3)=(-(params(96)*(-(y(6)*y(7)))));
  v2(55,1)=15;
  v2(55,2)=2;
  v2(55,3)=  v2(54,3);
  v2(56,1)=15;
  v2(56,2)=433;
  v2(56,3)=(-(params(96)*y(7)));
  v2(57,1)=15;
  v2(57,2)=4;
  v2(57,3)=  v2(56,3);
  v2(58,1)=15;
  v2(58,2)=721;
  v2(58,3)=(-(params(96)*(-(y(2)*y(7)))));
  v2(59,1)=15;
  v2(59,2)=6;
  v2(59,3)=  v2(58,3);
  v2(60,1)=15;
  v2(60,2)=722;
  v2(60,3)=(-(params(96)*(y(46)-y(7)*y(1))));
  v2(61,1)=15;
  v2(61,2)=150;
  v2(61,3)=  v2(60,3);
  v2(62,1)=15;
  v2(62,2)=6482;
  v2(62,3)=(-(y(6)*params(96)));
  v2(63,1)=15;
  v2(63,2)=190;
  v2(63,3)=  v2(62,3);
  v2(64,1)=15;
  v2(64,2)=6486;
  v2(64,3)=(-(y(2)*params(96)));
  v2(65,1)=15;
  v2(65,2)=766;
  v2(65,3)=  v2(64,3);
  v2(66,1)=15;
  v2(66,2)=865;
  v2(66,3)=(-(params(96)*(-(y(2)*y(6)-y(4)))));
  v2(67,1)=15;
  v2(67,2)=7;
  v2(67,3)=  v2(66,3);
  v2(68,1)=15;
  v2(68,2)=866;
  v2(68,3)=(-(params(96)*(-(y(6)*y(1)))));
  v2(69,1)=15;
  v2(69,2)=151;
  v2(69,3)=  v2(68,3);
  v2(70,1)=15;
  v2(70,2)=868;
  v2(70,3)=(-(params(96)*y(1)));
  v2(71,1)=15;
  v2(71,2)=439;
  v2(71,3)=  v2(70,3);
  v2(72,1)=15;
  v2(72,2)=870;
  v2(72,3)=(-(params(96)*(-(y(2)*y(1)))));
  v2(73,1)=15;
  v2(73,2)=727;
  v2(73,3)=  v2(72,3);
  v2(74,1)=17;
  v2(74,2)=180;
  v2(74,3)=(-(params(41)*T3/(y(2)*y(2))));
  v2(75,1)=17;
  v2(75,2)=5042;
  v2(75,3)=  v2(74,3);
  v2(76,1)=17;
  v2(76,2)=146;
  v2(76,3)=(-(params(41)*(-((-y(36))*(y(2)+y(2))))/(y(2)*y(2)*y(2)*y(2))));
  v2(77,1)=17;
  v2(77,2)=6948;
  v2(77,3)=(-((-(params(43)*T959))/(T254*T254)));
  v2(78,1)=17;
  v2(78,2)=5089;
  v2(78,3)=  v2(77,3);
  v2(79,1)=17;
  v2(79,2)=6957;
  v2(79,3)=1;
  v2(80,1)=17;
  v2(80,2)=6385;
  v2(80,3)=  v2(79,3);
  v2(81,1)=17;
  v2(81,2)=6961;
  v2(81,3)=(-((T254*T254*(-(y(36)*params(43)*T1356))-(-(y(36)*params(43)*T959))*(T254*T959+T254*T959))/(T254*T254*T254*T254)));
  v2(82,1)=18;
  v2(82,2)=6948;
  v2(82,3)=T3;
  v2(83,1)=18;
  v2(83,2)=5089;
  v2(83,3)=  v2(82,3);
  v2(84,1)=19;
  v2(84,2)=6948;
  v2(84,3)=(-T973);
  v2(85,1)=19;
  v2(85,2)=5089;
  v2(85,3)=  v2(84,3);
  v2(86,1)=19;
  v2(86,2)=6961;
  v2(86,3)=(-((2*T254*2*T254*(-(y(36)*params(43)*2*T1356))-(-(y(36)*params(43)*2*T959))*(2*T254*2*T959+2*T254*2*T959))/(2*T254*2*T254*2*T254*2*T254)));
  v2(87,1)=20;
  v2(87,2)=5076;
  v2(87,3)=(-(y(2)*T268*T737*2*T737+2*T273));
  v2(88,1)=20;
  v2(88,2)=180;
  v2(88,3)=T1392;
  v2(89,1)=20;
  v2(89,2)=5042;
  v2(89,3)=  v2(88,3);
  v2(90,1)=20;
  v2(90,2)=146;
  v2(90,3)=T1401;
  v2(91,1)=20;
  v2(91,2)=6948;
  v2(91,3)=(-(2*y(36)*T973));
  v2(92,1)=20;
  v2(92,2)=5089;
  v2(92,3)=  v2(91,3);
  v2(93,1)=20;
  v2(93,2)=6961;
  v2(93,3)=(-(T274*(2*T254*2*T254*(-(params(43)*2*T1356))-(-(params(43)*2*T959))*(2*T254*2*T959+2*T254*2*T959))/(2*T254*2*T254*2*T254*2*T254)));
  v2(94,1)=21;
  v2(94,2)=4931;
  v2(94,3)=(-(T282*y(51)*T1413));
  v2(95,1)=21;
  v2(95,2)=179;
  v2(95,3)=(-(y(51)*T713*T776));
  v2(96,1)=21;
  v2(96,2)=4898;
  v2(96,3)=  v2(95,3);
  v2(97,1)=21;
  v2(97,2)=146;
  v2(97,3)=(-(T281*T1419));
  v2(98,1)=21;
  v2(98,2)=7095;
  v2(98,3)=1;
  v2(99,1)=21;
  v2(99,2)=5522;
  v2(99,3)=  v2(98,3);
  v2(100,1)=21;
  v2(100,2)=7235;
  v2(100,3)=(-(T282*(T713+y(35)*T1413)));
  v2(101,1)=21;
  v2(101,2)=4947;
  v2(101,3)=  v2(100,3);
  v2(102,1)=21;
  v2(102,2)=7202;
  v2(102,3)=(-(T776*y(35)*T713));
  v2(103,1)=21;
  v2(103,2)=195;
  v2(103,3)=  v2(102,3);
  v2(104,1)=21;
  v2(104,2)=7251;
  v2(104,3)=(-(T282*y(35)*y(35)*T1412));
  v2(105,1)=22;
  v2(105,2)=4931;
  v2(105,3)=T1433-y(40)*y(38)*params(35)*(-((-(y(50)*y(39)))*(y(35)+y(35))))/(y(35)*y(35)*y(35)*y(35));
  v2(106,1)=22;
  v2(106,2)=5363;
  v2(106,3)=(-(y(40)*params(35)*(-(y(50)*y(39)))/(y(35)*y(35))));
  v2(107,1)=22;
  v2(107,2)=4934;
  v2(107,3)=  v2(106,3);
  v2(108,1)=22;
  v2(108,2)=5507;
  v2(108,3)=(-(y(40)*y(38)*params(35)*(-y(50))/(y(35)*y(35))));
  v2(109,1)=22;
  v2(109,2)=4935;
  v2(109,3)=  v2(108,3);
  v2(110,1)=22;
  v2(110,2)=5510;
  v2(110,3)=(-(y(40)*params(35)*y(50)/y(35)));
  v2(111,1)=22;
  v2(111,2)=5367;
  v2(111,3)=  v2(110,3);
  v2(112,1)=22;
  v2(112,2)=5651;
  v2(112,3)=(-(y(38)*params(35)*(-(y(50)*y(39)))/(y(35)*y(35))));
  v2(113,1)=22;
  v2(113,2)=4936;
  v2(113,3)=  v2(112,3);
  v2(114,1)=22;
  v2(114,2)=5654;
  v2(114,3)=(-(params(35)*y(50)*y(39)/y(35)));
  v2(115,1)=22;
  v2(115,2)=5368;
  v2(115,3)=  v2(114,3);
  v2(116,1)=22;
  v2(116,2)=5655;
  v2(116,3)=(-(y(38)*params(35)*y(50)/y(35)));
  v2(117,1)=22;
  v2(117,2)=5512;
  v2(117,3)=  v2(116,3);
  v2(118,1)=22;
  v2(118,2)=7091;
  v2(118,3)=(-(y(40)*y(38)*params(35)*(-y(39))/(y(35)*y(35))));
  v2(119,1)=22;
  v2(119,2)=4946;
  v2(119,3)=  v2(118,3);
  v2(120,1)=22;
  v2(120,2)=7094;
  v2(120,3)=(-(y(40)*params(35)*y(39)/y(35)));
  v2(121,1)=22;
  v2(121,2)=5378;
  v2(121,3)=  v2(120,3);
  v2(122,1)=22;
  v2(122,2)=7095;
  v2(122,3)=(-(y(40)*y(38)*params(35)*1/y(35)));
  v2(123,1)=22;
  v2(123,2)=5522;
  v2(123,3)=  v2(122,3);
  v2(124,1)=22;
  v2(124,2)=7096;
  v2(124,3)=(-(y(38)*params(35)*y(39)/y(35)));
  v2(125,1)=22;
  v2(125,2)=5666;
  v2(125,3)=  v2(124,3);
  v2(126,1)=23;
  v2(126,2)=5907;
  v2(126,3)=(-T876);
  v2(127,1)=23;
  v2(127,2)=330;
  v2(127,3)=  v2(126,3);
  v2(128,1)=23;
  v2(128,2)=5946;
  v2(128,3)=(-(y(3)*params(42)*getPowerDeriv(y(42),params(40),2)));
  v2(129,1)=23;
  v2(129,2)=6236;
  v2(129,3)=(-((1-params(42))*getPowerDeriv(y(44),(-params(40)),2)));
  v2(130,1)=24;
  v2(130,2)=7106;
  v2(130,3)=(-(T3/(y(50)*y(50))));
  v2(131,1)=25;
  v2(131,2)=5076;
  v2(131,3)=(-(T3/(y(36)*y(36))));
  v2(132,1)=26;
  v2(132,2)=4641;
  v2(132,3)=(-(T3/(y(33)*y(33))));
  v2(133,1)=26;
  v2(133,2)=4931;
  v2(133,3)=params(39)*getPowerDeriv(y(35),1+params(38),2)/(1+params(38));
  v2(134,1)=28;
  v2(134,2)=4931;
  v2(134,3)=(-((-((-(y(39)*y(50)*params(35)))*(y(35)+y(35))))/(y(35)*y(35)*y(35)*y(35))));
  v2(135,1)=28;
  v2(135,2)=5507;
  v2(135,3)=(-((-(y(50)*params(35)))/(y(35)*y(35))));
  v2(136,1)=28;
  v2(136,2)=4935;
  v2(136,3)=  v2(135,3);
  v2(137,1)=28;
  v2(137,2)=7091;
  v2(137,3)=(-((-(params(35)*y(39)))/(y(35)*y(35))));
  v2(138,1)=28;
  v2(138,2)=4946;
  v2(138,3)=  v2(137,3);
  v2(139,1)=28;
  v2(139,2)=7095;
  v2(139,3)=(-(params(35)/y(35)));
  v2(140,1)=28;
  v2(140,2)=5522;
  v2(140,3)=  v2(139,3);
  v2(141,1)=29;
  v2(141,2)=5221;
  v2(141,3)=(-((-((-(y(39)*y(50)*(1-params(35))))*(y(37)+y(37))))/(y(37)*y(37)*y(37)*y(37))));
  v2(142,1)=29;
  v2(142,2)=5509;
  v2(142,3)=(-((-(y(50)*(1-params(35))))/(y(37)*y(37))));
  v2(143,1)=29;
  v2(143,2)=5223;
  v2(143,3)=  v2(142,3);
  v2(144,1)=29;
  v2(144,2)=7093;
  v2(144,3)=(-((-((1-params(35))*y(39)))/(y(37)*y(37))));
  v2(145,1)=29;
  v2(145,2)=5234;
  v2(145,3)=  v2(144,3);
  v2(146,1)=29;
  v2(146,2)=7095;
  v2(146,3)=(-((1-params(35))/y(37)));
  v2(147,1)=29;
  v2(147,2)=5522;
  v2(147,3)=  v2(146,3);
  v2(148,1)=30;
  v2(148,2)=4931;
  v2(148,3)=(-((-((-y(37))*(y(35)+y(35))))/(y(35)*y(35)*y(35)*y(35))));
  v2(149,1)=30;
  v2(149,2)=5219;
  v2(149,3)=(-(T3/(y(35)*y(35))));
  v2(150,1)=30;
  v2(150,2)=4933;
  v2(150,3)=  v2(149,3);
  v2(151,1)=31;
  v2(151,2)=4931;
  v2(151,3)=(-(T1433/y(38)));
  v2(152,1)=31;
  v2(152,2)=5363;
  v2(152,3)=(-((-T718)/(y(38)*y(38))));
  v2(153,1)=31;
  v2(153,2)=4934;
  v2(153,3)=  v2(152,3);
  v2(154,1)=31;
  v2(154,2)=5366;
  v2(154,3)=(-((-((-T288)*(y(38)+y(38))))/(y(38)*y(38)*y(38)*y(38))));
  v2(155,1)=32;
  v2(155,2)=10441;
  v2(155,3)=(-((-((-y(76))*(y(73)+y(73))))/(y(73)*y(73)*y(73)*y(73))));
  v2(156,1)=32;
  v2(156,2)=10873;
  v2(156,3)=(-(T3/(y(73)*y(73))));
  v2(157,1)=32;
  v2(157,2)=10444;
  v2(157,3)=  v2(156,3);
  v2(158,1)=33;
  v2(158,2)=10702;
  v2(158,3)=(-(T3/(y(75)*y(75))));
  v2(159,1)=33;
  v2(159,2)=6555;
  v2(159,3)=  v2(158,3);
  v2(160,1)=33;
  v2(160,2)=10731;
  v2(160,3)=(-((-((-y(46))*(y(75)+y(75))))/(y(75)*y(75)*y(75)*y(75))));
  v2(161,1)=35;
  v2(161,2)=6526;
  v2(161,3)=(-(T357*params(35)/((1-params(35))*y(76))*T1543));
  v2(162,1)=35;
  v2(162,2)=10846;
  v2(162,3)=(-(params(35)/((1-params(35))*y(76))*T938*1/params(35)+T357*(T1141*T1543+T938*(-(params(35)*(1-params(35))))/((1-params(35))*y(76)*(1-params(35))*y(76)))));
  v2(163,1)=35;
  v2(163,2)=6556;
  v2(163,3)=  v2(162,3);
  v2(164,1)=35;
  v2(164,2)=10876;
  v2(164,3)=(-(1/params(35)*T1142+1/params(35)*T1142+T357*(T1141*T1141*T1542+T938*(-((-((1-params(35))*y(46)*params(35)))*((1-params(35))*(1-params(35))*y(76)+(1-params(35))*(1-params(35))*y(76))))/((1-params(35))*y(76)*(1-params(35))*y(76)*(1-params(35))*y(76)*(1-params(35))*y(76)))));
  v2(165,1)=36;
  v2(165,2)=10441;
  v2(165,3)=(-((-((-y(76))*(y(73)+y(73))))/(y(73)*y(73)*y(73)*y(73))));
  v2(166,1)=36;
  v2(166,2)=10873;
  v2(166,3)=(-(T3/(y(73)*y(73))));
  v2(167,1)=36;
  v2(167,2)=10444;
  v2(167,3)=  v2(166,3);
  v2(168,1)=37;
  v2(168,2)=146;
  v2(168,3)=(-((T368*T368*(-((y(46)*y(6)-(1-params(37))*y(45))*(-((-(y(39)*y(50)*(1-params(35))))*(y(2)+y(2))))/(y(2)*y(2)*y(2)*y(2))))-T782*T1578)/(T368*T368*T368*T368)));
  v2(169,1)=37;
  v2(169,2)=5474;
  v2(169,3)=(-((T368*T368*(-((y(46)*y(6)-(1-params(37))*y(45))*(-(y(50)*(1-params(35))))/(y(2)*y(2))))-T849*T1578)/(T368*T368*T368*T368)));
  v2(170,1)=37;
  v2(170,2)=183;
  v2(170,3)=  v2(169,3);
  v2(171,1)=37;
  v2(171,2)=5511;
  v2(171,3)=(-((-(T849*(T368*T847+T368*T847)))/(T368*T368*T368*T368)));
  v2(172,1)=37;
  v2(172,2)=722;
  v2(172,3)=(-((-(y(46)*T780))/(T368*T368)));
  v2(173,1)=37;
  v2(173,2)=150;
  v2(173,3)=  v2(172,3);
  v2(174,1)=37;
  v2(174,2)=759;
  v2(174,3)=(-((-(y(46)*T847))/(T368*T368)));
  v2(175,1)=37;
  v2(175,2)=5478;
  v2(175,3)=  v2(174,3);
  v2(176,1)=37;
  v2(176,2)=6338;
  v2(176,3)=(-((-((-(1-params(37)))*T780))/(T368*T368)));
  v2(177,1)=37;
  v2(177,2)=189;
  v2(177,3)=  v2(176,3);
  v2(178,1)=37;
  v2(178,2)=6375;
  v2(178,3)=(-((-((-(1-params(37)))*T847))/(T368*T368)));
  v2(179,1)=37;
  v2(179,2)=5517;
  v2(179,3)=  v2(178,3);
  v2(180,1)=37;
  v2(180,2)=6482;
  v2(180,3)=(-((-(y(6)*T780))/(T368*T368)));
  v2(181,1)=37;
  v2(181,2)=190;
  v2(181,3)=  v2(180,3);
  v2(182,1)=37;
  v2(182,2)=6519;
  v2(182,3)=(-((-(y(6)*T847))/(T368*T368)));
  v2(183,1)=37;
  v2(183,2)=5518;
  v2(183,3)=  v2(182,3);
  v2(184,1)=37;
  v2(184,2)=6486;
  v2(184,3)=(-(1/T368));
  v2(185,1)=37;
  v2(185,2)=766;
  v2(185,3)=  v2(184,3);
  v2(186,1)=37;
  v2(186,2)=7058;
  v2(186,3)=(-((T368*T368*(-((y(46)*y(6)-(1-params(37))*y(45))*(-((1-params(35))*y(39)))/(y(2)*y(2))))-T1005*T1578)/(T368*T368*T368*T368)));
  v2(187,1)=37;
  v2(187,2)=194;
  v2(187,3)=  v2(186,3);
  v2(188,1)=37;
  v2(188,2)=7095;
  v2(188,3)=(-((T368*T368*(-((y(46)*y(6)-(1-params(37))*y(45))*(1-params(35))/y(2)))-T1005*(T368*T847+T368*T847))/(T368*T368*T368*T368)));
  v2(189,1)=37;
  v2(189,2)=5522;
  v2(189,3)=  v2(188,3);
  v2(190,1)=37;
  v2(190,2)=7062;
  v2(190,3)=(-((-(y(46)*T1003))/(T368*T368)));
  v2(191,1)=37;
  v2(191,2)=770;
  v2(191,3)=  v2(190,3);
  v2(192,1)=37;
  v2(192,2)=7101;
  v2(192,3)=(-((-((-(1-params(37)))*T1003))/(T368*T368)));
  v2(193,1)=37;
  v2(193,2)=6386;
  v2(193,3)=  v2(192,3);
  v2(194,1)=37;
  v2(194,2)=7102;
  v2(194,3)=(-((-(y(6)*T1003))/(T368*T368)));
  v2(195,1)=37;
  v2(195,2)=6530;
  v2(195,3)=  v2(194,3);
  v2(196,1)=37;
  v2(196,2)=7106;
  v2(196,3)=(-((-(T1005*(T368*T1003+T368*T1003)))/(T368*T368*T368*T368)));
  v2(197,1)=38;
  v2(197,2)=146;
  v2(197,3)=(-((T282*T373*T282*T373*(-(y(50)*y(39)*T373*T1419))-(-(y(50)*y(39)*T373*T776))*(T282*T373*T373*T776+T282*T373*T373*T776))/(T282*T373*T282*T373*T282*T373*T282*T373)));
  v2(198,1)=38;
  v2(198,2)=5474;
  v2(198,3)=(-((-(y(50)*T373*T776))/(T282*T373*T282*T373)));
  v2(199,1)=38;
  v2(199,2)=183;
  v2(199,3)=  v2(198,3);
  v2(200,1)=38;
  v2(200,2)=7058;
  v2(200,3)=(-((-(y(39)*T373*T776))/(T282*T373*T282*T373)));
  v2(201,1)=38;
  v2(201,2)=194;
  v2(201,3)=  v2(200,3);
  v2(202,1)=38;
  v2(202,2)=7095;
  v2(202,3)=(-(1/(T282*T373)));
  v2(203,1)=38;
  v2(203,2)=5522;
  v2(203,3)=  v2(202,3);
  v2(204,1)=38;
  v2(204,2)=7202;
  v2(204,3)=(-((T282*T373*T282*T373*(-(y(50)*y(39)*T776*T1017))-(-(y(50)*y(39)*T282*T1017))*(T282*T373*T373*T776+T282*T373*T373*T776))/(T282*T373*T282*T373*T282*T373*T282*T373)));
  v2(205,1)=38;
  v2(205,2)=195;
  v2(205,3)=  v2(204,3);
  v2(206,1)=38;
  v2(206,2)=7239;
  v2(206,3)=(-((-(y(50)*T282*T1017))/(T282*T373*T282*T373)));
  v2(207,1)=38;
  v2(207,2)=5523;
  v2(207,3)=  v2(206,3);
  v2(208,1)=38;
  v2(208,2)=7250;
  v2(208,3)=(-((-(y(39)*T282*T1017))/(T282*T373*T282*T373)));
  v2(209,1)=38;
  v2(209,2)=7107;
  v2(209,3)=  v2(208,3);
  v2(210,1)=38;
  v2(210,2)=7251;
  v2(210,3)=(-((T282*T373*T282*T373*(-(y(50)*y(39)*T282*getPowerDeriv(y(51),params(35),2)))-(-(y(50)*y(39)*T282*T1017))*(T282*T373*T282*T1017+T282*T373*T282*T1017))/(T282*T373*T282*T373*T282*T373*T282*T373)));
  v2(211,1)=38;
  v2(211,2)=11891;
  v2(211,3)=getPowerDeriv(y(83),params(35),2);
  v2(212,1)=39;
  v2(212,2)=5221;
  v2(212,3)=(-(y(69)*y(68)*T797*T1709));
  v2(213,1)=39;
  v2(213,2)=5797;
  v2(213,3)=(-(y(69)*y(68)*(T864*T1709+T798*(-y(45))/(y(41)*y(41)))));
  v2(214,1)=39;
  v2(214,2)=5225;
  v2(214,3)=  v2(213,3);
  v2(215,1)=39;
  v2(215,2)=5801;
  v2(215,3)=(-(y(69)*y(68)*T1729));
  v2(216,1)=39;
  v2(216,2)=6373;
  v2(216,3)=(-(y(69)*y(68)*(T926*T1709+T798*1/y(41))));
  v2(217,1)=39;
  v2(217,2)=5229;
  v2(217,3)=  v2(216,3);
  v2(218,1)=39;
  v2(218,2)=6377;
  v2(218,3)=(-(y(69)*y(68)*(T926*T864*T1708+T798*(-y(37))/(y(41)*y(41)))));
  v2(219,1)=39;
  v2(219,2)=5805;
  v2(219,3)=  v2(218,3);
  v2(220,1)=39;
  v2(220,2)=6381;
  v2(220,3)=(-(y(69)*y(68)*T926*T926*T1708));
  v2(221,1)=39;
  v2(221,2)=9685;
  v2(221,3)=(-(y(69)*T797*T798));
  v2(222,1)=39;
  v2(222,2)=5252;
  v2(222,3)=  v2(221,3);
  v2(223,1)=39;
  v2(223,2)=9689;
  v2(223,3)=(-(y(69)*T798*T864));
  v2(224,1)=39;
  v2(224,2)=5828;
  v2(224,3)=  v2(223,3);
  v2(225,1)=39;
  v2(225,2)=9693;
  v2(225,3)=(-(y(69)*T798*T926));
  v2(226,1)=39;
  v2(226,2)=6404;
  v2(226,3)=  v2(225,3);
  v2(227,1)=39;
  v2(227,2)=9829;
  v2(227,3)=(-(y(68)*T797*T798));
  v2(228,1)=39;
  v2(228,2)=5253;
  v2(228,3)=  v2(227,3);
  v2(229,1)=39;
  v2(229,2)=9833;
  v2(229,3)=(-(y(68)*T798*T864));
  v2(230,1)=39;
  v2(230,2)=5829;
  v2(230,3)=  v2(229,3);
  v2(231,1)=39;
  v2(231,2)=9837;
  v2(231,3)=(-(y(68)*T798*T926));
  v2(232,1)=39;
  v2(232,2)=6405;
  v2(232,3)=  v2(231,3);
  v2(233,1)=39;
  v2(233,2)=9860;
  v2(233,3)=(-T381);
  v2(234,1)=39;
  v2(234,2)=9717;
  v2(234,3)=  v2(233,3);
  v2(235,1)=40;
  v2(235,2)=5221;
  v2(235,3)=(-(T797*T1709));
  v2(236,1)=40;
  v2(236,2)=5797;
  v2(236,3)=(-(T864*T1709+T798*(-y(45))/(y(41)*y(41))));
  v2(237,1)=40;
  v2(237,2)=5225;
  v2(237,3)=  v2(236,3);
  v2(238,1)=40;
  v2(238,2)=5801;
  v2(238,3)=(-T1729);
  v2(239,1)=40;
  v2(239,2)=6373;
  v2(239,3)=(-(T926*T1709+T798*1/y(41)));
  v2(240,1)=40;
  v2(240,2)=5229;
  v2(240,3)=  v2(239,3);
  v2(241,1)=40;
  v2(241,2)=6377;
  v2(241,3)=(-(T926*T864*T1708+T798*(-y(37))/(y(41)*y(41))));
  v2(242,1)=40;
  v2(242,2)=5805;
  v2(242,3)=  v2(241,3);
  v2(243,1)=40;
  v2(243,2)=6381;
  v2(243,3)=(-(T926*T926*T1708));
  v2(244,1)=41;
  v2(244,2)=6373;
  v2(244,3)=T3;
  v2(245,1)=41;
  v2(245,2)=5229;
  v2(245,3)=  v2(244,3);
  v2(246,1)=42;
  v2(246,2)=10151;
  v2(246,3)=(-(1/params(108)*1/params(108)*getPowerDeriv(y(71)/params(108),params(105),2)));
  v2(247,1)=43;
  v2(247,2)=1770;
  v2(247,3)=(-((-y(70))/(y(13)*y(13))));
  v2(248,1)=43;
  v2(248,2)=5917;
  v2(248,3)=  v2(247,3);
  v2(249,1)=43;
  v2(249,2)=1741;
  v2(249,3)=(-(y(42)*(-((-y(70))*(y(13)+y(13))))/(y(13)*y(13)*y(13)*y(13))));
  v2(250,1)=43;
  v2(250,2)=9978;
  v2(250,3)=(-(1/y(13)));
  v2(251,1)=43;
  v2(251,2)=5974;
  v2(251,3)=  v2(250,3);
  v2(252,1)=43;
  v2(252,2)=9949;
  v2(252,3)=(-(y(42)*T3/(y(13)*y(13))));
  v2(253,1)=43;
  v2(253,2)=1798;
  v2(253,3)=  v2(252,3);
  v2(254,1)=44;
  v2(254,2)=1741;
  v2(254,3)=(-((-((-y(70))*(y(13)+y(13))))/(y(13)*y(13)*y(13)*y(13))));
  v2(255,1)=44;
  v2(255,2)=9949;
  v2(255,3)=(-(T3/(y(13)*y(13))));
  v2(256,1)=44;
  v2(256,2)=1798;
  v2(256,3)=  v2(255,3);
  v2(257,1)=45;
  v2(257,2)=5797;
  v2(257,3)=(-((-y(45))/(y(41)*y(41))));
  v2(258,1)=45;
  v2(258,2)=5225;
  v2(258,3)=  v2(257,3);
  v2(259,1)=45;
  v2(259,2)=5801;
  v2(259,3)=(-((-((-(y(45)*y(37)))*(y(41)+y(41))))/(y(41)*y(41)*y(41)*y(41))));
  v2(260,1)=45;
  v2(260,2)=6373;
  v2(260,3)=(-(1/y(41)));
  v2(261,1)=45;
  v2(261,2)=5229;
  v2(261,3)=  v2(260,3);
  v2(262,1)=45;
  v2(262,2)=6377;
  v2(262,3)=(-((-y(37))/(y(41)*y(41))));
  v2(263,1)=45;
  v2(263,2)=5805;
  v2(263,3)=  v2(262,3);
  v2(264,1)=46;
  v2(264,2)=5076;
  v2(264,3)=(-(y(2)*T268*T737*2*T737));
  v2(265,1)=46;
  v2(265,2)=180;
  v2(265,3)=T1392;
  v2(266,1)=46;
  v2(266,2)=5042;
  v2(266,3)=  v2(265,3);
  v2(267,1)=46;
  v2(267,2)=146;
  v2(267,3)=T1401;
  v2(268,1)=47;
  v2(268,2)=4495;
  v2(268,3)=(-(T408*T3/(y(32)*y(32))));
  v2(269,1)=47;
  v2(269,2)=4352;
  v2(269,3)=  v2(268,3);
  v2(270,1)=47;
  v2(270,2)=4496;
  v2(270,3)=(-(T408*(-((-y(31))*(y(32)+y(32))))/(y(32)*y(32)*y(32)*y(32))));
  v2(271,1)=48;
  v2(271,2)=5946;
  v2(271,3)=(-(T884*(-(params(42)*getPowerDeriv(y(42),params(40)-1,2)))/(1-params(42))+T883*T883*getPowerDeriv(T417,1/(1-params(40)),2)));
  v2(272,1)=49;
  v2(272,2)=5448;
  v2(272,3)=(-(T426*T815));
  v2(273,1)=49;
  v2(273,2)=17174;
  v2(273,3)=  v2(272,3);
  v2(274,1)=49;
  v2(274,2)=5366;
  v2(274,3)=(-(y(120)*T426*params(34)*params(42)*(-((-y(122))*(y(38)+y(38))))/(y(38)*y(38)*y(38)*y(38))));
  v2(275,1)=49;
  v2(275,2)=17544;
  v2(275,3)=(-(T426*T826));
  v2(276,1)=49;
  v2(276,2)=17258;
  v2(276,3)=  v2(275,3);
  v2(277,1)=49;
  v2(277,2)=17462;
  v2(277,3)=(-(y(120)*T426*params(34)*params(42)*T3/(y(38)*y(38))));
  v2(278,1)=49;
  v2(278,2)=5450;
  v2(278,3)=  v2(277,3);
  v2(279,1)=49;
  v2(279,2)=17688;
  v2(279,3)=(-T895);
  v2(280,1)=49;
  v2(280,2)=17259;
  v2(280,3)=  v2(279,3);
  v2(281,1)=49;
  v2(281,2)=17606;
  v2(281,3)=(-(y(120)*T815*T894));
  v2(282,1)=49;
  v2(282,2)=5451;
  v2(282,3)=  v2(281,3);
  v2(283,1)=49;
  v2(283,2)=17690;
  v2(283,3)=(-(y(120)*T826*T894));
  v2(284,1)=49;
  v2(284,2)=17547;
  v2(284,3)=  v2(283,3);
  v2(285,1)=49;
  v2(285,2)=17691;
  v2(285,3)=(-(y(120)*T425*getPowerDeriv(y(123),params(40),2)));
  v2(286,1)=49;
  v2(286,2)=7096;
  v2(286,3)=T3;
  v2(287,1)=49;
  v2(287,2)=5666;
  v2(287,3)=  v2(286,3);
  v2(288,1)=50;
  v2(288,2)=5449;
  v2(288,3)=(-(T432*T815));
  v2(289,1)=50;
  v2(289,2)=17318;
  v2(289,3)=  v2(288,3);
  v2(290,1)=50;
  v2(290,2)=5366;
  v2(290,3)=(-(y(121)*T432*params(34)*params(42)*(-((-y(122))*(y(38)+y(38))))/(y(38)*y(38)*y(38)*y(38))));
  v2(291,1)=50;
  v2(291,2)=17545;
  v2(291,3)=(-(T432*T826));
  v2(292,1)=50;
  v2(292,2)=17402;
  v2(292,3)=  v2(291,3);
  v2(293,1)=50;
  v2(293,2)=17462;
  v2(293,3)=(-(y(121)*T432*params(34)*params(42)*T3/(y(38)*y(38))));
  v2(294,1)=50;
  v2(294,2)=5450;
  v2(294,3)=  v2(293,3);
  v2(295,1)=50;
  v2(295,2)=17689;
  v2(295,3)=(-T899);
  v2(296,1)=50;
  v2(296,2)=17403;
  v2(296,3)=  v2(295,3);
  v2(297,1)=50;
  v2(297,2)=17606;
  v2(297,3)=(-(y(121)*T815*T898));
  v2(298,1)=50;
  v2(298,2)=5451;
  v2(298,3)=  v2(297,3);
  v2(299,1)=50;
  v2(299,2)=17690;
  v2(299,3)=(-(y(121)*T826*T898));
  v2(300,1)=50;
  v2(300,2)=17547;
  v2(300,3)=  v2(299,3);
  v2(301,1)=50;
  v2(301,2)=17691;
  v2(301,3)=(-(y(121)*T425*getPowerDeriv(y(123),params(40)-1,2)));
  v2(302,1)=51;
  v2(302,2)=7531;
  v2(302,3)=(-(T3/(y(53)*y(53))));
  v2(303,1)=51;
  v2(303,2)=6101;
  v2(303,3)=  v2(302,3);
  v2(304,1)=51;
  v2(304,2)=7541;
  v2(304,3)=(-((-((-y(43))*(y(53)+y(53))))/(y(53)*y(53)*y(53)*y(53))));
  v2(305,1)=52;
  v2(305,2)=8909;
  v2(305,3)=(-(params(34)*T871));
  v2(306,1)=52;
  v2(306,2)=17918;
  v2(306,3)=  v2(305,3);
  v2(307,1)=53;
  v2(307,2)=1306;
  v2(307,3)=(-(y(59)*(1-params(35))*(-((-y(63))*(y(10)+y(10))))/(y(10)*y(10)*y(10)*y(10))/y(11)));
  v2(308,1)=53;
  v2(308,2)=8362;
  v2(308,3)=(-((1-params(35))*(-y(63))/(y(10)*y(10))/y(11)));
  v2(309,1)=53;
  v2(309,2)=1355;
  v2(309,3)=  v2(308,3);
  v2(310,1)=53;
  v2(310,2)=1450;
  v2(310,3)=(-((-(y(59)*(1-params(35))*(-y(63))/(y(10)*y(10))))/(y(11)*y(11))));
  v2(311,1)=53;
  v2(311,2)=1307;
  v2(311,3)=  v2(310,3);
  v2(312,1)=53;
  v2(312,2)=1499;
  v2(312,3)=(-((-((1-params(35))*y(63)/y(10)))/(y(11)*y(11))));
  v2(313,1)=53;
  v2(313,2)=8363;
  v2(313,3)=  v2(312,3);
  v2(314,1)=53;
  v2(314,2)=1451;
  v2(314,3)=(-((-((-((1-params(35))*y(63)/y(10)*y(59)+(1-params(37))*y(60)))*(y(11)+y(11))))/(y(11)*y(11)*y(11)*y(11))));
  v2(315,1)=53;
  v2(315,2)=8507;
  v2(315,3)=(-((-(1-params(37)))/(y(11)*y(11))));
  v2(316,1)=53;
  v2(316,2)=1500;
  v2(316,3)=  v2(315,3);
  v2(317,1)=53;
  v2(317,2)=8938;
  v2(317,3)=(-(y(59)*(1-params(35))*T3/(y(10)*y(10))/y(11)));
  v2(318,1)=53;
  v2(318,2)=1359;
  v2(318,3)=  v2(317,3);
  v2(319,1)=53;
  v2(319,2)=8987;
  v2(319,3)=(-((1-params(35))*T1043/y(11)));
  v2(320,1)=53;
  v2(320,2)=8415;
  v2(320,3)=  v2(319,3);
  v2(321,1)=53;
  v2(321,2)=8939;
  v2(321,3)=(-((-(y(59)*(1-params(35))*T1043))/(y(11)*y(11))));
  v2(322,1)=53;
  v2(322,2)=1503;
  v2(322,3)=  v2(321,3);
  v2(323,1)=56;
  v2(323,2)=1351;
  v2(323,3)=(-(params(41)*T3/(y(10)*y(10))));
  v2(324,1)=56;
  v2(324,2)=7786;
  v2(324,3)=  v2(323,3);
  v2(325,1)=56;
  v2(325,2)=1306;
  v2(325,3)=(-(params(41)*(-((-y(55))*(y(10)+y(10))))/(y(10)*y(10)*y(10)*y(10))));
  v2(326,1)=56;
  v2(326,2)=8545;
  v2(326,3)=1;
  v2(327,1)=56;
  v2(327,2)=6972;
  v2(327,3)=  v2(326,3);
  v2(328,1)=57;
  v2(328,2)=7825;
  v2(328,3)=T3;
  v2(329,1)=57;
  v2(329,2)=6967;
  v2(329,3)=  v2(328,3);
  v2(330,1)=58;
  v2(330,2)=7831;
  v2(330,3)=(-(y(10)*T268*T1043*2*T1043));
  v2(331,1)=58;
  v2(331,2)=1351;
  v2(331,3)=(-(T268*T1043*T1046+y(10)*T268*(T1059*2*T1043+T1046*T3/(y(10)*y(10)))));
  v2(332,1)=58;
  v2(332,2)=7786;
  v2(332,3)=  v2(331,3);
  v2(333,1)=58;
  v2(333,2)=1306;
  v2(333,3)=(-(T1063+T1063+y(10)*T268*(T1059*2*T1059+T1046*(-((-y(55))*(y(10)+y(10))))/(y(10)*y(10)*y(10)*y(10)))));
  v2(334,1)=59;
  v2(334,2)=7251;
  v2(334,3)=(-(T490*y(54)*T1927));
  v2(335,1)=59;
  v2(335,2)=7683;
  v2(335,3)=(-(T490*(T1023+y(51)*T1927)));
  v2(336,1)=59;
  v2(336,2)=7254;
  v2(336,3)=  v2(335,3);
  v2(337,1)=59;
  v2(337,2)=7686;
  v2(337,3)=(-(T490*y(51)*y(51)*T1926));
  v2(338,1)=59;
  v2(338,2)=1347;
  v2(338,3)=(-(y(54)*T1023*T1067));
  v2(339,1)=59;
  v2(339,2)=7210;
  v2(339,3)=  v2(338,3);
  v2(340,1)=59;
  v2(340,2)=1350;
  v2(340,3)=(-(y(51)*T1023*T1067));
  v2(341,1)=59;
  v2(341,2)=7642;
  v2(341,3)=  v2(340,3);
  v2(342,1)=59;
  v2(342,2)=1306;
  v2(342,3)=(-(T489*getPowerDeriv(y(10),1-params(35),2)));
  v2(343,1)=59;
  v2(343,2)=8986;
  v2(343,3)=1;
  v2(344,1)=59;
  v2(344,2)=8271;
  v2(344,3)=  v2(343,3);
  v2(345,1)=60;
  v2(345,2)=7686;
  v2(345,3)=params(39)*getPowerDeriv(y(54),params(38),2)-y(59)*params(35)*y(57)*(-((-(y(63)*y(58)))*(y(54)+y(54))))/(y(54)*y(54)*y(54)*y(54));
  v2(346,1)=60;
  v2(346,2)=8118;
  v2(346,3)=(-(y(59)*params(35)*(-(y(63)*y(58)))/(y(54)*y(54))));
  v2(347,1)=60;
  v2(347,2)=7689;
  v2(347,3)=  v2(346,3);
  v2(348,1)=60;
  v2(348,2)=8262;
  v2(348,3)=(-(y(59)*params(35)*y(57)*(-y(63))/(y(54)*y(54))));
  v2(349,1)=60;
  v2(349,2)=7690;
  v2(349,3)=  v2(348,3);
  v2(350,1)=60;
  v2(350,2)=8265;
  v2(350,3)=(-(y(59)*params(35)*y(63)/y(54)));
  v2(351,1)=60;
  v2(351,2)=8122;
  v2(351,3)=  v2(350,3);
  v2(352,1)=60;
  v2(352,2)=8406;
  v2(352,3)=(-(params(35)*y(57)*(-(y(63)*y(58)))/(y(54)*y(54))));
  v2(353,1)=60;
  v2(353,2)=7691;
  v2(353,3)=  v2(352,3);
  v2(354,1)=60;
  v2(354,2)=8409;
  v2(354,3)=(-(params(35)*y(63)*y(58)/y(54)));
  v2(355,1)=60;
  v2(355,2)=8123;
  v2(355,3)=  v2(354,3);
  v2(356,1)=60;
  v2(356,2)=8410;
  v2(356,3)=(-(params(35)*y(57)*y(63)/y(54)));
  v2(357,1)=60;
  v2(357,2)=8267;
  v2(357,3)=  v2(356,3);
  v2(358,1)=60;
  v2(358,2)=8982;
  v2(358,3)=(-(y(59)*params(35)*y(57)*(-y(58))/(y(54)*y(54))));
  v2(359,1)=60;
  v2(359,2)=7695;
  v2(359,3)=  v2(358,3);
  v2(360,1)=60;
  v2(360,2)=8985;
  v2(360,3)=(-(y(59)*params(35)*y(58)/y(54)));
  v2(361,1)=60;
  v2(361,2)=8127;
  v2(361,3)=  v2(360,3);
  v2(362,1)=60;
  v2(362,2)=8986;
  v2(362,3)=(-(y(59)*params(35)*y(57)*1/y(54)));
  v2(363,1)=60;
  v2(363,2)=8271;
  v2(363,3)=  v2(362,3);
  v2(364,1)=60;
  v2(364,2)=8987;
  v2(364,3)=(-(params(35)*y(57)*y(58)/y(54)));
  v2(365,1)=60;
  v2(365,2)=8415;
  v2(365,3)=  v2(364,3);
  v2(366,1)=63;
  v2(366,2)=9136;
  v2(366,3)=(-(T3/(y(64)*y(64))));
  v2(367,1)=66;
  v2(367,2)=1161;
  v2(367,3)=(-(params(47)*T3/(y(9)*y(9))));
  v2(368,1)=66;
  v2(368,2)=7251;
  v2(368,3)=T3/(y(51)*y(51));
  v2(369,1)=67;
  v2(369,2)=1016;
  v2(369,3)=(-(params(49)*T3/(y(8)*y(8))));
  v2(370,1)=67;
  v2(370,2)=6961;
  v2(370,3)=T3/(y(49)*y(49));
  v2(371,1)=68;
  v2(371,2)=581;
  v2(371,3)=(-(params(52)*T3/(y(5)*y(5))));
  v2(372,1)=68;
  v2(372,2)=6091;
  v2(372,3)=T3/(y(43)*y(43));
  v2(373,1)=69;
  v2(373,2)=1596;
  v2(373,3)=(-(params(103)*T3/(y(12)*y(12))));
  v2(374,1)=69;
  v2(374,2)=9716;
  v2(374,3)=T3/(y(68)*y(68));
  g2 = sparse(v2(:,1),v2(:,2),v2(:,3),90,20736);
if nargout >= 4,
  %
  % Third order derivatives
  %

  g3 = sparse([],[],[],90,2985984);
end
end
end
end
