function out = JR_slow
out{1} = @init;
out{2} = @fun_eval;
out{3} = @jacobian;
out{4} = @jacobianp;
out{5} = @hessians;
out{6} = @hessiansp;
out{7} = @der3;
out{8} = [];
out{9} = [];

% --------------------------------------------------------------------------
function dydt = fun_eval(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
C1=par_CC;
C2=0.8*par_CC;
C3=0.25*par_CC;
C4=0.25*par_CC;
S4=2*par_e0/(1+exp(par_r*(par_v0-(kmrgd(3)-par_theta*kmrgd(5)-(1-par_theta)*kmrgd(7)))));
S5=2*par_e0/(1+exp(par_r*(par_v0-(C1*kmrgd(1)))));
S6=2*par_e0/(1+exp(par_r*(par_v0-(C3*kmrgd(1)))));
S7=S6;
dydt=[kmrgd(2);
par_A*par_a*S4-2*par_a*kmrgd(2)-par_a*par_a*kmrgd(1);
kmrgd(4);
par_A*par_a*(par_II+C2*S5)-2*par_a*kmrgd(4)-par_a*par_a*kmrgd(3);
kmrgd(6);
par_Bf*par_bf*C4*S6-2*par_bf*kmrgd(6)-par_bf*par_bf*kmrgd(5);
kmrgd(8);
par_Bs*par_bs*C4*S7-2*par_bs*kmrgd(8)-par_bs*par_bs*kmrgd(7);];

% --------------------------------------------------------------------------
function [tspan,y0,options] = init
handles = feval(JR_slow);
y0=[0,0,0,0,0,0,0,0];
options = odeset('Jacobian',handles(3),'JacobianP',handles(4),'Hessians',handles(5),'HessiansP',handles(6));
tspan = [0 10];

% --------------------------------------------------------------------------
function jac = jacobian(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
jac=[ 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 ; -par_a^2 , -2*par_a , (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , -(2*par_A*par_a*par_e0*par_r*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 ; (8*par_A*par_a*par_e0*par_r*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , -par_a^2 , -2*par_a , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 ; (par_e0*par_r*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , -par_bf^2 , -2*par_bf , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 ; (par_e0*par_r*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , -par_bs^2 , -2*par_bs ];
% --------------------------------------------------------------------------
function jacp = jacobianp(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
jacp=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , (2*par_a*par_e0)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1) , (2*par_A*par_e0)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1) - 2*kmrgd(1)*par_a - 2*kmrgd(2) , 0 , 0 , 0 , 0 , (2*par_A*par_a)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1) , (2*par_A*par_a*par_e0*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , -(2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , -(2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; par_A*par_a , par_a*(par_II + (8*par_e0*par_CC)/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1))) , par_A*(par_II + (8*par_e0*par_CC)/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1))) - 2*kmrgd(3)*par_a - 2*kmrgd(4) , 0 , 0 , 0 , 0 , (8*par_A*par_a*par_CC)/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)) , -(8*par_A*par_a*par_e0*par_CC*exp(par_r*(par_v0 - kmrgd(1)*par_CC))*(par_v0 - kmrgd(1)*par_CC))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , -(8*par_A*par_a*par_e0*par_r*par_CC*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , par_A*par_a*((8*par_e0)/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)) + (8*kmrgd(1)*par_e0*par_r*par_CC*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2)) ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , (par_e0*par_CC*par_bf)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) , (par_e0*par_CC*par_Bf)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) - 2*kmrgd(5)*par_bf - 2*kmrgd(6) , (par_CC*par_Bf*par_bf)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) , -(par_e0*par_CC*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , -(par_e0*par_r*par_CC*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , (par_e0*par_Bf*par_bf)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) + (kmrgd(1)*par_e0*par_r*par_CC*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , (par_e0*par_CC*par_bs)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) , (par_e0*par_CC*par_Bs)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) - 2*kmrgd(7)*par_bs - 2*kmrgd(8) , 0 , 0 , (par_CC*par_Bs*par_bs)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) , -(par_e0*par_CC*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , -(par_e0*par_r*par_CC*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , (par_e0*par_Bs*par_bs)/(2*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)) + (kmrgd(1)*par_e0*par_r*par_CC*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) ];
% --------------------------------------------------------------------------
function hess = hessians(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
hess1=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (16*par_A*par_a*par_e0*par_r^2*par_CC^3*exp(2*par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^3) - (8*par_A*par_a*par_e0*par_r^2*par_CC^3*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^2*par_CC^3*par_Bf*par_bf*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(16*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) - (par_e0*par_r^2*par_CC^3*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(32*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^2*par_CC^3*par_Bs*par_bs*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(16*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) - (par_e0*par_r^2*par_CC^3*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(32*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess2=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess3=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess4=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess5=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (4*par_A*par_a*par_e0*par_r^2*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess6=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess7=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess8=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hess(:,:,1) =hess1;
hess(:,:,2) =hess2;
hess(:,:,3) =hess3;
hess(:,:,4) =hess4;
hess(:,:,5) =hess5;
hess(:,:,6) =hess6;
hess(:,:,7) =hess7;
hess(:,:,8) =hess8;
% --------------------------------------------------------------------------
function hessp = hessiansp(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
hessp1=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp2=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , -(2*par_a*par_e0*par_r*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_a*par_e0*par_r*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp3=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; -2*par_a , -2 , (2*par_A*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , -(2*par_A*par_e0*par_r*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_A*par_e0*par_r*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , -2*par_a , -2 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp4=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC^2*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp5=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC^2*par_Bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , -2*par_bs , -2 ];
hessp6=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC^2*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp7=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC^2*par_Bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , -2*par_bf , -2 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp8=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , -(2*par_A*par_a*par_r*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_A*par_a*par_r*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_r*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_r*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp9=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 + (4*par_A*par_a*par_e0*par_r*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (2*par_A*par_a*par_e0*par_r*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (2*par_A*par_a*par_e0*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (2*par_A*par_a*par_e0*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 + (4*par_A*par_a*par_e0*par_r*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_A*par_a*par_e0*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) + (8*par_A*par_a*par_e0*par_r*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC))*(par_v0 - kmrgd(1)*par_CC))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) - (16*par_A*par_a*par_e0*par_r*par_CC^2*exp(2*par_r*(par_v0 - kmrgd(1)*par_CC))*(par_v0 - kmrgd(1)*par_CC))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) + (par_e0*par_r*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (par_e0*par_r*par_CC^2*par_Bf*par_bf*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) + (par_e0*par_r*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (par_e0*par_r*par_CC^2*par_Bs*par_bs*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4))*(par_v0 - (kmrgd(1)*par_CC)/4))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp10=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_A*par_a*par_e0*par_r^2*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) - (16*par_A*par_a*par_e0*par_r^2*par_CC^2*exp(2*par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^2*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (par_e0*par_r^2*par_CC^2*par_Bf*par_bf*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^2*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(8*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (par_e0*par_r^2*par_CC^2*par_Bs*par_bs*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp11=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 , (4*par_A*par_a*par_e0*par_r^2*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^2*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7)))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 , 0 , (2*par_A*par_a*par_e0*par_r*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 + (2*par_A*par_a*par_e0*par_r^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (4*par_A*par_a*par_e0*par_r^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(kmrgd(5) - kmrgd(7))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp12=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (16*par_A*par_a*par_e0*par_r*par_CC*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) - (8*kmrgd(1)*par_A*par_a*par_e0*par_r^2*par_CC^2*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) + (16*kmrgd(1)*par_A*par_a*par_e0*par_r^2*par_CC^2*exp(2*par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (kmrgd(1)*par_e0*par_r^2*par_CC^2*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(32*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) + (kmrgd(1)*par_e0*par_r^2*par_CC^2*par_Bf*par_bf*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(16*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r*par_CC*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(4*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (kmrgd(1)*par_e0*par_r^2*par_CC^2*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(32*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) + (kmrgd(1)*par_e0*par_r^2*par_CC^2*par_Bs*par_bs*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(16*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
hessp(:,:,1) =hessp1;
hessp(:,:,2) =hessp2;
hessp(:,:,3) =hessp3;
hessp(:,:,4) =hessp4;
hessp(:,:,5) =hessp5;
hessp(:,:,6) =hessp6;
hessp(:,:,7) =hessp7;
hessp(:,:,8) =hessp8;
hessp(:,:,9) =hessp9;
hessp(:,:,10) =hessp10;
hessp(:,:,11) =hessp11;
hessp(:,:,12) =hessp12;
%---------------------------------------------------------------------------
function tens3  = der3(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
tens31=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (8*par_A*par_a*par_e0*par_r^3*par_CC^4*exp(par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^2) - (48*par_A*par_a*par_e0*par_r^3*par_CC^4*exp(2*par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^3) + (48*par_A*par_a*par_e0*par_r^3*par_CC^4*exp(3*par_r*(par_v0 - kmrgd(1)*par_CC)))/(5*(exp(par_r*(par_v0 - kmrgd(1)*par_CC)) + 1)^4) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^3*par_CC^4*par_Bf*par_bf*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(128*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (3*par_e0*par_r^3*par_CC^4*par_Bf*par_bf*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(64*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) + (3*par_e0*par_r^3*par_CC^4*par_Bf*par_bf*exp(3*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(64*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^4) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; (par_e0*par_r^3*par_CC^4*par_Bs*par_bs*exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(128*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^2) - (3*par_e0*par_r^3*par_CC^4*par_Bs*par_bs*exp(2*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(64*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^3) + (3*par_e0*par_r^3*par_CC^4*par_Bs*par_bs*exp(3*par_r*(par_v0 - (kmrgd(1)*par_CC)/4)))/(64*(exp(par_r*(par_v0 - (kmrgd(1)*par_CC)/4)) + 1)^4) , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens32=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens33=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens34=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens35=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens36=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens37=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens38=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens39=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens310=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens311=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens312=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens313=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens314=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens315=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens316=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens317=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens318=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens319=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens320=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens321=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens322=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens323=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens324=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens325=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens326=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens327=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens328=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens329=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens330=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens331=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens332=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens333=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens334=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens335=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens336=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens337=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens338=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens339=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens340=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens341=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens342=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens343=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens344=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens345=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens346=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens347=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens348=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens349=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens350=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens351=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens352=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens353=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*par_theta^2*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1))/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens354=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens355=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 - (2*par_A*par_a*par_e0*par_r^3*par_theta*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*par_theta*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^2)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 , (2*par_A*par_a*par_e0*par_r^3*exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^3)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^2 - (12*par_A*par_a*par_e0*par_r^3*exp(-2*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^3)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^3 + (12*par_A*par_a*par_e0*par_r^3*exp(-3*par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1)))*(par_theta - 1)^3)/(exp(-par_r*(kmrgd(3) - par_v0 - kmrgd(5)*par_theta + kmrgd(7)*(par_theta - 1))) + 1)^4 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens356=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens357=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens358=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens359=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens360=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens361=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens362=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens363=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens364=[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ];
tens3(:,:,1,1) =tens31;
tens3(:,:,1,2) =tens32;
tens3(:,:,1,3) =tens33;
tens3(:,:,1,4) =tens34;
tens3(:,:,1,5) =tens35;
tens3(:,:,1,6) =tens36;
tens3(:,:,1,7) =tens37;
tens3(:,:,1,8) =tens38;
tens3(:,:,2,1) =tens39;
tens3(:,:,2,2) =tens310;
tens3(:,:,2,3) =tens311;
tens3(:,:,2,4) =tens312;
tens3(:,:,2,5) =tens313;
tens3(:,:,2,6) =tens314;
tens3(:,:,2,7) =tens315;
tens3(:,:,2,8) =tens316;
tens3(:,:,3,1) =tens317;
tens3(:,:,3,2) =tens318;
tens3(:,:,3,3) =tens319;
tens3(:,:,3,4) =tens320;
tens3(:,:,3,5) =tens321;
tens3(:,:,3,6) =tens322;
tens3(:,:,3,7) =tens323;
tens3(:,:,3,8) =tens324;
tens3(:,:,4,1) =tens325;
tens3(:,:,4,2) =tens326;
tens3(:,:,4,3) =tens327;
tens3(:,:,4,4) =tens328;
tens3(:,:,4,5) =tens329;
tens3(:,:,4,6) =tens330;
tens3(:,:,4,7) =tens331;
tens3(:,:,4,8) =tens332;
tens3(:,:,5,1) =tens333;
tens3(:,:,5,2) =tens334;
tens3(:,:,5,3) =tens335;
tens3(:,:,5,4) =tens336;
tens3(:,:,5,5) =tens337;
tens3(:,:,5,6) =tens338;
tens3(:,:,5,7) =tens339;
tens3(:,:,5,8) =tens340;
tens3(:,:,6,1) =tens341;
tens3(:,:,6,2) =tens342;
tens3(:,:,6,3) =tens343;
tens3(:,:,6,4) =tens344;
tens3(:,:,6,5) =tens345;
tens3(:,:,6,6) =tens346;
tens3(:,:,6,7) =tens347;
tens3(:,:,6,8) =tens348;
tens3(:,:,7,1) =tens349;
tens3(:,:,7,2) =tens350;
tens3(:,:,7,3) =tens351;
tens3(:,:,7,4) =tens352;
tens3(:,:,7,5) =tens353;
tens3(:,:,7,6) =tens354;
tens3(:,:,7,7) =tens355;
tens3(:,:,7,8) =tens356;
tens3(:,:,8,1) =tens357;
tens3(:,:,8,2) =tens358;
tens3(:,:,8,3) =tens359;
tens3(:,:,8,4) =tens360;
tens3(:,:,8,5) =tens361;
tens3(:,:,8,6) =tens362;
tens3(:,:,8,7) =tens363;
tens3(:,:,8,8) =tens364;
%---------------------------------------------------------------------------
function tens4  = der4(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
%---------------------------------------------------------------------------
function tens5  = der5(t,kmrgd,par_II,par_A,par_a,par_Bs,par_bs,par_Bf,par_bf,par_e0,par_r,par_v0,par_theta,par_CC)
