Jac_Sarco.m

Go to the documentation of this file.

#include "SHDynamics.m"
namespace models{
namespace motorunit{


/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

function J = models.motorunit.SHDynamics.Jac_Sarco(matrix<double> y,rowvec<double> t) {


c = this.sarcoconst;

c36c37 = c(36)*c(37);

/*  Conditional statements
 * See dytd_sarcomere; */
alg72 = y(55)*0.001*y(32) - c(98);
h73 = alg72 > 0;
h74 = alg72 <= 0;
h75 = y(1)*((y(4) -  y(5)*(exp(((-c(7)*y(1))/c36c37))))/(1 - (exp(((-c(7)*y(1))/c36c37))))) > 0;

J = sparse(56,56);
J(1,1) = ((c(38)*(313/(2*(y(4)/y(5) + 5)) + exp((c(7)*y(1))/(c(36)*c(37)))*(1565/(2*(y(4)/y(5) + 5)) - 313/2)))/(45*(exp(-(c(24) - y(1))/c(27)) + 1)^4*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)) - 1000/c(3) + (c(39)*y(9)^4*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(40)*y(11)^3*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(48)*((3*c(7)*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (4*c(38)*y(1)*exp(-(c(24) - y(1))/c(27))*(313/(2*(y(4)/y(5) + 5)) + exp((c(7)*y(1))/(c(36)*c(37)))*(1565/(2*(y(4)/y(5) + 5)) - 313/2)))/(45*c(27)*(exp(-(c(24) - y(1))/c(27)) + 1)^5*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(38)*y(1)*exp((c(7)*y(1))/(c(36)*c(37)))*(1565/(2*(y(4)/y(5) + 5)) - 313/2))/(45*c(36)*c(37)*(exp(-(c(24) - y(1))/c(27)) + 1)^4*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(38)*y(1)*exp((c(7)*y(1))/(c(36)*c(37)))*(313/(2*(y(4)/y(5) + 5)) + exp((c(7)*y(1))/(c(36)*c(37)))*(1565/(2*(y(4)/y(5) + 5)) - 313/2)))/(45*c(36)*c(37)*(exp(-(c(24) - y(1))/c(27)) + 1)^4*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(39)*y(1)*y(5)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(40)*y(1)*y(8)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) - (c(7)*c(41)*h75*y(1)*y(5)^3*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(1,2) = 1000/(c(1)*c(3));
J(1,4) = ((c(39)*y(1)*y(9)^4*y(10))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(38)*y(1)*(313/(2*y(5)*(y(4)/y(5) + 5)^2) + (1565*exp((c(7)*y(1))/(c(36)*c(37))))/(2*y(5)*(y(4)/y(5) + 5)^2)))/(45*(exp(-(c(24) - y(1))/c(27)) + 1)^4*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(44)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(1,5) = ((c(38)*y(1)*((313*y(4))/(2*y(5)^2*(y(4)/y(5) + 5)^2) + (1565*y(4)*exp((c(7)*y(1))/(c(36)*c(37))))/(2*y(5)^2*(y(4)/y(5) + 5)^2)))/(45*(exp(-(c(24) - y(1))/c(27)) + 1)^4*(exp((c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (2*c(7)*c(45)*c(48))/(y(5)^2*(c(45)/y(5) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)) - (c(41)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37))))*(2*y(5)*exp(-2*c(42)*log(y(5)/y(4))) - 2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43) - (2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(1,6) = ((c(40)*y(1)*y(11)^3*y(12)*y(13))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (3*c(7)*c(46)*c(48))/(y(6)^2*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)))/c(1);
J(1,8) = ((2*c(7)*c(48)*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp((10*y(8))/673))/(23555*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)))/c(1);
J(1,9) = (2*c(39)*y(1)*y(9)^3*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*c(1)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(1,10) = (c(39)*y(1)*y(9)^4*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(1)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(1,11) = (c(40)*y(1)*y(11)^2*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*c(1)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(1,12) = (c(40)*y(1)*y(11)^3*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*c(1)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(1,13) = (c(40)*y(1)*y(11)^3*y(12)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*c(1)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(2,1) = 1000/(c(1)*c(2)*c(3));
J(2,2) = ((c(38)*c(50)*(313/(2*(y(4)/y(3) + 5)) + exp((c(7)*y(2))/(c(36)*c(37)))*(1565/(2*(y(4)/y(3) + 5)) - 313/2)))/(45*(exp(-(c(24) - y(2))/c(27)) + 1)^4*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)) - 1000/(c(2)*c(3)) + (c(39)*c(52)*y(14)^4*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(48)*c(54)*((3*c(7)*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) + (c(40)*c(53)*y(16)^3*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(51)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (4*c(38)*c(50)*y(2)*exp(-(c(24) - y(2))/c(27))*(313/(2*(y(4)/y(3) + 5)) + exp((c(7)*y(2))/(c(36)*c(37)))*(1565/(2*(y(4)/y(3) + 5)) - 313/2)))/(45*c(27)*(exp(-(c(24) - y(2))/c(27)) + 1)^5*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(38)*c(50)*y(2)*exp((c(7)*y(2))/(c(36)*c(37)))*(313/(2*(y(4)/y(3) + 5)) + exp((c(7)*y(2))/(c(36)*c(37)))*(1565/(2*(y(4)/y(3) + 5)) - 313/2)))/(45*c(36)*c(37)*(exp(-(c(24) - y(2))/c(27)) + 1)^4*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(38)*c(50)*y(2)*exp((c(7)*y(2))/(c(36)*c(37)))*(1565/(2*(y(4)/y(3) + 5)) - 313/2))/(45*c(36)*c(37)*(exp(-(c(24) - y(2))/c(27)) + 1)^4*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(3)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(7)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*c(51)*y(2)*y(3)^3*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(2,3) = ((c(38)*c(50)*y(2)*((313*y(4))/(2*y(3)^2*(y(4)/y(3) + 5)^2) + (1565*y(4)*exp((c(7)*y(2))/(c(36)*c(37))))/(2*y(3)^2*(y(4)/y(3) + 5)^2)))/(45*(exp(-(c(24) - y(2))/c(27)) + 1)^4*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (2*c(7)*c(45)*c(48)*c(54))/(y(3)^2*(c(45)/y(3) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)) - (c(41)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37))))*(2*y(3)*exp(-2*c(42)*log(y(3)/y(4))) - 2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43) - (2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(2,4) = -((c(38)*c(50)*y(2)*(313/(2*y(3)*(y(4)/y(3) + 5)^2) + (1565*exp((c(7)*y(2))/(c(36)*c(37))))/(2*y(3)*(y(4)/y(3) + 5)^2)))/(45*(exp(-(c(24) - y(2))/c(27)) + 1)^4*(exp((c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(39)*c(52)*y(2)*y(14)^4*y(15))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*c(44)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/c(1);
J(2,6) = ((c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (3*c(7)*c(46)*c(48)*c(54))/(y(6)^2*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)))/c(1);
J(2,7) = ((2*c(7)*c(48)*c(54)*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp((10*y(7))/673))/(23555*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) - (c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)))/c(1);
J(2,14) = (2*c(39)*c(52)*y(2)*y(14)^3*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*c(1)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(2,15) = (c(39)*c(52)*y(2)*y(14)^4*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(1)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(2,16) = (c(40)*c(53)*y(2)*y(16)^2*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*c(1)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(2,17) = (c(40)*c(53)*y(2)*y(16)^3*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*c(1)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(2,18) = (c(40)*c(53)*y(2)*y(16)^3*y(17)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*c(1)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(3,2) = -((c(39)*c(52)*y(14)^4*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (2*c(7)*c(48)*c(54)*((3*c(7)*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) - (c(41)*c(51)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(3)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) - (c(7)*c(41)*c(51)*y(2)*y(3)^3*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(4)*c(7));
J(3,3) = - 1/c(8) - ((4*c(7)*c(45)*c(48)*c(54))/(y(3)^2*(c(45)/y(3) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)) - (c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37))))*(2*y(3)*exp(-2*c(42)*log(y(3)/y(4))) - 2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43) - (2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(4)*c(7));
J(3,4) = -((c(39)*c(52)*y(2)*y(14)^4*y(15))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(44)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(4)*c(7));
J(3,5) = 1/c(8);
J(3,6) = -(3*c(46)*c(48)*c(54))/(500*c(4)*y(6)^2*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1));
J(3,7) = (c(48)*c(54)*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp((10*y(7))/673))/(5888750*c(4)*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2);
J(3,14) = -(c(39)*c(52)*y(2)*y(14)^3*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(12500*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(3,15) = -(c(39)*c(52)*y(2)*y(14)^4*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(4,1) = ((c(39)*y(9)^4*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (2*c(7)*c(48)*((3*c(7)*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (c(41)*h75*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(39)*y(1)*y(5)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*h75*y(1)*y(5)^3*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(5)*c(7));
J(4,2) = (c(10)*((c(39)*c(52)*y(14)^4*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (2*c(7)*c(48)*c(54)*((3*c(7)*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) - (c(41)*c(51)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(3)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) - (c(7)*c(41)*c(51)*y(2)*y(3)^3*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2)))/(1000*c(4)*c(7));
J(4,3) = (c(10)*((4*c(7)*c(45)*c(48)*c(54))/(y(3)^2*(c(45)/y(3) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)) - (c(39)*c(52)*y(2)*y(14)^4*y(15)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37))))*(2*y(3)*exp(-2*c(42)*log(y(3)/y(4))) - 2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43) - (2*c(42)*y(3)*exp(-2*c(42)*log(y(3)/y(4))))/c(43))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2)))/(1000*c(4)*c(7));
J(4,4) = ((c(39)*y(1)*y(9)^4*y(10))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(44)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(5)*c(7)) + (c(10)*((c(39)*c(52)*y(2)*y(14)^4*y(15))/(50*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*c(51)*y(2)*y(3)^2*exp(-2*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*(1/((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))^2*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(44)*c(51)*y(2)*y(3)^4*exp(-4*c(42)*log(y(3)/y(4)))*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(2)*(2*c(42) - 2))/(c(36)*c(37)))*((y(3)^2*exp(-2*c(42)*log(y(3)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2)))/(1000*c(4)*c(7));
J(4,5) = ((4*c(7)*c(45)*c(48))/(y(5)^2*(c(45)/y(5) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)) - (c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37))))*(2*y(5)*exp(-2*c(42)*log(y(5)/y(4))) - 2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43) - (2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(5)*c(7));
J(4,6) = (3*c(46)*c(48))/(500*c(5)*y(6)^2*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)) + (3*c(10)*c(46)*c(48)*c(54))/(500*c(4)*y(6)^2*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1));
J(4,7) = -(c(10)*c(48)*c(54)*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp((10*y(7))/673))/(5888750*c(4)*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2);
J(4,8) = -(c(48)*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp((10*y(8))/673))/(5888750*c(5)*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2);
J(4,9) = (c(39)*y(1)*y(9)^3*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(12500*c(5)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(4,10) = (c(39)*y(1)*y(9)^4*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50000*c(5)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(4,14) = (c(10)*c(39)*c(52)*y(2)*y(14)^3*y(15)*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(12500*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(4,15) = (c(10)*c(39)*c(52)*y(2)*y(14)^4*(y(4) - y(3)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(50000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(5,1) = -((c(39)*y(9)^4*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (2*c(7)*c(48)*((3*c(7)*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (c(41)*h75*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(39)*y(1)*y(5)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*h75*y(1)*y(5)^3*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(7)*c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2) + (c(7)*c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(2*c(42) - 2)*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(36)*c(37)*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(6)*c(7));
J(5,3) = 1/c(11);
J(5,4) = -((c(39)*y(1)*y(9)^4*y(10))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*c(42)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*c(44)*h75*y(1)*y(5)^4*exp(-4*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*c(43)*c(47)^2*y(4)*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(6)*c(7));
J(5,5) = - 1/c(11) - ((4*c(7)*c(45)*c(48))/(y(5)^2*(c(45)/y(5) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)) - (c(39)*y(1)*y(9)^4*y(10)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(50*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (c(41)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37))))*(2*y(5)*exp(-2*c(42)*log(y(5)/y(4))) - 2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))))/(50*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))^2*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(42)*h75*y(1)*y(5)*exp(-2*c(42)*log(y(5)/y(4)))*(1/((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1) - 1)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(41)*c(44)*h75*y(1)*y(5)^2*exp(-2*c(42)*log(y(5)/y(4)))*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((2*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43) - (2*c(42)*y(5)*exp(-2*c(42)*log(y(5)/y(4))))/c(43))*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50*c(47)^2*(c(43) + y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)*((c(44)*exp((c(7)*y(1)*(2*c(42) - 2))/(c(36)*c(37)))*((y(5)^2*exp(-2*c(42)*log(y(5)/y(4))))/c(43) + 1))/c(47)^2 + 1)^2))/(1000*c(6)*c(7));
J(5,6) = -(3*c(46)*c(48))/(500*c(6)*y(6)^2*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1));
J(5,8) = (c(48)*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp((10*y(8))/673))/(5888750*c(6)*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2);
J(5,9) = -(c(39)*y(1)*y(9)^3*y(10)*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(12500*c(6)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(5,10) = -(c(39)*y(1)*y(9)^4*(y(4) - y(5)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(50000*c(6)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(6,1) = ((c(40)*y(11)^3*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (3*c(7)*c(48)*((3*c(7)*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) + (c(7)*c(40)*y(1)*y(8)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2))/(1000*c(5)*c(7));
J(6,2) = (c(10)*((c(40)*c(53)*y(16)^3*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (3*c(7)*c(48)*c(54)*((3*c(7)*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(7)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1))))/(1000*c(4)*c(7));
J(6,3) = -(3*c(10)*c(45)*c(48)*c(54))/(500*c(4)*y(3)^2*(c(45)/y(3) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1));
J(6,5) = -(3*c(45)*c(48))/(500*c(5)*y(5)^2*(c(45)/y(5) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1));
J(6,6) = ((c(40)*y(1)*y(11)^3*y(12)*y(13))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (9*c(7)*c(46)*c(48))/(y(6)^2*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)))/(1000*c(5)*c(7)) + (c(10)*((c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (9*c(7)*c(46)*c(48)*c(54))/(y(6)^2*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1))))/(1000*c(4)*c(7));
J(6,7) = (c(10)*((6*c(7)*c(48)*c(54)*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp((10*y(7))/673))/(23555*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) - (c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1))))/(1000*c(4)*c(7));
J(6,8) = ((6*c(7)*c(48)*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp((10*y(8))/673))/(23555*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)))/(1000*c(5)*c(7));
J(6,11) = (c(40)*y(1)*y(11)^2*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25000*c(5)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(6,12) = (c(40)*y(1)*y(11)^3*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75000*c(5)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(6,13) = (c(40)*y(1)*y(11)^3*y(12)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75000*c(5)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(6,16) = (c(10)*c(40)*c(53)*y(2)*y(16)^2*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(6,17) = (c(10)*c(40)*c(53)*y(2)*y(16)^3*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(6,18) = (c(10)*c(40)*c(53)*y(2)*y(16)^3*y(17)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(7,2) = -((c(40)*c(53)*y(16)^3*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (3*c(7)*c(48)*c(54)*((3*c(7)*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37)))*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)^2) + (c(7)*c(40)*c(53)*y(2)*y(7)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)))/(1000*c(4)*c(7));
J(7,3) = (3*c(45)*c(48)*c(54))/(500*c(4)*y(3)^2*(c(45)/y(3) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1));
J(7,6) = -((c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)) - (9*c(7)*c(46)*c(48)*c(54))/(y(6)^2*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)))/(1000*c(4)*c(7));
J(7,7) = - 1/c(9) - ((6*c(7)*c(48)*c(54)*exp(-(c(7)*y(2))/(c(36)*c(37)))*exp((10*y(7))/673))/(23555*(c(45)/y(3) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(2))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(2))/(c(36)*c(37)))*(exp((10*y(7))/673)/175 - 1/175) + 1)^2) - (c(40)*c(53)*y(2)*y(16)^3*y(17)*y(18)*exp(-(c(7)*y(2))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1)))/(1000*c(4)*c(7));
J(7,8) = 1/c(9);
J(7,16) = -(c(40)*c(53)*y(2)*y(16)^2*y(17)*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(25000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(7,17) = -(c(40)*c(53)*y(2)*y(16)^3*y(18)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(7,18) = -(c(40)*c(53)*y(2)*y(16)^3*y(17)*(y(6) - y(7)*exp(-(c(7)*y(2))/(c(36)*c(37)))))/(75000*c(4)*c(7)*(exp(-(c(7)*y(2))/(c(36)*c(37))) - 1));
J(8,1) = -((c(40)*y(11)^3*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (3*c(7)*c(48)*((3*c(7)*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/(250*c(36)*c(37)) + (c(7)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175))/(c(36)*c(37))))/((c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) + (c(7)*c(40)*y(1)*y(8)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) + (c(7)*c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37)))*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75*c(36)*c(37)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)^2))/(1000*c(6)*c(7));
J(8,5) = (3*c(45)*c(48))/(500*c(6)*y(5)^2*(c(45)/y(5) + 1)^3*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1));
J(8,6) = -((c(40)*y(1)*y(11)^3*y(12)*y(13))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)) - (9*c(7)*c(46)*c(48))/(y(6)^2*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^4*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)))/(1000*c(6)*c(7));
J(8,7) = 1/c(12);
J(8,8) = - 1/c(12) - ((6*c(7)*c(48)*exp(-(c(7)*y(1))/(c(36)*c(37)))*exp((10*y(8))/673))/(23555*(c(45)/y(5) + 1)^2*(c(46)/y(6) + 1)^3*((3*exp(-(c(7)*y(1))/(10*c(36)*c(37))))/25 + exp(-(c(7)*y(1))/(c(36)*c(37)))*(exp((10*y(8))/673)/175 - 1/175) + 1)^2) - (c(40)*y(1)*y(11)^3*y(12)*y(13)*exp(-(c(7)*y(1))/(c(36)*c(37))))/(75*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1)))/(1000*c(6)*c(7));
J(8,11) = -(c(40)*y(1)*y(11)^2*y(12)*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(25000*c(6)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(8,12) = -(c(40)*y(1)*y(11)^3*y(13)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75000*c(6)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(8,13) = -(c(40)*y(1)*y(11)^3*y(12)*(y(6) - y(8)*exp(-(c(7)*y(1))/(c(36)*c(37)))))/(75000*c(6)*c(7)*(exp(-(c(7)*y(1))/(c(36)*c(37))) - 1));
J(9,1) = (c(17)*(y(9) - 1))/(exp((c(22) - y(1))/c(33)) - 1) + (c(20)*y(9)*exp((c(22) - y(1))/c(35)))/c(35) - (c(17)*exp((c(22) - y(1))/c(33))*(c(22) - y(1))*(y(9) - 1))/(c(33)*(exp((c(22) - y(1))/c(33)) - 1)^2);
J(9,9) = - c(20)*exp((c(22) - y(1))/c(35)) - (c(17)*(c(22) - y(1)))/(exp((c(22) - y(1))/c(33)) - 1);
J(10,1) = - (exp((4*y(1))/103 + 160/103)*(y(10) - 1/(exp(-(c(26) - y(1))/c(29)) + 1)))/25750 - (exp(-(c(26) - y(1))/c(29))*exp((4*y(1))/103 + 160/103))/(1000*c(29)*(exp(-(c(26) - y(1))/c(29)) + 1)^2);
J(10,10) = -exp((4*y(1))/103 + 160/103)/1000;
J(11,1) = (c(16)*(y(11) - 1))/(exp((c(21) - y(1))/c(32)) - 1) + (c(19)*y(11)*exp((c(21) - y(1))/c(34)))/c(34) - (c(16)*exp((c(21) - y(1))/c(32))*(c(21) - y(1))*(y(11) - 1))/(c(32)*(exp((c(21) - y(1))/c(32)) - 1)^2);
J(11,11) = - c(19)*exp((c(21) - y(1))/c(34)) - (c(16)*(c(21) - y(1)))/(exp((c(21) - y(1))/c(32)) - 1);
J(12,1) = (c(15)*exp((c(23) - y(1))/c(30))*(y(12) - 1))/c(30) - (c(18)*y(12)*exp((c(23) - y(1))/c(31)))/(c(31)*(exp((c(23) - y(1))/c(31)) + 1)^2);
J(12,12) = - c(18)/(exp((c(23) - y(1))/c(31)) + 1) - c(15)*exp((c(23) - y(1))/c(30));
J(13,1) = - (y(13) - 1/(exp(-(c(25) - y(1))/c(28)) + 1))*((113*c(49))/857100000 + (113*y(1))/857100000) - (exp(-(c(25) - y(1))/c(28))*((113*(c(49)/100 + y(1)/100)^2)/171420 + 1/42855))/(c(28)*(exp(-(c(25) - y(1))/c(28)) + 1)^2);
J(13,13) = - (113*(c(49)/100 + y(1)/100)^2)/171420 - 1/42855;
J(14,2) = (c(17)*(y(14) - 1))/(exp((c(22) - y(2))/c(33)) - 1) + (c(20)*y(14)*exp((c(22) - y(2))/c(35)))/c(35) - (c(17)*exp((c(22) - y(2))/c(33))*(c(22) - y(2))*(y(14) - 1))/(c(33)*(exp((c(22) - y(2))/c(33)) - 1)^2);
J(14,14) = - c(20)*exp((c(22) - y(2))/c(35)) - (c(17)*(c(22) - y(2)))/(exp((c(22) - y(2))/c(33)) - 1);
J(15,2) = - (4*exp((4*y(2))/103 + 160/103)*(y(15) - 1/(exp(-(c(26) - y(2))/c(29)) + 1)))/103 - (exp(-(c(26) - y(2))/c(29))*exp((4*y(2))/103 + 160/103))/(c(29)*(exp(-(c(26) - y(2))/c(29)) + 1)^2);
J(15,15) = -exp((4*y(2))/103 + 160/103);
J(16,2) = (c(16)*(y(16) - 1))/(exp((c(21) - y(2))/c(32)) - 1) + (c(19)*y(16)*exp((c(21) - y(2))/c(34)))/c(34) - (c(16)*exp((c(21) - y(2))/c(32))*(c(21) - y(2))*(y(16) - 1))/(c(32)*(exp((c(21) - y(2))/c(32)) - 1)^2);
J(16,16) = - c(19)*exp((c(21) - y(2))/c(34)) - (c(16)*(c(21) - y(2)))/(exp((c(21) - y(2))/c(32)) - 1);
J(17,2) = (c(15)*exp((c(23) - y(2))/c(30))*(y(17) - 1))/c(30) - (c(18)*y(17)*exp((c(23) - y(2))/c(31)))/(c(31)*(exp((c(23) - y(2))/c(31)) + 1)^2);
J(17,17) = - c(18)/(exp((c(23) - y(2))/c(31)) + 1) - c(15)*exp((c(23) - y(2))/c(30));
J(18,2) = - (y(18) - 1/(exp(-(c(25) - y(2))/c(28)) + 1))*((113*c(49))/857100000 + (113*y(2))/857100000) - (exp(-(c(25) - y(2))/c(28))*((113*(c(49)/100 + y(2)/100)^2)/171420 + 1/42855))/(c(28)*(exp(-(c(25) - y(2))/c(28)) + 1)^2);
J(18,18) = - (113*(c(49)/100 + y(2)/100)^2)/171420 - 1/42855;
J(19,2) = - (c(57)*c(58)*y(20)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(58)*y(19)*exp(-(c(60) - y(2))/(8*c(59))))/(4*c(57)*c(59));
J(19,19) = - c(56) - (2*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/c(57);
J(19,20) = (c(57)*c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(19,24) = c(55);
J(20,2) = (c(57)*c(58)*y(20)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(57)*c(58)*y(21)*exp((c(60) - y(2))/(8*c(59))))/(8*c(59)) + (c(58)*y(19)*exp(-(c(60) - y(2))/(8*c(59))))/(4*c(57)*c(59)) - (3*c(58)*y(20)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(57)*c(59));
J(20,19) = (2*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/c(57);
J(20,20) = - c(56)*c(57) - (c(57)*c(58)*exp((c(60) - y(2))/(8*c(59))))/2 - (3*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/(2*c(57));
J(20,21) = c(57)*c(58)*exp((c(60) - y(2))/(8*c(59)));
J(20,25) = c(55)/c(57);
J(21,2) = (c(57)*c(58)*y(21)*exp((c(60) - y(2))/(8*c(59))))/(8*c(59)) - (3*c(57)*c(58)*y(22)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) + (3*c(58)*y(20)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(57)*c(59)) - (c(58)*y(21)*exp(-(c(60) - y(2))/(8*c(59))))/(8*c(57)*c(59));
J(21,20) = (3*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/(2*c(57));
J(21,21) = - c(56)*c(57)^2 - c(57)*c(58)*exp((c(60) - y(2))/(8*c(59))) - (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/c(57);
J(21,22) = (3*c(57)*c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(21,26) = c(55)/c(57)^2;
J(22,2) = (3*c(57)*c(58)*y(22)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(57)*c(58)*y(23)*exp((c(60) - y(2))/(8*c(59))))/(4*c(59)) + (c(58)*y(21)*exp(-(c(60) - y(2))/(8*c(59))))/(8*c(57)*c(59)) - (c(58)*y(22)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(57)*c(59));
J(22,21) = (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/c(57);
J(22,22) = - c(56)*c(57)^3 - (3*c(57)*c(58)*exp((c(60) - y(2))/(8*c(59))))/2 - (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/(2*c(57));
J(22,23) = 2*c(57)*c(58)*exp((c(60) - y(2))/(8*c(59)));
J(22,27) = c(55)/c(57)^3;
J(23,2) = (c(57)*c(58)*y(23)*exp((c(60) - y(2))/(8*c(59))))/(4*c(59)) + (c(58)*y(22)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(57)*c(59));
J(23,22) = (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/(2*c(57));
J(23,23) = - c(56)*c(57)^4 - 2*c(57)*c(58)*exp((c(60) - y(2))/(8*c(59)));
J(23,28) = c(55)/c(57)^4;
J(24,2) = - (c(58)*y(24)*exp(-(c(60) - y(2))/(8*c(59))))/(4*c(59)) - (c(58)*y(25)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59));
J(24,19) = c(56);
J(24,24) = - c(55) - 2*c(58)*exp(-(c(60) - y(2))/(8*c(59)));
J(24,25) = (c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(25,2) = (c(58)*y(24)*exp(-(c(60) - y(2))/(8*c(59))))/(4*c(59)) - (3*c(58)*y(25)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(59)) + (c(58)*y(25)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(58)*y(26)*exp((c(60) - y(2))/(8*c(59))))/(8*c(59));
J(25,20) = c(56)*c(57);
J(25,24) = 2*c(58)*exp(-(c(60) - y(2))/(8*c(59)));
J(25,25) = - c(55)/c(57) - (3*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/2 - (c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(25,26) = c(58)*exp((c(60) - y(2))/(8*c(59)));
J(26,2) = (3*c(58)*y(25)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(58)*y(26)*exp(-(c(60) - y(2))/(8*c(59))))/(8*c(59)) + (c(58)*y(26)*exp((c(60) - y(2))/(8*c(59))))/(8*c(59)) - (3*c(58)*y(27)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59));
J(26,21) = c(56)*c(57)^2;
J(26,25) = (3*c(58)*exp(-(c(60) - y(2))/(8*c(59))))/2;
J(26,26) = - c(55)/c(57)^2 - c(58)*exp(-(c(60) - y(2))/(8*c(59))) - c(58)*exp((c(60) - y(2))/(8*c(59)));
J(26,27) = (3*c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(27,2) = (c(58)*y(26)*exp(-(c(60) - y(2))/(8*c(59))))/(8*c(59)) - (c(58)*y(27)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(59)) + (3*c(58)*y(27)*exp((c(60) - y(2))/(8*c(59))))/(16*c(59)) - (c(58)*y(28)*exp((c(60) - y(2))/(8*c(59))))/(4*c(59));
J(27,22) = c(56)*c(57)^3;
J(27,26) = c(58)*exp(-(c(60) - y(2))/(8*c(59)));
J(27,27) = - c(55)/c(57)^3 - (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/2 - (3*c(58)*exp((c(60) - y(2))/(8*c(59))))/2;
J(27,28) = 2*c(58)*exp((c(60) - y(2))/(8*c(59)));
J(28,2) = (c(58)*y(27)*exp(-(c(60) - y(2))/(8*c(59))))/(16*c(59)) + (c(58)*y(28)*exp((c(60) - y(2))/(8*c(59))))/(4*c(59));
J(28,23) = c(56)*c(57)^4;
J(28,27) = (c(58)*exp(-(c(60) - y(2))/(8*c(59))))/2;
J(28,28) = - c(55)/c(57)^4 - 2*c(58)*exp((c(60) - y(2))/(8*c(59)));
J(29,19) = -(c(99)*(y(29) - y(30)))/c(102);
J(29,20) = -(c(99)*(y(29) - y(30)))/c(102);
J(29,21) = -(c(99)*(y(29) - y(30)))/c(102);
J(29,22) = -(c(99)*(y(29) - y(30)))/c(102);
J(29,23) = -(c(99)*(y(29) - y(30)))/c(102);
J(29,29) = c(71)*(y(34) - c(73) + y(36)) - c(63)/c(102) - c(64)/c(102) - c(79)*y(44) - c(61)/(c(102)*(c(62) + y(29))) - (c(99)*(y(19) + y(20) + y(21) + y(22) + y(23)))/c(102) + (c(61)*y(29))/(c(102)*(c(62) + y(29))^2);
J(29,30) = c(63)/c(102) + (c(99)*(y(19) + y(20) + y(21) + y(22) + y(23)))/c(102);
J(29,31) = c(64)/c(102);
J(29,34) = c(72) + c(71)*y(29);
J(29,36) = c(71)*y(29);
J(29,40) = c(80);
J(29,44) = -c(79)*y(29);
J(30,19) = (c(99)*(y(29) - y(30)))/c(104);
J(30,20) = (c(99)*(y(29) - y(30)))/c(104);
J(30,21) = (c(99)*(y(29) - y(30)))/c(104);
J(30,22) = (c(99)*(y(29) - y(30)))/c(104);
J(30,23) = (c(99)*(y(29) - y(30)))/c(104);
J(30,29) = c(63)/c(104) + c(61)/(c(104)*(c(62) + y(29))) + (c(99)*(y(19) + y(20) + y(21) + y(22) + y(23)))/c(104) - (c(61)*y(29))/(c(104)*(c(62) + y(29))^2);
J(30,30) = - c(76)*(c(78) - y(38)) - c(63)/c(104) - c(65)/c(104) - (c(99)*(y(19) + y(20) + y(21) + y(22) + y(23)))/c(104);
J(30,32) = c(65)/c(104);
J(30,38) = c(77) + c(76)*y(30);
J(31,29) = c(64)/c(103);
J(31,31) = c(68)*(y(33) - c(70) + y(48) + y(49) + y(50) + y(51) + y(52) + y(53)) - c(68)*y(33) - c(68)*y(49) - c(68)*y(50) - c(79)*y(45) - c(63)/c(103) - c(64)/c(103) + c(71)*(y(35) - c(73) + y(37)) - c(61)/(c(103)*(c(62) + y(31))) + (c(61)*y(31))/(c(103)*(c(62) + y(31))^2);
J(31,32) = c(63)/c(103);
J(31,33) = c(69);
J(31,35) = c(72) + c(71)*y(31);
J(31,37) = c(71)*y(31);
J(31,41) = c(80);
J(31,45) = -c(79)*y(31);
J(31,48) = c(69) + c(68)*y(31);
J(31,50) = c(69);
J(31,51) = c(69) + c(68)*y(31);
J(31,52) = c(68)*y(31);
J(31,53) = c(68)*y(31);
J(32,30) = c(65)/c(105);
J(32,31) = c(63)/c(105) + c(61)/(c(105)*(c(62) + y(31))) - (c(61)*y(31))/(c(105)*(c(62) + y(31))^2);
J(32,32) = c(96)*h73*y(55)*(c(98) - (y(32)*y(55))/1000) - c(63)/c(105) - c(65)/c(105) - c(97)*h74*y(55)*y(56) - (c(96)*h73*y(32)*y(55)^2)/1000 - c(76)*(c(78) - y(39));
J(32,39) = c(77) + c(76)*y(32);
J(32,55) = c(96)*h73*y(32)*(c(98) - (y(32)*y(55))/1000) - (c(96)*h73*y(32)^2*y(55))/1000 - c(97)*h74*y(32)*y(56);
J(32,56) = 1000*c(97)*h74*(c(98) - (y(32)*y(55))/1000);
J(33,31) = - c(68)*(y(33) - c(70) + y(48) + y(49) + y(50) + y(51) + y(52) + y(53)) - c(68)*y(33);
J(33,33) = - c(69) - c(85) - 2*c(68)*y(31);
J(33,48) = c(69) - c(68)*y(31);
J(33,49) = -c(68)*y(31);
J(33,50) = c(86) - c(68)*y(31);
J(33,51) = -c(68)*y(31);
J(33,52) = -c(68)*y(31);
J(33,53) = -c(68)*y(31);
J(34,29) = -c(71)*(y(34) - c(73) + y(36));
J(34,34) = - c(72) - c(71)*y(29);
J(34,36) = -c(71)*y(29);
J(35,31) = -c(71)*(y(35) - c(73) + y(37));
J(35,35) = - c(72) - c(71)*y(31);
J(35,37) = -c(71)*y(31);
J(36,34) = -c(74)*y(46);
J(36,36) = - c(75) - c(74)*y(46);
J(36,46) = -c(74)*(y(34) - c(73) + y(36));
J(37,35) = -c(74)*y(47);
J(37,37) = - c(75) - c(74)*y(47);
J(37,47) = -c(74)*(y(35) - c(73) + y(37));
J(38,30) = c(76)*(c(78) - y(38));
J(38,38) = - c(77) - c(76)*y(30);
J(39,32) = c(76)*(c(78) - y(39));
J(39,39) = - c(77) - c(76)*y(32);
J(40,29) = c(79)*y(44);
J(40,40) = - c(80) - c(83)/c(102);
J(40,41) = c(83)/c(102);
J(40,44) = c(79)*y(29);
J(41,31) = c(79)*y(45);
J(41,40) = c(83)/c(103);
J(41,41) = - c(80) - c(83)/c(103);
J(41,45) = c(79)*y(31);
J(42,42) = - c(82) - c(83)/c(102);
J(42,43) = c(83)/c(102);
J(42,44) = c(81)*y(46);
J(42,46) = c(81)*y(44);
J(43,42) = c(83)/c(103);
J(43,43) = - c(82) - c(83)/c(103);
J(43,45) = c(81)*y(47);
J(43,47) = c(81)*y(45);
J(44,29) = -c(79)*y(44);
J(44,40) = c(80);
J(44,42) = c(82);
J(44,44) = - c(79)*y(29) - c(81)*y(46) - c(83)/c(102);
J(44,45) = c(83)/c(102);
J(44,46) = -c(81)*y(44);
J(45,31) = -c(79)*y(45);
J(45,41) = c(80);
J(45,43) = c(82);
J(45,44) = c(83)/c(103);
J(45,45) = - c(79)*y(31) - c(81)*y(47) - c(83)/c(103);
J(45,47) = -c(81)*y(45);
J(46,34) = c(74)*y(46);
J(46,36) = c(75) + c(74)*y(46);
J(46,42) = c(82);
J(46,44) = -c(81)*y(46);
J(46,46) = c(74)*(y(34) - c(73) + y(36)) - c(84)/c(102) - c(81)*y(44);
J(46,47) = c(84)/c(102);
J(47,35) = c(74)*y(47);
J(47,37) = c(75) + c(74)*y(47);
J(47,43) = c(82);
J(47,45) = -c(81)*y(47);
J(47,46) = c(84)/c(103);
J(47,47) = c(74)*(y(35) - c(73) + y(37)) - c(84)/c(103) - c(81)*y(45);
J(48,31) = c(68)*y(33);
J(48,33) = c(68)*y(31);
J(48,48) = - c(69) - c(87);
J(48,51) = c(88);
J(49,31) = -c(68)*y(49);
J(49,33) = -c(85);
J(49,48) = -c(85);
J(49,49) = - c(85) - c(86) - c(68)*y(31);
J(49,50) = c(69) - c(85);
J(49,51) = -c(85);
J(49,52) = -c(85);
J(49,53) = -c(85);
J(50,31) = c(68)*y(49) - c(68)*y(50);
J(50,33) = c(85);
J(50,49) = c(68)*y(31);
J(50,50) = - c(69) - c(86) - c(68)*y(31);
J(50,51) = c(69);
J(51,31) = c(68)*y(50);
J(51,48) = c(87);
J(51,50) = c(68)*y(31);
J(51,51) = - c(69) - c(88) - c(89);
J(51,52) = c(90);
J(51,53) = c(93);
J(52,51) = c(89);
J(52,52) = - c(90) - c(91);
J(52,53) = c(92);
J(53,52) = c(91);
J(53,53) = - c(92) - c(93);
J(54,52) = c(91)/1000;
J(54,53) = -c(92)/1000;
J(54,54) = - c(94) - c(95)/c(103);
J(54,55) = c(95)/c(103);
J(55,32) = (c(96)*h73*y(55)*(c(98) - (y(32)*y(55))/1000))/1000 - (c(96)*h73*y(32)*y(55)^2)/1000000 - (c(97)*h74*y(55)*y(56))/1000;
J(55,54) = c(95)/c(105);
J(55,55) = (c(96)*h73*y(32)*(c(98) - (y(32)*y(55))/1000))/1000 - (c(97)*h74*y(32)*y(56))/1000 - (c(96)*h73*y(32)^2*y(55))/1000000 - c(95)/c(105);
J(55,56) = c(97)*h74*(c(98) - (y(32)*y(55))/1000);
J(56,32) = (c(97)*h74*y(55)*y(56))/1000 + (c(96)*h73*y(32)*y(55)^2)/1000000 - (c(96)*h73*y(55)*(c(98) - (y(32)*y(55))/1000))/1000;
J(56,55) = (c(97)*h74*y(32)*y(56))/1000 + (c(96)*h73*y(32)^2*y(55))/1000000 - (c(96)*h73*y(32)*(c(98) - (y(32)*y(55))/1000))/1000;
J(56,56) = -c(97)*h74*(c(98) - (y(32)*y(55))/1000);
end
}
};
};