function [residual, g1, g2] = quest1_dynamic(y, x)
%
% Status : Computes dynamic model for Dynare
%
%Warning : this file is generated automatically by Dynare
%			from model file (.modglobal M_ 
if M_.param_nbr > 0
  params = M_.params;
end
  residual = zeros( 102, 1);

	%
	% Model equations
	%

T25=exp(y(44))-((params(46)*exp(y(44)))/(1.0+y(34))*(1.0+y(64)))/(1.0+y(63));
lhs =exp(y(56));
rhs =((1.0-params(46))*exp(y(84)))/T25;
residual(1)= lhs-rhs;
lhs =y(79);
rhs =exp(y(92))*params(61)*((1.0-exp(y(40))-params(47)*(1.0-exp(y(40))))/(1.0-params(47)))^(-params(55));
residual(2)= lhs-rhs;
lhs =y(56)-y(56);
rhs =y(63)+(y(34)+y(31))-y(64);
residual(3)= lhs-rhs;
lhs =(1.0+params(125))*exp(y(43));
rhs =(exp(y(40))*(1.0-y(77)-params(109))/exp(y(60))+exp(y(55)))-y(73);
residual(4)= lhs-rhs;
lhs =y(82);
rhs =exp(y(40))/exp(y(60));
residual(5)= lhs-rhs;
lhs =y(3);
rhs =exp(y(43));
residual(6)= lhs-rhs;
lhs =y(76);
rhs =exp(y(55));
residual(7)= lhs-rhs;
lhs =y(75);
rhs =exp(y(55))-y(73);
residual(8)= lhs-rhs;
lhs =y(1);
rhs =exp(y(41));
residual(9)= lhs-rhs;
lhs =y(26);
rhs =exp(y(46));
residual(10)= lhs-rhs;
lhs =y(83);
rhs =(1.0-y(77)-params(109))*y(82);
residual(11)= lhs-rhs;
lhs =exp(y(42));
rhs =exp(y(43))*params(106)+exp(y(44))*(1.0-params(106));
residual(12)= lhs-rhs;
lhs =exp(y(49));
rhs =(1.0-params(107))*((-(params(31)/params(32)))*(y(18)-params(44)-params(38))+(params(31)*params(6))/params(32)*(y(18)-params(44)-params(38))+exp(y(89))/params(32)*(y(71)-1.0))+params(107)*params(62)*(1.0-params(114))*(exp(y(59))-exp(y(59))/exp(y(60))-params(103)*(y(38)+params(96)+y(99)));
residual(13)= lhs-rhs;
lhs =exp(y(59))*(1.0-params(114))*y(8)*(1.0-params(5));
rhs =y(71)-((1.0-y(72)-params(14)-params(96)-y(99)+y(64))-y(63)-y(65))*y(71)+params(2)*(y(78)-1.0)+params(3)*(y(78)-1.0)^2.0;
residual(14)= lhs-rhs;
lhs =1.0/params(6)-1.0;
rhs =(y(38)+y(31))-y(64);
residual(15)= lhs-rhs;
lhs =y(62);
rhs =(1.0+y(38))^(-params(134));
residual(16)= lhs-rhs;
lhs =y(23);
rhs =exp(y(49))-(params(14)+y(25))/(1.0+y(25));
residual(17)= lhs-rhs;
lhs =y(46);
rhs =(1.0-params(41))*log(params(42))+params(41)*y(46)+params(28)*y(57)+params(29)*y(57)+x(6);
residual(18)= lhs-rhs;
lhs =y(47);
rhs =(1.0-params(48))*log(params(53))+params(48)*y(47)+y(57)*params(52)+y(90);
residual(19)= lhs-rhs;
lhs =y(90);
rhs =params(71)*y(90)+x(8);
residual(20)= lhs-rhs;
lhs =y(19)-y(18);
rhs =y(47)-y(50)-y(47)+y(50);
residual(21)= lhs-rhs;
lhs =y(88);
rhs =x(6);
residual(22)= lhs-rhs;
lhs =y(55);
rhs =log(params(121))+params(118)*(y(40)-log(params(56)))+params(119)*(y(40)-log(params(56)))+params(120)*(y(4)-params(16))+y(85);
residual(23)= lhs-rhs;
lhs =y(85);
rhs =params(122)*y(85)+x(23);
residual(24)= lhs-rhs;
lhs =1.0;
rhs =(exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47))+exp(y(45)))-exp(y(48));
residual(25)= lhs-rhs;
lhs =y(74);
rhs =exp(y(45))-exp(y(48));
residual(26)= lhs-rhs;
lhs =y(24);
rhs =y(40)-y(40);
residual(27)= lhs-rhs;
lhs =(y(64)+y(16)-y(34))-y(63);
rhs =y(46)-y(46);
residual(28)= lhs-rhs;
lhs =(y(64)+y(9)-y(34))-y(63);
rhs =y(42)-y(42);
residual(29)= lhs-rhs;
lhs =(y(64)+y(11)-y(34))-y(63);
rhs =y(43)-y(43);
residual(30)= lhs-rhs;
lhs =(y(64)+y(12)-y(34))-y(63);
rhs =y(44)-y(44);
residual(31)= lhs-rhs;
lhs =y(50);
rhs =(y(49)+(-y(59)))-y(23);
residual(32)= lhs-rhs;
lhs =y(34);
rhs =(1.0-params(5))*(y(23)+y(32))+params(5)*(y(28)+y(24)*(1.0+params(58)));
residual(33)= lhs-rhs;
lhs =y(28)-params(43);
rhs =y(100)+x(25);
residual(34)= lhs-rhs;
lhs =y(100);
rhs =params(83)*y(100)+x(22);
residual(35)= lhs-rhs;
lhs =exp(y(101))/exp(y(60));
rhs =(y(8)*params(5))/exp(y(40))*(1.0+y(51))-(y(40)-y(40))*(1.0-params(113))*params(33)+((1.0-params(113))*params(33))/(1.0+y(72))*(y(40)-y(40));
residual(36)= lhs-rhs;
lhs =exp(y(59))*y(8)*(1.0-params(5));
rhs =y(78)*(params(2)+(y(78)-1.0)*params(3)*2.0);
residual(37)= lhs-rhs;
lhs =y(8);
rhs =1.0-(params(113)+y(86))-params(34)*(params(6)*((y(63)*params(100)+(1.0-params(100))*y(63))-params(37))-(y(63)-params(37)));
residual(38)= lhs-rhs;
lhs =y(80);
rhs =y(63)+y(81);
residual(39)= lhs-rhs;
lhs =y(72);
rhs =y(38)-y(63);
residual(40)= lhs-rhs;
lhs =y(57);
rhs =(1.0-params(5))*log(y(78))+params(5)*(y(40)-log(params(56)));
residual(41)= lhs-rhs;
lhs =y(38);
rhs =params(49)*y(38)+(1.0-params(49))*(params(44)+params(24)+y(102))+params(116)*(y(64)-y(102))+y(57)*params(123)+params(124)*y(58)+y(93);
residual(42)= lhs-rhs;
lhs =y(93);
rhs =params(75)*y(93)+x(13);
residual(43)= lhs-rhs;
lhs =y(57)-y(57);
rhs =y(34)-y(36);
residual(44)= lhs-rhs;
T448=exp(y(56))/(1.0+params(125));
lhs =y(79);
rhs =T448*((1.0-y(77)-params(109))*(1.0-params(115))/(-params(115))*exp((-y(60)))+params(36)/params(115)*(y(80)-params(44)-params(37)))-(params(36)*params(6)*exp(y(56)))/params(115)*((params(102)*y(80)+(1.0-params(102))*y(80))-params(44)-params(37));
residual(45)= lhs-rhs;
lhs =y(52);
rhs =(1.0-params(97))*y(53)-(params(97)*(1.0-params(97))*y(53)^2.0*(params(105)-1.0))/2.0+y(94);
residual(46)= lhs-rhs;
lhs =y(52)-y(52);
rhs =y(64)-y(63);
residual(47)= lhs-rhs;
lhs =y(53)-y(53);
rhs =y(66)-y(63);
residual(48)= lhs-rhs;
lhs =y(54)-y(54);
rhs =y(69)-y(63);
residual(49)= lhs-rhs;
T518=(1.0-params(97))*(exp(y(52))/exp(y(53)))^(params(105)-1.0);
T532=((1.0+y(21))-params(44))^(1.0-params(99));
T533=((1.0+y(21))-params(44))^params(99)/T532;
T536=T533^(params(6)*params(50));
lhs =exp(y(48));
rhs =(exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*T518*exp(y(91))*T536;
residual(50)= lhs-rhs;
T550=exp(y(7))^(params(110)*params(4))/exp(y(54));
T557=exp(y(61))^params(4);
T558=exp(y(54))*(1.0-params(110))*exp(y(87))*T550^params(104)*T557;
T568=((1.0+y(14))-params(44))^(1.0-params(98));
T569=((1.0+y(14))-params(44))^params(98)/T568;
T572=T569^(params(6)*params(26));
lhs =exp(y(45));
rhs =T558*T572;
residual(51)= lhs-rhs;
lhs =y(38);
rhs =(y(39)+y(13))-params(95)*y(2)+y(98);
residual(52)= lhs-rhs;
lhs =y(2);
rhs =(exp(y(45))+((1.0+y(72))-y(34)-y(25))*y(2))-exp(y(48));
residual(53)= lhs-rhs;
lhs =y(45)-y(45);
rhs =(y(69)-y(63)+y(14))-y(34);
residual(54)= lhs-rhs;
lhs =y(48)-y(48);
rhs =(y(66)-y(63)+y(21))-y(34);
residual(55)= lhs-rhs;
T606=(exp(y(53))/(1.0+y(63)))^(1.0-params(108));
T608=exp(y(7))^params(4)^params(108);
lhs =exp(y(53));
rhs =T606*T608*exp(y(96));
residual(56)= lhs-rhs;
lhs =y(7)-y(7);
rhs =(y(68)+y(13))-y(63);
residual(57)= lhs-rhs;
T639=exp(y(52))^(1.0-params(111));
lhs =exp(y(54));
rhs =exp(y(97))*(1.0-x(4)+params(35)*(params(37)+params(6)*((params(101)*y(69)+(1.0-params(101))*y(69))-params(37))-y(69)))*T639;
residual(58)= lhs-rhs;
lhs =exp(y(41));
rhs =(exp(y(55))+exp(y(47))+exp(y(46))+((1.0+y(72))-y(34)-y(25))*exp(y(41)))-exp(y(40))*(y(77)+params(109))/exp(y(60))-params(114)*(1.0-exp(y(40))/exp(y(60)))-params(125)*exp(y(42))-y(73);
residual(59)= lhs-rhs;
lhs =y(73)-y(73);
rhs =params(7)*(exp(y(41))-params(9))+params(8)*(exp(y(41))-exp(y(41)));
residual(60)= lhs-rhs;
lhs =y(39);
rhs =(1.0-params(85))*params(23)+params(85)*y(39)+params(86)*(y(68)-params(40))+params(87)*(y(37)-params(45))+x(10);
residual(61)= lhs-rhs;
lhs =y(68)-params(40);
rhs =params(88)*(y(39)-params(23))+(y(68)-params(40))*params(89)+(y(37)-params(45))*params(90)+x(18);
residual(62)= lhs-rhs;
lhs =y(37)-params(45);
rhs =(y(39)-params(23))*params(91)+(y(68)-params(40))*params(92)+(y(37)-params(45))*params(93)+params(94)*(y(61)-params(60))+x(26);
residual(63)= lhs-rhs;
lhs =y(61)-y(61);
rhs =y(37)-y(34);
residual(64)= lhs-rhs;
lhs =y(84);
rhs =params(63)*y(84)+x(1);
residual(65)= lhs-rhs;
lhs =y(86);
rhs =params(67)*y(86)+x(3);
residual(66)= lhs-rhs;
lhs =y(87);
rhs =params(68)*y(87)+x(5);
residual(67)= lhs-rhs;
lhs =y(89);
rhs =params(70)*y(89)+x(7);
residual(68)= lhs-rhs;
lhs =y(91);
rhs =params(72)*y(91)+x(9);
residual(69)= lhs-rhs;
lhs =y(51)-params(58);
rhs =params(74)*(y(51)-params(58))+x(12);
residual(70)= lhs-rhs;
lhs =y(92);
rhs =params(73)*y(92)+x(11);
residual(71)= lhs-rhs;
lhs =y(98);
rhs =params(81)*y(98)+x(20);
residual(72)= lhs-rhs;
lhs =y(99);
rhs =params(82)*y(99)+x(21);
residual(73)= lhs-rhs;
lhs =y(101);
rhs =x(24);
residual(74)= lhs-rhs;
lhs =y(102)-params(37);
rhs =params(51)*(y(102)-params(37));
residual(75)= lhs-rhs;
lhs =y(59)-y(59);
rhs =(y(63)+y(34)-y(23))-y(64)+y(65);
residual(76)= lhs-rhs;
lhs =y(65);
rhs =params(38)+y(95);
residual(77)= lhs-rhs;
lhs =y(95);
rhs =params(77)*y(95)+x(15);
residual(78)= lhs-rhs;
lhs =y(49)-y(49);
rhs =y(18)-y(23);
residual(79)= lhs-rhs;
lhs =y(60)-y(60);
rhs =y(63)+y(34)+(-y(80));
residual(80)= lhs-rhs;
lhs =log(y(78))-log(y(78));
rhs =y(32);
residual(81)= lhs-rhs;
lhs =y(94);
rhs =params(76)*y(94)+x(14);
residual(82)= lhs-rhs;
lhs =y(96);
rhs =params(78)*y(96)+x(16);
residual(83)= lhs-rhs;
lhs =y(97);
rhs =params(80)*y(97)+x(19);
residual(84)= lhs-rhs;
lhs =y(29);
rhs =(1.0-params(5))*y(32)+params(5)*y(28);
residual(85)= lhs-rhs;
lhs =y(33);
rhs =(-y(60))+y(60);
residual(86)= lhs-rhs;
lhs =y(77);
rhs =params(126)*(1.0+y(57)*params(1));
residual(87)= lhs-rhs;
lhs =y(30)-y(34)-y(63);
rhs =y(55)-y(55);
residual(88)= lhs-rhs;
lhs =y(25)-params(39);
rhs =(y(25)-params(39))*params(79)+x(17);
residual(89)= lhs-rhs;
lhs =y(4)-params(16);
rhs =params(64)*(y(4)-params(16))+params(65)*(y(5)-params(16))+params(66)*(y(6)-params(16))+x(2);
residual(90)= lhs-rhs;
lhs =log(y(73))-log(y(73));
rhs =y(27)-y(34)-y(63);
residual(91)= lhs-rhs;
lhs =y(35);
rhs =y(34)+y(25);
residual(92)= lhs-rhs;
lhs =y(10);
rhs =y(25)+y(9);
residual(93)= lhs-rhs;
lhs =y(20);
rhs =y(18)+y(25);
residual(94)= lhs-rhs;
lhs =y(17);
rhs =y(25)+y(16);
residual(95)= lhs-rhs;
lhs =y(15);
rhs =y(25)+y(14)+params(17);
residual(96)= lhs-rhs;
lhs =y(22);
rhs =y(25)+y(21)+params(18);
residual(97)= lhs-rhs;
lhs =y(67);
rhs =y(66)+params(19);
residual(98)= lhs-rhs;
lhs =y(70);
rhs =y(69)+params(20);
residual(99)= lhs-rhs;
lhs =y(58);
rhs =y(57);
residual(100)= lhs-rhs;
lhs =y(5);
rhs =y(4);
residual(101)= lhs-rhs;
lhs =y(6);
rhs =y(5);
residual(102)= lhs-rhs;
T938=exp(y(60))*exp(y(60));
T952=T532*T532;
T954=T568*T568;
if ~isreal(residual)
  residual = real(residual)+imag(residual).^2;
end
if nargout >= 2,
  g1 = zeros(102, 102);

	%
	% Jacobian matrix
	%

T1088=exp(y(40))*exp(y(60))/T938;
  g1(1, 56)=  g1(1, 56)+exp(y(56));
  g1(3, 56)=  g1(3, 56)+1.0;
  g1(45, 56)=  g1(45, 56)+(-(((1.0-y(77)-params(109))*(1.0-params(115))/(-params(115))*exp((-y(60)))+params(36)/params(115)*(y(80)-params(44)-params(37)))*(exp(y(56))*(1.0+params(125)))/((1.0+params(125))*(1.0+params(125)))));
  g1(1, 84)=  g1(1, 84)+(-(((1.0-params(46))*exp(y(84))*T25)/(T25*T25)));
  g1(65, 84)=  g1(65, 84)+1.0;
  g1(1, 44)=  g1(1, 44)+(-((-((1.0-params(46))*exp(y(84))*exp(y(44))))/(T25*T25)));
  g1(12, 44)=  g1(12, 44)+(-(exp(y(44))*(1.0-params(106))));
  g1(31, 44)=  g1(31, 44)+(-1.0);
  g1(1, 44)=  g1(1, 44)+(-((-((1.0-params(46))*exp(y(84))*(-(((1.0+y(63))*(1.0+y(64))*(params(46)*exp(y(44))*(1.0+y(34)))/((1.0+y(34))*(1.0+y(34))))/((1.0+y(63))*(1.0+y(63)))))))/(T25*T25)));
  g1(31, 44)=  g1(31, 44)+1.0;
  g1(1, 34)=  g1(1, 34)+(-((-((1.0-params(46))*exp(y(84))*(-(((1.0+y(63))*(1.0+y(64))*(-(params(46)*exp(y(44))))/((1.0+y(34))*(1.0+y(34))))/((1.0+y(63))*(1.0+y(63)))))))/(T25*T25)));
  g1(3, 34)=  g1(3, 34)+(-1.0);
  g1(28, 34)=  g1(28, 34)+(-1.0);
  g1(29, 34)=  g1(29, 34)+(-1.0);
  g1(30, 34)=  g1(30, 34)+(-1.0);
  g1(31, 34)=  g1(31, 34)+(-1.0);
  g1(33, 34)=  g1(33, 34)+1.0;
  g1(44, 34)=  g1(44, 34)+(-1.0);
  g1(53, 34)=  g1(53, 34)+y(2);
  g1(54, 34)=  g1(54, 34)+1.0;
  g1(55, 34)=  g1(55, 34)+1.0;
  g1(59, 34)=  g1(59, 34)+exp(y(41));
  g1(64, 34)=  g1(64, 34)+1.0;
  g1(76, 34)=  g1(76, 34)+(-1.0);
  g1(80, 34)=  g1(80, 34)+(-1.0);
  g1(88, 34)=  g1(88, 34)+(-1.0);
  g1(91, 34)=  g1(91, 34)+1.0;
  g1(92, 34)=  g1(92, 34)+(-1.0);
  g1(1, 64)=  g1(1, 64)+(-((-((1.0-params(46))*exp(y(84))*(-(((params(46)*exp(y(44)))/(1.0+y(34))*(1.0+y(63)))/((1.0+y(63))*(1.0+y(63)))))))/(T25*T25)));
  g1(3, 64)=  g1(3, 64)+1.0;
  g1(28, 64)=  g1(28, 64)+1.0;
  g1(29, 64)=  g1(29, 64)+1.0;
  g1(30, 64)=  g1(30, 64)+1.0;
  g1(31, 64)=  g1(31, 64)+1.0;
  g1(42, 64)=  g1(42, 64)+(-params(116));
  g1(47, 64)=  g1(47, 64)+(-1.0);
  g1(76, 64)=  g1(76, 64)+1.0;
  g1(1, 63)=  g1(1, 63)+(-((-((1.0-params(46))*exp(y(84))*(-((-((params(46)*exp(y(44)))/(1.0+y(34))*(1.0+y(64))))/((1.0+y(63))*(1.0+y(63)))))))/(T25*T25)));
  g1(3, 63)=  g1(3, 63)+(-1.0);
  g1(28, 63)=  g1(28, 63)+(-1.0);
  g1(29, 63)=  g1(29, 63)+(-1.0);
  g1(30, 63)=  g1(30, 63)+(-1.0);
  g1(31, 63)=  g1(31, 63)+(-1.0);
  g1(38, 63)=  g1(38, 63)+(-params(34));
  g1(39, 63)=  g1(39, 63)+(-1.0);
  g1(47, 63)=  g1(47, 63)+1.0;
  g1(48, 63)=  g1(48, 63)+1.0;
  g1(49, 63)=  g1(49, 63)+1.0;
  g1(54, 63)=  g1(54, 63)+1.0;
  g1(55, 63)=  g1(55, 63)+1.0;
  g1(56, 63)=  g1(56, 63)+(-(exp(y(96))*T608*(-exp(y(53)))/((1.0+y(63))*(1.0+y(63)))*(1.0-params(108))*(exp(y(53))/(1.0+y(63)))^(1.0-params(108)-1.0)));
  g1(57, 63)=  g1(57, 63)+1.0;
  g1(76, 63)=  g1(76, 63)+(-1.0);
  g1(80, 63)=  g1(80, 63)+(-1.0);
  g1(88, 63)=  g1(88, 63)+(-1.0);
  g1(91, 63)=  g1(91, 63)+1.0;
  g1(2, 79)=  g1(2, 79)+1.0;
  g1(45, 79)=  g1(45, 79)+1.0;
  g1(2, 92)=  g1(2, 92)+(-(exp(y(92))*params(61)*((1.0-exp(y(40))-params(47)*(1.0-exp(y(40))))/(1.0-params(47)))^(-params(55))));
  g1(71, 92)=  g1(71, 92)+1.0;
  g1(2, 40)=  g1(2, 40)+(-(exp(y(92))*params(61)*((1.0-params(47))*(-exp(y(40))))/((1.0-params(47))*(1.0-params(47)))*(-params(55))*((1.0-exp(y(40))-params(47)*(1.0-exp(y(40))))/(1.0-params(47)))^((-params(55))-1.0)));
  g1(4, 40)=  g1(4, 40)+(-(exp(y(40))*(1.0-y(77)-params(109))/exp(y(60))));
  g1(5, 40)=  g1(5, 40)+(-((exp(y(40))*exp(y(60)))/T938));
  g1(23, 40)=  g1(23, 40)+(-params(118));
  g1(27, 40)=  g1(27, 40)+(-1.0);
  g1(36, 40)=  g1(36, 40)+(-((1.0+y(51))*(-(exp(y(40))*y(8)*params(5)))/(exp(y(40))*exp(y(40)))-(1.0-params(113))*params(33)+(-(((1.0-params(113))*params(33))/(1.0+y(72))))));
  g1(41, 40)=  g1(41, 40)+(-params(5));
  g1(59, 40)=  g1(59, 40)+(-((-(exp(y(40))*(y(77)+params(109))/exp(y(60))))-params(114)*(-((exp(y(40))*exp(y(60)))/T938))));
  g1(2, 40)=  g1(2, 40)+(-(exp(y(92))*params(61)*(-params(55))*((1.0-exp(y(40))-params(47)*(1.0-exp(y(40))))/(1.0-params(47)))^((-params(55))-1.0)*((1.0-params(47))*(-(params(47)*(-exp(y(40))))))/((1.0-params(47))*(1.0-params(47)))));
  g1(23, 40)=  g1(23, 40)+(-params(119));
  g1(27, 40)=  g1(27, 40)+1.0;
  g1(36, 40)=  g1(36, 40)+(-((1.0-params(113))*params(33)));
  g1(3, 56)=  g1(3, 56)+(-1.0);
  g1(3, 31)=  g1(3, 31)+(-1.0);
  g1(4, 43)=  g1(4, 43)+(1.0+params(125))*exp(y(43));
  g1(6, 43)=  g1(6, 43)+(-exp(y(43)));
  g1(12, 43)=  g1(12, 43)+(-(exp(y(43))*params(106)));
  g1(30, 43)=  g1(30, 43)+(-1.0);
  g1(4, 77)=  g1(4, 77)+(-(exp(y(40))*(-exp(y(60)))/T938));
  g1(11, 77)=  g1(11, 77)+y(82);
  g1(45, 77)=  g1(45, 77)+(-(T448*(-((1.0-params(115))/(-params(115))*exp((-y(60)))))));
  g1(59, 77)=  g1(59, 77)+T1088;
  g1(87, 77)=  g1(87, 77)+1.0;
  g1(4, 60)=  g1(4, 60)+(-(exp(y(40))*(-((1.0-y(77)-params(109))*exp(y(60))))/T938));
  g1(5, 60)=  g1(5, 60)+(-((-(exp(y(40))*exp(y(60))))/T938));
  g1(13, 60)=  g1(13, 60)+(-(params(107)*params(62)*(1.0-params(114))*(-((-(exp(y(60))*exp(y(59))))/T938))));
  g1(36, 60)=  g1(36, 60)+(-(exp(y(60))*exp(y(101))))/T938;
  g1(45, 60)=  g1(45, 60)+(-(T448*(1.0-y(77)-params(109))*(1.0-params(115))/(-params(115))*(-exp((-y(60))))));
  g1(59, 60)=  g1(59, 60)+(-((-(exp(y(40))*(-(exp(y(60))*(y(77)+params(109))))/T938))-params(114)*(-((-(exp(y(40))*exp(y(60))))/T938))));
  g1(80, 60)=  g1(80, 60)+1.0;
  g1(86, 60)=  g1(86, 60)+1.0;
  g1(4, 55)=  g1(4, 55)+(-exp(y(55)));
  g1(7, 55)=  g1(7, 55)+(-exp(y(55)));
  g1(8, 55)=  g1(8, 55)+(-exp(y(55)));
  g1(23, 55)=  g1(23, 55)+1.0;
  g1(59, 55)=  g1(59, 55)+(-exp(y(55)));
  g1(88, 55)=  g1(88, 55)+(-1.0);
  g1(4, 73)=  g1(4, 73)+1.0;
  g1(8, 73)=  g1(8, 73)+1.0;
  g1(59, 73)=  g1(59, 73)+1.0;
  g1(60, 73)=  g1(60, 73)+1.0;
  g1(91, 73)=  g1(91, 73)+1.0/y(73);
  g1(5, 82)=  g1(5, 82)+1.0;
  g1(11, 82)=  g1(11, 82)+(-(1.0-y(77)-params(109)));
  g1(6, 3)=  g1(6, 3)+1.0;
  g1(7, 76)=  g1(7, 76)+1.0;
  g1(8, 75)=  g1(8, 75)+1.0;
  g1(9, 1)=  g1(9, 1)+1.0;
  g1(9, 41)=  g1(9, 41)+(-exp(y(41)));
  g1(59, 41)=  g1(59, 41)+exp(y(41));
  g1(60, 41)=  g1(60, 41)+(-(exp(y(41))*params(7)+exp(y(41))*params(8)));
  g1(10, 26)=  g1(10, 26)+1.0;
  g1(10, 46)=  g1(10, 46)+(-exp(y(46)));
  g1(18, 46)=  g1(18, 46)+1.0;
  g1(25, 46)=  g1(25, 46)+(-exp(y(46)));
  g1(28, 46)=  g1(28, 46)+(-1.0);
  g1(50, 46)=  g1(50, 46)+(-(T536*exp(y(91))*exp(y(46))*T518));
  g1(59, 46)=  g1(59, 46)+(-exp(y(46)));
  g1(11, 83)=  g1(11, 83)+1.0;
  g1(12, 42)=  g1(12, 42)+exp(y(42));
  g1(25, 42)=  g1(25, 42)+(-exp(y(42)));
  g1(29, 42)=  g1(29, 42)+(-1.0);
  g1(50, 42)=  g1(50, 42)+(-(T536*exp(y(91))*exp(y(42))*T518));
  g1(59, 42)=  g1(59, 42)+params(125)*exp(y(42));
  g1(13, 49)=  g1(13, 49)+exp(y(49));
  g1(17, 49)=  g1(17, 49)+(-exp(y(49)));
  g1(32, 49)=  g1(32, 49)+(-1.0);
  g1(79, 49)=  g1(79, 49)+1.0;
  g1(13, 18)=  g1(13, 18)+(-((1.0-params(107))*(-(params(31)/params(32)))));
  g1(21, 18)=  g1(21, 18)+(-1.0);
  g1(79, 18)=  g1(79, 18)+(-1.0);
  g1(94, 18)=  g1(94, 18)+(-1.0);
  g1(13, 18)=  g1(13, 18)+(-((1.0-params(107))*(params(31)*params(6))/params(32)));
  g1(13, 89)=  g1(13, 89)+(-((1.0-params(107))*(y(71)-1.0)*(params(32)*exp(y(89)))/(params(32)*params(32))));
  g1(68, 89)=  g1(68, 89)+1.0;
  g1(13, 71)=  g1(13, 71)+(-((1.0-params(107))*exp(y(89))/params(32)));
  g1(14, 71)=  g1(14, 71)+(-1.0);
  g1(13, 59)=  g1(13, 59)+(-(params(107)*params(62)*(1.0-params(114))*(exp(y(59))-(exp(y(60))*exp(y(59)))/T938)));
  g1(14, 59)=  g1(14, 59)+exp(y(59))*(1.0-params(114))*y(8)*(1.0-params(5));
  g1(32, 59)=  g1(32, 59)+1.0;
  g1(37, 59)=  g1(37, 59)+exp(y(59))*y(8)*(1.0-params(5));
  g1(76, 59)=  g1(76, 59)+1.0;
  g1(13, 38)=  g1(13, 38)+(-(params(107)*params(62)*(1.0-params(114))*(-params(103))));
  g1(15, 38)=  g1(15, 38)+(-1.0);
  g1(16, 38)=  g1(16, 38)+(-((-params(134))*(1.0+y(38))^((-params(134))-1.0)));
  g1(40, 38)=  g1(40, 38)+(-1.0);
  g1(42, 38)=  g1(42, 38)+1.0;
  g1(52, 38)=  g1(52, 38)+1.0;
  g1(13, 99)=  g1(13, 99)+(-(params(107)*params(62)*(1.0-params(114))*(-params(103))));
  g1(14, 99)=  g1(14, 99)+(-y(71));
  g1(73, 99)=  g1(73, 99)+1.0;
  g1(14, 8)=  g1(14, 8)+exp(y(59))*(1.0-params(114))*(1.0-params(5));
  g1(36, 8)=  g1(36, 8)+(-((1.0+y(51))*(exp(y(40))*params(5))/(exp(y(40))*exp(y(40)))));
  g1(37, 8)=  g1(37, 8)+exp(y(59))*(1.0-params(5));
  g1(38, 8)=  g1(38, 8)+1.0;
  g1(14, 72)=  g1(14, 72)+(-y(71));
  g1(36, 72)=  g1(36, 72)+(-((y(40)-y(40))*(-((1.0-params(113))*params(33)))/((1.0+y(72))*(1.0+y(72)))));
  g1(40, 72)=  g1(40, 72)+1.0;
  g1(53, 72)=  g1(53, 72)+(-y(2));
  g1(59, 72)=  g1(59, 72)+(-exp(y(41)));
  g1(14, 64)=  g1(14, 64)+y(71);
  g1(15, 64)=  g1(15, 64)+1.0;
  g1(14, 63)=  g1(14, 63)+(-y(71));
  g1(38, 63)=  g1(38, 63)+params(34)*params(6)*params(100);
  g1(40, 63)=  g1(40, 63)+1.0;
  g1(14, 65)=  g1(14, 65)+(-y(71));
  g1(14, 71)=  g1(14, 71)+(1.0-y(72)-params(14)-params(96)-y(99)+y(64))-y(63)-y(65);
  g1(14, 78)=  g1(14, 78)+(-(params(2)+params(3)*2.0*(y(78)-1.0)^(2.0-1.0)));
  g1(37, 78)=  g1(37, 78)+(-(params(2)+(y(78)-1.0)*params(3)*2.0+y(78)*params(3)*2.0));
  g1(41, 78)=  g1(41, 78)+(-((1.0-params(5))*1.0/y(78)));
  g1(81, 78)=  g1(81, 78)+1.0/y(78);
  g1(15, 31)=  g1(15, 31)+(-1.0);
  g1(16, 62)=  g1(16, 62)+1.0;
  g1(17, 23)=  g1(17, 23)+1.0;
  g1(32, 23)=  g1(32, 23)+1.0;
  g1(33, 23)=  g1(33, 23)+(-(1.0-params(5)));
  g1(76, 23)=  g1(76, 23)+1.0;
  g1(17, 25)=  g1(17, 25)+((1.0+y(25))-(params(14)+y(25)))/((1.0+y(25))*(1.0+y(25)));
  g1(53, 25)=  g1(53, 25)+y(2);
  g1(59, 25)=  g1(59, 25)+exp(y(41));
  g1(89, 25)=  g1(89, 25)+1.0;
  g1(92, 25)=  g1(92, 25)+(-1.0);
  g1(93, 25)=  g1(93, 25)+(-1.0);
  g1(94, 25)=  g1(94, 25)+(-1.0);
  g1(95, 25)=  g1(95, 25)+(-1.0);
  g1(96, 25)=  g1(96, 25)+(-1.0);
  g1(97, 25)=  g1(97, 25)+(-1.0);
  g1(18, 46)=  g1(18, 46)+(-params(41));
  g1(28, 46)=  g1(28, 46)+1.0;
  g1(18, 57)=  g1(18, 57)+(-params(28));
  g1(41, 57)=  g1(41, 57)+1.0;
  g1(44, 57)=  g1(44, 57)+1.0;
  g1(87, 57)=  g1(87, 57)+(-(params(126)*params(1)));
  g1(18, 57)=  g1(18, 57)+(-params(29));
  g1(19, 57)=  g1(19, 57)+(-params(52));
  g1(42, 57)=  g1(42, 57)+(-params(123));
  g1(44, 57)=  g1(44, 57)+(-1.0);
  g1(100, 57)=  g1(100, 57)+(-1.0);
  g1(19, 47)=  g1(19, 47)+1.0;
  g1(21, 47)=  g1(21, 47)+(-1.0);
  g1(25, 47)=  g1(25, 47)+(-exp(y(47)));
  g1(50, 47)=  g1(50, 47)+(-(T536*exp(y(91))*exp(y(47))*T518));
  g1(59, 47)=  g1(59, 47)+(-exp(y(47)));
  g1(19, 47)=  g1(19, 47)+(-params(48));
  g1(21, 47)=  g1(21, 47)+1.0;
  g1(19, 90)=  g1(19, 90)+(-1.0);
  g1(20, 90)=  g1(20, 90)+1.0;
  g1(20, 90)=  g1(20, 90)+(-params(71));
  g1(21, 19)=  g1(21, 19)+1.0;
  g1(21, 50)=  g1(21, 50)+1.0;
  g1(25, 50)=  g1(25, 50)+(-exp(y(50)));
  g1(32, 50)=  g1(32, 50)+1.0;
  g1(50, 50)=  g1(50, 50)+(-(T536*exp(y(91))*exp(y(50))*T518));
  g1(21, 50)=  g1(21, 50)+(-1.0);
  g1(22, 88)=  g1(22, 88)+1.0;
  g1(23, 4)=  g1(23, 4)+(-params(120));
  g1(90, 4)=  g1(90, 4)+1.0;
  g1(23, 85)=  g1(23, 85)+(-1.0);
  g1(24, 85)=  g1(24, 85)+1.0;
  g1(24, 85)=  g1(24, 85)+(-params(122));
  g1(25, 45)=  g1(25, 45)+(-exp(y(45)));
  g1(26, 45)=  g1(26, 45)+(-exp(y(45)));
  g1(51, 45)=  g1(51, 45)+exp(y(45));
  g1(53, 45)=  g1(53, 45)+(-exp(y(45)));
  g1(54, 45)=  g1(54, 45)+1.0;
  g1(25, 48)=  g1(25, 48)+exp(y(48));
  g1(26, 48)=  g1(26, 48)+exp(y(48));
  g1(50, 48)=  g1(50, 48)+exp(y(48));
  g1(53, 48)=  g1(53, 48)+exp(y(48));
  g1(55, 48)=  g1(55, 48)+1.0;
  g1(26, 74)=  g1(26, 74)+1.0;
  g1(27, 24)=  g1(27, 24)+1.0;
  g1(33, 24)=  g1(33, 24)+(-(params(5)*(1.0+params(58))));
  g1(28, 16)=  g1(28, 16)+1.0;
  g1(95, 16)=  g1(95, 16)+(-1.0);
  g1(29, 9)=  g1(29, 9)+1.0;
  g1(93, 9)=  g1(93, 9)+(-1.0);
  g1(29, 42)=  g1(29, 42)+1.0;
  g1(30, 11)=  g1(30, 11)+1.0;
  g1(30, 43)=  g1(30, 43)+1.0;
  g1(31, 12)=  g1(31, 12)+1.0;
  g1(33, 32)=  g1(33, 32)+(-(1.0-params(5)));
  g1(81, 32)=  g1(81, 32)+(-1.0);
  g1(85, 32)=  g1(85, 32)+(-(1.0-params(5)));
  g1(33, 28)=  g1(33, 28)+(-params(5));
  g1(34, 28)=  g1(34, 28)+1.0;
  g1(85, 28)=  g1(85, 28)+(-params(5));
  g1(34, 100)=  g1(34, 100)+(-1.0);
  g1(35, 100)=  g1(35, 100)+1.0;
  g1(35, 100)=  g1(35, 100)+(-params(83));
  g1(36, 101)=  g1(36, 101)+(exp(y(60))*exp(y(101)))/T938;
  g1(74, 101)=  g1(74, 101)+1.0;
  g1(36, 51)=  g1(36, 51)+(-((y(8)*params(5))/exp(y(40))));
  g1(70, 51)=  g1(70, 51)+1.0;
  g1(36, 40)=  g1(36, 40)+(-(((1.0-params(113))*params(33))/(1.0+y(72))));
  g1(38, 86)=  g1(38, 86)+1.0;
  g1(66, 86)=  g1(66, 86)+1.0;
  g1(38, 63)=  g1(38, 63)+params(34)*params(6)*(1.0-params(100));
  g1(39, 80)=  g1(39, 80)+1.0;
  g1(45, 80)=  g1(45, 80)+(-(T448*params(36)/params(115)));
  g1(80, 80)=  g1(80, 80)+1.0;
  g1(39, 81)=  g1(39, 81)+(-1.0);
  g1(42, 38)=  g1(42, 38)+(-params(49));
  g1(42, 102)=  g1(42, 102)+(-(1.0-params(49)+(-params(116))));
  g1(75, 102)=  g1(75, 102)+1.0;
  g1(42, 58)=  g1(42, 58)+(-params(124));
  g1(42, 93)=  g1(42, 93)+(-1.0);
  g1(43, 93)=  g1(43, 93)+1.0;
  g1(43, 93)=  g1(43, 93)+(-params(75));
  g1(44, 36)=  g1(44, 36)+1.0;
  g1(45, 56)=  g1(45, 56)+((params(102)*y(80)+(1.0-params(102))*y(80))-params(44)-params(37))*(params(115)*params(36)*params(6)*exp(y(56)))/(params(115)*params(115));
  g1(45, 80)=  g1(45, 80)+(params(36)*params(6)*exp(y(56)))/params(115)*params(102);
  g1(45, 80)=  g1(45, 80)+(params(36)*params(6)*exp(y(56)))/params(115)*(1.0-params(102));
  g1(46, 52)=  g1(46, 52)+1.0;
  g1(47, 52)=  g1(47, 52)+1.0;
  g1(50, 52)=  g1(50, 52)+(-(T536*exp(y(91))*(exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*(1.0-params(97))*(exp(y(52))*exp(y(53)))/(exp(y(53))*exp(y(53)))*(params(105)-1.0)*(exp(y(52))/exp(y(53)))^(params(105)-1.0-1.0)));
  g1(58, 52)=  g1(58, 52)+(-(exp(y(97))*(1.0-x(4)+params(35)*(params(37)+params(6)*((params(101)*y(69)+(1.0-params(101))*y(69))-params(37))-y(69)))*exp(y(52))*(1.0-params(111))*exp(y(52))^(1.0-params(111)-1.0)));
  g1(46, 53)=  g1(46, 53)+(-(1.0-params(97)-(2.0*(params(105)-1.0)*params(97)*(1.0-params(97))*2.0*y(53)^(2.0-1.0))/(2.0*2.0)));
  g1(48, 53)=  g1(48, 53)+1.0;
  g1(50, 53)=  g1(50, 53)+(-(T536*exp(y(91))*(exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*(1.0-params(97))*(params(105)-1.0)*(exp(y(52))/exp(y(53)))^(params(105)-1.0-1.0)*(-(exp(y(52))*exp(y(53))))/(exp(y(53))*exp(y(53)))));
  g1(56, 53)=  g1(56, 53)+exp(y(53));
  g1(46, 94)=  g1(46, 94)+(-1.0);
  g1(82, 94)=  g1(82, 94)+1.0;
  g1(47, 52)=  g1(47, 52)+(-1.0);
  g1(48, 53)=  g1(48, 53)+(-1.0);
  g1(56, 53)=  g1(56, 53)+(-(exp(y(96))*T608*(1.0-params(108))*(exp(y(53))/(1.0+y(63)))^(1.0-params(108)-1.0)*((1.0+y(63))*exp(y(53)))/((1.0+y(63))*(1.0+y(63)))));
  g1(48, 66)=  g1(48, 66)+(-1.0);
  g1(55, 66)=  g1(55, 66)+(-1.0);
  g1(98, 66)=  g1(98, 66)+(-1.0);
  g1(49, 54)=  g1(49, 54)+1.0;
  g1(51, 54)=  g1(51, 54)+(-(T572*T557*(exp(y(54))*(1.0-params(110))*exp(y(87))*T550^params(104)+exp(y(54))*(1.0-params(110))*exp(y(87))*(-(exp(y(7))^(params(110)*params(4))*exp(y(54))))/(exp(y(54))*exp(y(54)))*params(104)*T550^(params(104)-1.0))));
  g1(58, 54)=  g1(58, 54)+exp(y(54));
  g1(49, 54)=  g1(49, 54)+(-1.0);
  g1(49, 69)=  g1(49, 69)+(-1.0);
  g1(54, 69)=  g1(54, 69)+(-1.0);
  g1(58, 69)=  g1(58, 69)+(-(exp(y(97))*T639*(-params(35))));
  g1(99, 69)=  g1(99, 69)+(-1.0);
  g1(50, 91)=  g1(50, 91)+(-((exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*T518*exp(y(91))*T536));
  g1(69, 91)=  g1(69, 91)+1.0;
  g1(50, 21)=  g1(50, 21)+(-((exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*T518*exp(y(91))*(T532*params(99)*((1.0+y(21))-params(44))^(params(99)-1.0))/T952*params(6)*params(50)*T533^(params(6)*params(50)-1.0)));
  g1(50, 21)=  g1(50, 21)+(-((exp(y(46))+exp(y(42))+exp(y(50))+exp(y(47)))*T518*exp(y(91))*params(6)*params(50)*T533^(params(6)*params(50)-1.0)*(-(((1.0+y(21))-params(44))^params(99)*(1.0-params(99))*((1.0+y(21))-params(44))^(1.0-params(99)-1.0)))/T952));
  g1(55, 21)=  g1(55, 21)+(-1.0);
  g1(97, 21)=  g1(97, 21)+(-1.0);
  g1(51, 87)=  g1(51, 87)+(-(T558*T572));
  g1(67, 87)=  g1(67, 87)+1.0;
  g1(51, 7)=  g1(51, 7)+(-(T572*T557*exp(y(54))*(1.0-params(110))*exp(y(87))*params(104)*T550^(params(104)-1.0)*(exp(y(54))*exp(y(7))*params(110)*params(4)*exp(y(7))^(params(110)*params(4)-1.0))/(exp(y(54))*exp(y(54)))));
  g1(56, 7)=  g1(56, 7)+(-(exp(y(96))*T606*exp(y(7))*params(4)*exp(y(7))^(params(4)-1.0)*params(108)*exp(y(7))^params(4)^(params(108)-1.0)));
  g1(57, 7)=  g1(57, 7)+1.0;
  g1(51, 61)=  g1(51, 61)+(-(T572*exp(y(54))*(1.0-params(110))*exp(y(87))*T550^params(104)*exp(y(61))*params(4)*exp(y(61))^(params(4)-1.0)));
  g1(63, 61)=  g1(63, 61)+(-params(94));
  g1(64, 61)=  g1(64, 61)+1.0;
  g1(51, 14)=  g1(51, 14)+(-(T558*(T568*params(98)*((1.0+y(14))-params(44))^(params(98)-1.0))/T954*params(6)*params(26)*T569^(params(6)*params(26)-1.0)));
  g1(51, 14)=  g1(51, 14)+(-(T558*params(6)*params(26)*T569^(params(6)*params(26)-1.0)*(-(((1.0+y(14))-params(44))^params(98)*(1.0-params(98))*((1.0+y(14))-params(44))^(1.0-params(98)-1.0)))/T954));
  g1(54, 14)=  g1(54, 14)+(-1.0);
  g1(96, 14)=  g1(96, 14)+(-1.0);
  g1(52, 39)=  g1(52, 39)+(-1.0);
  g1(61, 39)=  g1(61, 39)+1.0;
  g1(52, 13)=  g1(52, 13)+(-1.0);
  g1(52, 2)=  g1(52, 2)+params(95);
  g1(53, 2)=  g1(53, 2)+1.0;
  g1(52, 98)=  g1(52, 98)+(-1.0);
  g1(72, 98)=  g1(72, 98)+1.0;
  g1(53, 2)=  g1(53, 2)+(-((1.0+y(72))-y(34)-y(25)));
  g1(54, 45)=  g1(54, 45)+(-1.0);
  g1(55, 48)=  g1(55, 48)+(-1.0);
  g1(56, 96)=  g1(56, 96)+(-(T606*T608*exp(y(96))));
  g1(83, 96)=  g1(83, 96)+1.0;
  g1(57, 7)=  g1(57, 7)+(-1.0);
  g1(57, 68)=  g1(57, 68)+(-1.0);
  g1(62, 68)=  g1(62, 68)+1.0;
  g1(57, 13)=  g1(57, 13)+(-1.0);
  g1(58, 69)=  g1(58, 69)+(-(exp(y(97))*T639*params(35)*params(6)*params(101)));
  g1(58, 69)=  g1(58, 69)+(-(exp(y(97))*T639*params(35)*params(6)*(1.0-params(101))));
  g1(58, 97)=  g1(58, 97)+(-(exp(y(97))*(1.0-x(4)+params(35)*(params(37)+params(6)*((params(101)*y(69)+(1.0-params(101))*y(69))-params(37))-y(69)))*T639));
  g1(84, 97)=  g1(84, 97)+1.0;
  g1(59, 41)=  g1(59, 41)+(-(((1.0+y(72))-y(34)-y(25))*exp(y(41))));
  g1(60, 41)=  g1(60, 41)+(-(params(8)*(-exp(y(41)))));
  g1(60, 73)=  g1(60, 73)+(-1.0);
  g1(91, 73)=  g1(91, 73)+(-(1.0/y(73)));
  g1(61, 39)=  g1(61, 39)+(-params(85));
  g1(62, 39)=  g1(62, 39)+(-params(88));
  g1(63, 39)=  g1(63, 39)+(-params(91));
  g1(61, 68)=  g1(61, 68)+(-params(86));
  g1(62, 68)=  g1(62, 68)+(-params(89));
  g1(63, 68)=  g1(63, 68)+(-params(92));
  g1(61, 37)=  g1(61, 37)+(-params(87));
  g1(62, 37)=  g1(62, 37)+(-params(90));
  g1(63, 37)=  g1(63, 37)+(-params(93));
  g1(63, 37)=  g1(63, 37)+1.0;
  g1(64, 37)=  g1(64, 37)+(-1.0);
  g1(64, 61)=  g1(64, 61)+(-1.0);
  g1(65, 84)=  g1(65, 84)+(-params(63));
  g1(66, 86)=  g1(66, 86)+(-params(67));
  g1(67, 87)=  g1(67, 87)+(-params(68));
  g1(68, 89)=  g1(68, 89)+(-params(70));
  g1(69, 91)=  g1(69, 91)+(-params(72));
  g1(70, 51)=  g1(70, 51)+(-params(74));
  g1(71, 92)=  g1(71, 92)+(-params(73));
  g1(72, 98)=  g1(72, 98)+(-params(81));
  g1(73, 99)=  g1(73, 99)+(-params(82));
  g1(75, 102)=  g1(75, 102)+(-params(51));
  g1(76, 59)=  g1(76, 59)+(-1.0);
  g1(76, 65)=  g1(76, 65)+(-1.0);
  g1(77, 65)=  g1(77, 65)+1.0;
  g1(77, 95)=  g1(77, 95)+(-1.0);
  g1(78, 95)=  g1(78, 95)+1.0;
  g1(78, 95)=  g1(78, 95)+(-params(77));
  g1(79, 49)=  g1(79, 49)+(-1.0);
  g1(79, 23)=  g1(79, 23)+1.0;
  g1(80, 60)=  g1(80, 60)+(-1.0);
  g1(86, 60)=  g1(86, 60)+(-1.0);
  g1(81, 78)=  g1(81, 78)+(-(1.0/y(78)));
  g1(82, 94)=  g1(82, 94)+(-params(76));
  g1(83, 96)=  g1(83, 96)+(-params(78));
  g1(84, 97)=  g1(84, 97)+(-params(80));
  g1(85, 29)=  g1(85, 29)+1.0;
  g1(86, 33)=  g1(86, 33)+1.0;
  g1(88, 30)=  g1(88, 30)+1.0;
  g1(88, 55)=  g1(88, 55)+1.0;
  g1(89, 25)=  g1(89, 25)+(-params(79));
  g1(90, 4)=  g1(90, 4)+(-params(64));
  g1(101, 4)=  g1(101, 4)+(-1.0);
  g1(90, 5)=  g1(90, 5)+(-params(65));
  g1(102, 5)=  g1(102, 5)+(-1.0);
  g1(90, 6)=  g1(90, 6)+(-params(66));
  g1(91, 27)=  g1(91, 27)+(-1.0);
  g1(92, 35)=  g1(92, 35)+1.0;
  g1(93, 10)=  g1(93, 10)+1.0;
  g1(94, 20)=  g1(94, 20)+1.0;
  g1(95, 17)=  g1(95, 17)+1.0;
  g1(96, 15)=  g1(96, 15)+1.0;
  g1(97, 22)=  g1(97, 22)+1.0;
  g1(98, 67)=  g1(98, 67)+1.0;
  g1(99, 70)=  g1(99, 70)+1.0;
  g1(100, 58)=  g1(100, 58)+1.0;
  g1(101, 5)=  g1(101, 5)+1.0;
  g1(102, 6)=  g1(102, 6)+1.0;
  if ~isreal(g1)
    g1 = real(g1)+2*imag(g1);
  end
end