/* * mg_2ode_avf.c * * AUTO C file for the vector field named: mg_2ode * * This file was generated by the program VFGEN (Version:2.4.0) * Generated on 10-Jul-2008 at 11:44 */ #include #include "auto_f2c.h" /* * FUNC Defines the vector field and its derivatives */ int func(integer ndim_, const doublereal *u_, const integer *icp_, const doublereal *par_, integer ijac_, doublereal *f_, doublereal *dfdu_, doublereal *dfdp_) { integer dfdu__dim1, dfdp__dim1; double x, symbol9_1_1, symbol9_2_1, symbol9_1_2, symbol9_2_2, symbol9_1_3, symbol9_2_3, symbol9_1_4, symbol9_2_4, symbol9_1_5, symbol9_2_5, symbol9_1_6, symbol9_2_6, symbol9_1_7, symbol9_2_7, symbol9_1_8, symbol9_2_8, symbol9_1_9, symbol9_2_9, symbol9_1_10, symbol9_2_10, symbol9_1_11, symbol9_2_11, symbol9_1_12, symbol9_2_12, symbol9_1_13, symbol9_2_13, symbol9_1_14, symbol9_2_14, symbol9_1_15, symbol9_2_15, symbol9_1_16, symbol9_2_16, symbol9_1_17, symbol9_2_17, symbol9_1_18, symbol9_2_18, symbol9_1_19, symbol9_2_19, symbol9_1_20, symbol9_2_20; double a, b, tau; double delayedx; dfdu__dim1 = ndim_; dfdp__dim1 = ndim_; x = u_[0]; symbol9_1_1 = u_[1]; symbol9_2_1 = u_[2]; symbol9_1_2 = u_[3]; symbol9_2_2 = u_[4]; symbol9_1_3 = u_[5]; symbol9_2_3 = u_[6]; symbol9_1_4 = u_[7]; symbol9_2_4 = u_[8]; symbol9_1_5 = u_[9]; symbol9_2_5 = u_[10]; symbol9_1_6 = u_[11]; symbol9_2_6 = u_[12]; symbol9_1_7 = u_[13]; symbol9_2_7 = u_[14]; symbol9_1_8 = u_[15]; symbol9_2_8 = u_[16]; symbol9_1_9 = u_[17]; symbol9_2_9 = u_[18]; symbol9_1_10 = u_[19]; symbol9_2_10 = u_[20]; symbol9_1_11 = u_[21]; symbol9_2_11 = u_[22]; symbol9_1_12 = u_[23]; symbol9_2_12 = u_[24]; symbol9_1_13 = u_[25]; symbol9_2_13 = u_[26]; symbol9_1_14 = u_[27]; symbol9_2_14 = u_[28]; symbol9_1_15 = u_[29]; symbol9_2_15 = u_[30]; symbol9_1_16 = u_[31]; symbol9_2_16 = u_[32]; symbol9_1_17 = u_[33]; symbol9_2_17 = u_[34]; symbol9_1_18 = u_[35]; symbol9_2_18 = u_[36]; symbol9_1_19 = u_[37]; symbol9_2_19 = u_[38]; symbol9_1_20 = u_[39]; symbol9_2_20 = u_[40]; a = par_[0]; b = par_[1]; tau = par_[2]; delayedx = symbol9_1_20; /* * The vector field */ f_[0] = delayedx/( ((delayedx*(delayedx*delayedx)*(delayedx*delayedx))*(delayedx*(delayedx*delayedx)*(delayedx*delayedx)))+1.0)*a-b*x; f_[1] = symbol9_2_1; f_[2] = -40.0*1.0/tau*( symbol9_2_1+20.0*( symbol9_1_1-x)/tau); f_[3] = symbol9_2_2; f_[4] = -40.0*( symbol9_2_2+-20.0*1.0/tau*( symbol9_1_1-symbol9_1_2))/tau; f_[5] = symbol9_2_3; f_[6] = -40.0*1.0/tau*( 20.0*( symbol9_1_3-symbol9_1_2)/tau+symbol9_2_3); f_[7] = symbol9_2_4; f_[8] = -40.0*1.0/tau*( symbol9_2_4+-20.0*1.0/tau*( symbol9_1_3-symbol9_1_4)); f_[9] = symbol9_2_5; f_[10] = -40.0*( 20.0*1.0/tau*( symbol9_1_5-symbol9_1_4)+symbol9_2_5)/tau; f_[11] = symbol9_2_6; f_[12] = 40.0*1.0/tau*( 20.0*( symbol9_1_5-symbol9_1_6)/tau-symbol9_2_6); f_[13] = symbol9_2_7; f_[14] = 40.0*1.0/tau*( 20.0*1.0/tau*( symbol9_1_6-symbol9_1_7)-symbol9_2_7); f_[15] = symbol9_2_8; f_[16] = -40.0*1.0/tau*( symbol9_2_8+20.0*1.0/tau*( symbol9_1_8-symbol9_1_7)); f_[17] = symbol9_2_9; f_[18] = -40.0*1.0/tau*( symbol9_2_9+-20.0*( symbol9_1_8-symbol9_1_9)/tau); f_[19] = symbol9_2_10; f_[20] = 40.0*( 20.0*( symbol9_1_9-symbol9_1_10)/tau-symbol9_2_10)/tau; f_[21] = symbol9_2_11; f_[22] = -40.0*( 20.0*( symbol9_1_11-symbol9_1_10)/tau+symbol9_2_11)/tau; f_[23] = symbol9_2_12; f_[24] = 40.0*( 20.0*1.0/tau*( symbol9_1_11-symbol9_1_12)-symbol9_2_12)/tau; f_[25] = symbol9_2_13; f_[26] = -40.0*1.0/tau*( 20.0*1.0/tau*( symbol9_1_13-symbol9_1_12)+symbol9_2_13); f_[27] = symbol9_2_14; f_[28] = -40.0*( symbol9_2_14+20.0*1.0/tau*( symbol9_1_14-symbol9_1_13))/tau; f_[29] = symbol9_2_15; f_[30] = -40.0*1.0/tau*( symbol9_2_15+-20.0*( symbol9_1_14-symbol9_1_15)/tau); f_[31] = symbol9_2_16; f_[32] = -40.0*( symbol9_2_16+-20.0*1.0/tau*( symbol9_1_15-symbol9_1_16))/tau; f_[33] = symbol9_2_17; f_[34] = -40.0*( 20.0*1.0/tau*( symbol9_1_17-symbol9_1_16)+symbol9_2_17)/tau; f_[35] = symbol9_2_18; f_[36] = -40.0*( symbol9_2_18+-20.0*1.0/tau*( symbol9_1_17-symbol9_1_18))/tau; f_[37] = symbol9_2_19; f_[38] = -40.0*1.0/tau*( 20.0*( symbol9_1_19-symbol9_1_18)/tau+symbol9_2_19); f_[39] = symbol9_2_20; f_[40] = -40.0*1.0/tau*( symbol9_2_20+-20.0*1.0/tau*( symbol9_1_19-symbol9_1_20)); if (ijac_ == 0) return 0; /* * The Jacobian */ ARRAY2D(dfdu_,0,0) = -b; ARRAY2D(dfdu_,0,1) = 0.0; ARRAY2D(dfdu_,0,2) = 0.0; ARRAY2D(dfdu_,0,3) = 0.0; ARRAY2D(dfdu_,0,4) = 0.0; ARRAY2D(dfdu_,0,5) = 0.0; ARRAY2D(dfdu_,0,6) = 0.0; ARRAY2D(dfdu_,0,7) = 0.0; ARRAY2D(dfdu_,0,8) = 0.0; ARRAY2D(dfdu_,0,9) = 0.0; ARRAY2D(dfdu_,0,10) = 0.0; ARRAY2D(dfdu_,0,11) = 0.0; ARRAY2D(dfdu_,0,12) = 0.0; ARRAY2D(dfdu_,0,13) = 0.0; ARRAY2D(dfdu_,0,14) = 0.0; ARRAY2D(dfdu_,0,15) = 0.0; ARRAY2D(dfdu_,0,16) = 0.0; ARRAY2D(dfdu_,0,17) = 0.0; ARRAY2D(dfdu_,0,18) = 0.0; ARRAY2D(dfdu_,0,19) = 0.0; ARRAY2D(dfdu_,0,20) = 0.0; ARRAY2D(dfdu_,0,21) = 0.0; ARRAY2D(dfdu_,0,22) = 0.0; ARRAY2D(dfdu_,0,23) = 0.0; ARRAY2D(dfdu_,0,24) = 0.0; ARRAY2D(dfdu_,0,25) = 0.0; ARRAY2D(dfdu_,0,26) = 0.0; ARRAY2D(dfdu_,0,27) = 0.0; ARRAY2D(dfdu_,0,28) = 0.0; ARRAY2D(dfdu_,0,29) = 0.0; ARRAY2D(dfdu_,0,30) = 0.0; ARRAY2D(dfdu_,0,31) = 0.0; ARRAY2D(dfdu_,0,32) = 0.0; ARRAY2D(dfdu_,0,33) = 0.0; ARRAY2D(dfdu_,0,34) = 0.0; ARRAY2D(dfdu_,0,35) = 0.0; ARRAY2D(dfdu_,0,36) = 0.0; ARRAY2D(dfdu_,0,37) = 0.0; ARRAY2D(dfdu_,0,38) = 0.0; ARRAY2D(dfdu_,0,39) = -10.0*a*((symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20))*(symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20)))/pow( ((symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20))*(symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20)))+1.0,2.0)+a/( ((symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20))*(symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20)))+1.0); ARRAY2D(dfdu_,0,40) = 0.0; ARRAY2D(dfdu_,1,0) = 0.0; ARRAY2D(dfdu_,1,1) = 0.0; ARRAY2D(dfdu_,1,2) = 1.0; ARRAY2D(dfdu_,1,3) = 0.0; ARRAY2D(dfdu_,1,4) = 0.0; ARRAY2D(dfdu_,1,5) = 0.0; ARRAY2D(dfdu_,1,6) = 0.0; ARRAY2D(dfdu_,1,7) = 0.0; ARRAY2D(dfdu_,1,8) = 0.0; ARRAY2D(dfdu_,1,9) = 0.0; ARRAY2D(dfdu_,1,10) = 0.0; ARRAY2D(dfdu_,1,11) = 0.0; ARRAY2D(dfdu_,1,12) = 0.0; ARRAY2D(dfdu_,1,13) = 0.0; ARRAY2D(dfdu_,1,14) = 0.0; ARRAY2D(dfdu_,1,15) = 0.0; ARRAY2D(dfdu_,1,16) = 0.0; ARRAY2D(dfdu_,1,17) = 0.0; ARRAY2D(dfdu_,1,18) = 0.0; ARRAY2D(dfdu_,1,19) = 0.0; ARRAY2D(dfdu_,1,20) = 0.0; ARRAY2D(dfdu_,1,21) = 0.0; ARRAY2D(dfdu_,1,22) = 0.0; ARRAY2D(dfdu_,1,23) = 0.0; ARRAY2D(dfdu_,1,24) = 0.0; ARRAY2D(dfdu_,1,25) = 0.0; ARRAY2D(dfdu_,1,26) = 0.0; ARRAY2D(dfdu_,1,27) = 0.0; ARRAY2D(dfdu_,1,28) = 0.0; ARRAY2D(dfdu_,1,29) = 0.0; ARRAY2D(dfdu_,1,30) = 0.0; ARRAY2D(dfdu_,1,31) = 0.0; ARRAY2D(dfdu_,1,32) = 0.0; ARRAY2D(dfdu_,1,33) = 0.0; ARRAY2D(dfdu_,1,34) = 0.0; ARRAY2D(dfdu_,1,35) = 0.0; ARRAY2D(dfdu_,1,36) = 0.0; ARRAY2D(dfdu_,1,37) = 0.0; ARRAY2D(dfdu_,1,38) = 0.0; ARRAY2D(dfdu_,1,39) = 0.0; ARRAY2D(dfdu_,1,40) = 0.0; ARRAY2D(dfdu_,2,0) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,2,1) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,2,2) = -40.0*1.0/tau; ARRAY2D(dfdu_,2,3) = 0.0; ARRAY2D(dfdu_,2,4) = 0.0; ARRAY2D(dfdu_,2,5) = 0.0; ARRAY2D(dfdu_,2,6) = 0.0; ARRAY2D(dfdu_,2,7) = 0.0; ARRAY2D(dfdu_,2,8) = 0.0; ARRAY2D(dfdu_,2,9) = 0.0; ARRAY2D(dfdu_,2,10) = 0.0; ARRAY2D(dfdu_,2,11) = 0.0; ARRAY2D(dfdu_,2,12) = 0.0; ARRAY2D(dfdu_,2,13) = 0.0; ARRAY2D(dfdu_,2,14) = 0.0; ARRAY2D(dfdu_,2,15) = 0.0; ARRAY2D(dfdu_,2,16) = 0.0; ARRAY2D(dfdu_,2,17) = 0.0; ARRAY2D(dfdu_,2,18) = 0.0; ARRAY2D(dfdu_,2,19) = 0.0; ARRAY2D(dfdu_,2,20) = 0.0; ARRAY2D(dfdu_,2,21) = 0.0; ARRAY2D(dfdu_,2,22) = 0.0; ARRAY2D(dfdu_,2,23) = 0.0; ARRAY2D(dfdu_,2,24) = 0.0; ARRAY2D(dfdu_,2,25) = 0.0; ARRAY2D(dfdu_,2,26) = 0.0; ARRAY2D(dfdu_,2,27) = 0.0; ARRAY2D(dfdu_,2,28) = 0.0; ARRAY2D(dfdu_,2,29) = 0.0; ARRAY2D(dfdu_,2,30) = 0.0; ARRAY2D(dfdu_,2,31) = 0.0; ARRAY2D(dfdu_,2,32) = 0.0; ARRAY2D(dfdu_,2,33) = 0.0; ARRAY2D(dfdu_,2,34) = 0.0; ARRAY2D(dfdu_,2,35) = 0.0; ARRAY2D(dfdu_,2,36) = 0.0; ARRAY2D(dfdu_,2,37) = 0.0; ARRAY2D(dfdu_,2,38) = 0.0; ARRAY2D(dfdu_,2,39) = 0.0; ARRAY2D(dfdu_,2,40) = 0.0; ARRAY2D(dfdu_,3,0) = 0.0; ARRAY2D(dfdu_,3,1) = 0.0; ARRAY2D(dfdu_,3,2) = 0.0; ARRAY2D(dfdu_,3,3) = 0.0; ARRAY2D(dfdu_,3,4) = 1.0; ARRAY2D(dfdu_,3,5) = 0.0; ARRAY2D(dfdu_,3,6) = 0.0; ARRAY2D(dfdu_,3,7) = 0.0; ARRAY2D(dfdu_,3,8) = 0.0; ARRAY2D(dfdu_,3,9) = 0.0; ARRAY2D(dfdu_,3,10) = 0.0; ARRAY2D(dfdu_,3,11) = 0.0; ARRAY2D(dfdu_,3,12) = 0.0; ARRAY2D(dfdu_,3,13) = 0.0; ARRAY2D(dfdu_,3,14) = 0.0; ARRAY2D(dfdu_,3,15) = 0.0; ARRAY2D(dfdu_,3,16) = 0.0; ARRAY2D(dfdu_,3,17) = 0.0; ARRAY2D(dfdu_,3,18) = 0.0; ARRAY2D(dfdu_,3,19) = 0.0; ARRAY2D(dfdu_,3,20) = 0.0; ARRAY2D(dfdu_,3,21) = 0.0; ARRAY2D(dfdu_,3,22) = 0.0; ARRAY2D(dfdu_,3,23) = 0.0; ARRAY2D(dfdu_,3,24) = 0.0; ARRAY2D(dfdu_,3,25) = 0.0; ARRAY2D(dfdu_,3,26) = 0.0; ARRAY2D(dfdu_,3,27) = 0.0; ARRAY2D(dfdu_,3,28) = 0.0; ARRAY2D(dfdu_,3,29) = 0.0; ARRAY2D(dfdu_,3,30) = 0.0; ARRAY2D(dfdu_,3,31) = 0.0; ARRAY2D(dfdu_,3,32) = 0.0; ARRAY2D(dfdu_,3,33) = 0.0; ARRAY2D(dfdu_,3,34) = 0.0; ARRAY2D(dfdu_,3,35) = 0.0; ARRAY2D(dfdu_,3,36) = 0.0; ARRAY2D(dfdu_,3,37) = 0.0; ARRAY2D(dfdu_,3,38) = 0.0; ARRAY2D(dfdu_,3,39) = 0.0; ARRAY2D(dfdu_,3,40) = 0.0; ARRAY2D(dfdu_,4,0) = 0.0; ARRAY2D(dfdu_,4,1) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,4,2) = 0.0; ARRAY2D(dfdu_,4,3) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,4,4) = -40.0*1.0/tau; ARRAY2D(dfdu_,4,5) = 0.0; ARRAY2D(dfdu_,4,6) = 0.0; ARRAY2D(dfdu_,4,7) = 0.0; ARRAY2D(dfdu_,4,8) = 0.0; ARRAY2D(dfdu_,4,9) = 0.0; ARRAY2D(dfdu_,4,10) = 0.0; ARRAY2D(dfdu_,4,11) = 0.0; ARRAY2D(dfdu_,4,12) = 0.0; ARRAY2D(dfdu_,4,13) = 0.0; ARRAY2D(dfdu_,4,14) = 0.0; ARRAY2D(dfdu_,4,15) = 0.0; ARRAY2D(dfdu_,4,16) = 0.0; ARRAY2D(dfdu_,4,17) = 0.0; ARRAY2D(dfdu_,4,18) = 0.0; ARRAY2D(dfdu_,4,19) = 0.0; ARRAY2D(dfdu_,4,20) = 0.0; ARRAY2D(dfdu_,4,21) = 0.0; ARRAY2D(dfdu_,4,22) = 0.0; ARRAY2D(dfdu_,4,23) = 0.0; ARRAY2D(dfdu_,4,24) = 0.0; ARRAY2D(dfdu_,4,25) = 0.0; ARRAY2D(dfdu_,4,26) = 0.0; ARRAY2D(dfdu_,4,27) = 0.0; ARRAY2D(dfdu_,4,28) = 0.0; ARRAY2D(dfdu_,4,29) = 0.0; ARRAY2D(dfdu_,4,30) = 0.0; ARRAY2D(dfdu_,4,31) = 0.0; ARRAY2D(dfdu_,4,32) = 0.0; ARRAY2D(dfdu_,4,33) = 0.0; ARRAY2D(dfdu_,4,34) = 0.0; ARRAY2D(dfdu_,4,35) = 0.0; ARRAY2D(dfdu_,4,36) = 0.0; ARRAY2D(dfdu_,4,37) = 0.0; ARRAY2D(dfdu_,4,38) = 0.0; ARRAY2D(dfdu_,4,39) = 0.0; ARRAY2D(dfdu_,4,40) = 0.0; ARRAY2D(dfdu_,5,0) = 0.0; ARRAY2D(dfdu_,5,1) = 0.0; ARRAY2D(dfdu_,5,2) = 0.0; ARRAY2D(dfdu_,5,3) = 0.0; ARRAY2D(dfdu_,5,4) = 0.0; ARRAY2D(dfdu_,5,5) = 0.0; ARRAY2D(dfdu_,5,6) = 1.0; ARRAY2D(dfdu_,5,7) = 0.0; ARRAY2D(dfdu_,5,8) = 0.0; ARRAY2D(dfdu_,5,9) = 0.0; ARRAY2D(dfdu_,5,10) = 0.0; ARRAY2D(dfdu_,5,11) = 0.0; ARRAY2D(dfdu_,5,12) = 0.0; ARRAY2D(dfdu_,5,13) = 0.0; ARRAY2D(dfdu_,5,14) = 0.0; ARRAY2D(dfdu_,5,15) = 0.0; ARRAY2D(dfdu_,5,16) = 0.0; ARRAY2D(dfdu_,5,17) = 0.0; ARRAY2D(dfdu_,5,18) = 0.0; ARRAY2D(dfdu_,5,19) = 0.0; ARRAY2D(dfdu_,5,20) = 0.0; ARRAY2D(dfdu_,5,21) = 0.0; ARRAY2D(dfdu_,5,22) = 0.0; ARRAY2D(dfdu_,5,23) = 0.0; ARRAY2D(dfdu_,5,24) = 0.0; ARRAY2D(dfdu_,5,25) = 0.0; ARRAY2D(dfdu_,5,26) = 0.0; ARRAY2D(dfdu_,5,27) = 0.0; ARRAY2D(dfdu_,5,28) = 0.0; ARRAY2D(dfdu_,5,29) = 0.0; ARRAY2D(dfdu_,5,30) = 0.0; ARRAY2D(dfdu_,5,31) = 0.0; ARRAY2D(dfdu_,5,32) = 0.0; ARRAY2D(dfdu_,5,33) = 0.0; ARRAY2D(dfdu_,5,34) = 0.0; ARRAY2D(dfdu_,5,35) = 0.0; ARRAY2D(dfdu_,5,36) = 0.0; ARRAY2D(dfdu_,5,37) = 0.0; ARRAY2D(dfdu_,5,38) = 0.0; ARRAY2D(dfdu_,5,39) = 0.0; ARRAY2D(dfdu_,5,40) = 0.0; ARRAY2D(dfdu_,6,0) = 0.0; ARRAY2D(dfdu_,6,1) = 0.0; ARRAY2D(dfdu_,6,2) = 0.0; ARRAY2D(dfdu_,6,3) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,6,4) = 0.0; ARRAY2D(dfdu_,6,5) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,6,6) = -40.0*1.0/tau; ARRAY2D(dfdu_,6,7) = 0.0; ARRAY2D(dfdu_,6,8) = 0.0; ARRAY2D(dfdu_,6,9) = 0.0; ARRAY2D(dfdu_,6,10) = 0.0; ARRAY2D(dfdu_,6,11) = 0.0; ARRAY2D(dfdu_,6,12) = 0.0; ARRAY2D(dfdu_,6,13) = 0.0; ARRAY2D(dfdu_,6,14) = 0.0; ARRAY2D(dfdu_,6,15) = 0.0; ARRAY2D(dfdu_,6,16) = 0.0; ARRAY2D(dfdu_,6,17) = 0.0; ARRAY2D(dfdu_,6,18) = 0.0; ARRAY2D(dfdu_,6,19) = 0.0; ARRAY2D(dfdu_,6,20) = 0.0; ARRAY2D(dfdu_,6,21) = 0.0; ARRAY2D(dfdu_,6,22) = 0.0; ARRAY2D(dfdu_,6,23) = 0.0; ARRAY2D(dfdu_,6,24) = 0.0; ARRAY2D(dfdu_,6,25) = 0.0; ARRAY2D(dfdu_,6,26) = 0.0; ARRAY2D(dfdu_,6,27) = 0.0; ARRAY2D(dfdu_,6,28) = 0.0; ARRAY2D(dfdu_,6,29) = 0.0; ARRAY2D(dfdu_,6,30) = 0.0; ARRAY2D(dfdu_,6,31) = 0.0; ARRAY2D(dfdu_,6,32) = 0.0; ARRAY2D(dfdu_,6,33) = 0.0; ARRAY2D(dfdu_,6,34) = 0.0; ARRAY2D(dfdu_,6,35) = 0.0; ARRAY2D(dfdu_,6,36) = 0.0; ARRAY2D(dfdu_,6,37) = 0.0; ARRAY2D(dfdu_,6,38) = 0.0; ARRAY2D(dfdu_,6,39) = 0.0; ARRAY2D(dfdu_,6,40) = 0.0; ARRAY2D(dfdu_,7,0) = 0.0; ARRAY2D(dfdu_,7,1) = 0.0; ARRAY2D(dfdu_,7,2) = 0.0; ARRAY2D(dfdu_,7,3) = 0.0; ARRAY2D(dfdu_,7,4) = 0.0; ARRAY2D(dfdu_,7,5) = 0.0; ARRAY2D(dfdu_,7,6) = 0.0; ARRAY2D(dfdu_,7,7) = 0.0; ARRAY2D(dfdu_,7,8) = 1.0; ARRAY2D(dfdu_,7,9) = 0.0; ARRAY2D(dfdu_,7,10) = 0.0; ARRAY2D(dfdu_,7,11) = 0.0; ARRAY2D(dfdu_,7,12) = 0.0; ARRAY2D(dfdu_,7,13) = 0.0; ARRAY2D(dfdu_,7,14) = 0.0; ARRAY2D(dfdu_,7,15) = 0.0; ARRAY2D(dfdu_,7,16) = 0.0; ARRAY2D(dfdu_,7,17) = 0.0; ARRAY2D(dfdu_,7,18) = 0.0; ARRAY2D(dfdu_,7,19) = 0.0; ARRAY2D(dfdu_,7,20) = 0.0; ARRAY2D(dfdu_,7,21) = 0.0; ARRAY2D(dfdu_,7,22) = 0.0; ARRAY2D(dfdu_,7,23) = 0.0; ARRAY2D(dfdu_,7,24) = 0.0; ARRAY2D(dfdu_,7,25) = 0.0; ARRAY2D(dfdu_,7,26) = 0.0; ARRAY2D(dfdu_,7,27) = 0.0; ARRAY2D(dfdu_,7,28) = 0.0; ARRAY2D(dfdu_,7,29) = 0.0; ARRAY2D(dfdu_,7,30) = 0.0; ARRAY2D(dfdu_,7,31) = 0.0; ARRAY2D(dfdu_,7,32) = 0.0; ARRAY2D(dfdu_,7,33) = 0.0; ARRAY2D(dfdu_,7,34) = 0.0; ARRAY2D(dfdu_,7,35) = 0.0; ARRAY2D(dfdu_,7,36) = 0.0; ARRAY2D(dfdu_,7,37) = 0.0; ARRAY2D(dfdu_,7,38) = 0.0; ARRAY2D(dfdu_,7,39) = 0.0; ARRAY2D(dfdu_,7,40) = 0.0; ARRAY2D(dfdu_,8,0) = 0.0; ARRAY2D(dfdu_,8,1) = 0.0; ARRAY2D(dfdu_,8,2) = 0.0; ARRAY2D(dfdu_,8,3) = 0.0; ARRAY2D(dfdu_,8,4) = 0.0; ARRAY2D(dfdu_,8,5) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,8,6) = 0.0; ARRAY2D(dfdu_,8,7) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,8,8) = -40.0*1.0/tau; ARRAY2D(dfdu_,8,9) = 0.0; ARRAY2D(dfdu_,8,10) = 0.0; ARRAY2D(dfdu_,8,11) = 0.0; ARRAY2D(dfdu_,8,12) = 0.0; ARRAY2D(dfdu_,8,13) = 0.0; ARRAY2D(dfdu_,8,14) = 0.0; ARRAY2D(dfdu_,8,15) = 0.0; ARRAY2D(dfdu_,8,16) = 0.0; ARRAY2D(dfdu_,8,17) = 0.0; ARRAY2D(dfdu_,8,18) = 0.0; ARRAY2D(dfdu_,8,19) = 0.0; ARRAY2D(dfdu_,8,20) = 0.0; ARRAY2D(dfdu_,8,21) = 0.0; ARRAY2D(dfdu_,8,22) = 0.0; ARRAY2D(dfdu_,8,23) = 0.0; ARRAY2D(dfdu_,8,24) = 0.0; ARRAY2D(dfdu_,8,25) = 0.0; ARRAY2D(dfdu_,8,26) = 0.0; ARRAY2D(dfdu_,8,27) = 0.0; ARRAY2D(dfdu_,8,28) = 0.0; ARRAY2D(dfdu_,8,29) = 0.0; ARRAY2D(dfdu_,8,30) = 0.0; ARRAY2D(dfdu_,8,31) = 0.0; ARRAY2D(dfdu_,8,32) = 0.0; ARRAY2D(dfdu_,8,33) = 0.0; ARRAY2D(dfdu_,8,34) = 0.0; ARRAY2D(dfdu_,8,35) = 0.0; ARRAY2D(dfdu_,8,36) = 0.0; ARRAY2D(dfdu_,8,37) = 0.0; ARRAY2D(dfdu_,8,38) = 0.0; ARRAY2D(dfdu_,8,39) = 0.0; ARRAY2D(dfdu_,8,40) = 0.0; ARRAY2D(dfdu_,9,0) = 0.0; ARRAY2D(dfdu_,9,1) = 0.0; ARRAY2D(dfdu_,9,2) = 0.0; ARRAY2D(dfdu_,9,3) = 0.0; ARRAY2D(dfdu_,9,4) = 0.0; ARRAY2D(dfdu_,9,5) = 0.0; ARRAY2D(dfdu_,9,6) = 0.0; ARRAY2D(dfdu_,9,7) = 0.0; ARRAY2D(dfdu_,9,8) = 0.0; ARRAY2D(dfdu_,9,9) = 0.0; ARRAY2D(dfdu_,9,10) = 1.0; ARRAY2D(dfdu_,9,11) = 0.0; ARRAY2D(dfdu_,9,12) = 0.0; ARRAY2D(dfdu_,9,13) = 0.0; ARRAY2D(dfdu_,9,14) = 0.0; ARRAY2D(dfdu_,9,15) = 0.0; ARRAY2D(dfdu_,9,16) = 0.0; ARRAY2D(dfdu_,9,17) = 0.0; ARRAY2D(dfdu_,9,18) = 0.0; ARRAY2D(dfdu_,9,19) = 0.0; ARRAY2D(dfdu_,9,20) = 0.0; ARRAY2D(dfdu_,9,21) = 0.0; ARRAY2D(dfdu_,9,22) = 0.0; ARRAY2D(dfdu_,9,23) = 0.0; ARRAY2D(dfdu_,9,24) = 0.0; ARRAY2D(dfdu_,9,25) = 0.0; ARRAY2D(dfdu_,9,26) = 0.0; ARRAY2D(dfdu_,9,27) = 0.0; ARRAY2D(dfdu_,9,28) = 0.0; ARRAY2D(dfdu_,9,29) = 0.0; ARRAY2D(dfdu_,9,30) = 0.0; ARRAY2D(dfdu_,9,31) = 0.0; ARRAY2D(dfdu_,9,32) = 0.0; ARRAY2D(dfdu_,9,33) = 0.0; ARRAY2D(dfdu_,9,34) = 0.0; ARRAY2D(dfdu_,9,35) = 0.0; ARRAY2D(dfdu_,9,36) = 0.0; ARRAY2D(dfdu_,9,37) = 0.0; ARRAY2D(dfdu_,9,38) = 0.0; ARRAY2D(dfdu_,9,39) = 0.0; ARRAY2D(dfdu_,9,40) = 0.0; ARRAY2D(dfdu_,10,0) = 0.0; ARRAY2D(dfdu_,10,1) = 0.0; ARRAY2D(dfdu_,10,2) = 0.0; ARRAY2D(dfdu_,10,3) = 0.0; ARRAY2D(dfdu_,10,4) = 0.0; ARRAY2D(dfdu_,10,5) = 0.0; ARRAY2D(dfdu_,10,6) = 0.0; ARRAY2D(dfdu_,10,7) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,10,8) = 0.0; ARRAY2D(dfdu_,10,9) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,10,10) = -40.0*1.0/tau; ARRAY2D(dfdu_,10,11) = 0.0; ARRAY2D(dfdu_,10,12) = 0.0; ARRAY2D(dfdu_,10,13) = 0.0; ARRAY2D(dfdu_,10,14) = 0.0; ARRAY2D(dfdu_,10,15) = 0.0; ARRAY2D(dfdu_,10,16) = 0.0; ARRAY2D(dfdu_,10,17) = 0.0; ARRAY2D(dfdu_,10,18) = 0.0; ARRAY2D(dfdu_,10,19) = 0.0; ARRAY2D(dfdu_,10,20) = 0.0; ARRAY2D(dfdu_,10,21) = 0.0; ARRAY2D(dfdu_,10,22) = 0.0; ARRAY2D(dfdu_,10,23) = 0.0; ARRAY2D(dfdu_,10,24) = 0.0; ARRAY2D(dfdu_,10,25) = 0.0; ARRAY2D(dfdu_,10,26) = 0.0; ARRAY2D(dfdu_,10,27) = 0.0; ARRAY2D(dfdu_,10,28) = 0.0; ARRAY2D(dfdu_,10,29) = 0.0; ARRAY2D(dfdu_,10,30) = 0.0; ARRAY2D(dfdu_,10,31) = 0.0; ARRAY2D(dfdu_,10,32) = 0.0; ARRAY2D(dfdu_,10,33) = 0.0; ARRAY2D(dfdu_,10,34) = 0.0; ARRAY2D(dfdu_,10,35) = 0.0; ARRAY2D(dfdu_,10,36) = 0.0; ARRAY2D(dfdu_,10,37) = 0.0; ARRAY2D(dfdu_,10,38) = 0.0; ARRAY2D(dfdu_,10,39) = 0.0; ARRAY2D(dfdu_,10,40) = 0.0; ARRAY2D(dfdu_,11,0) = 0.0; ARRAY2D(dfdu_,11,1) = 0.0; ARRAY2D(dfdu_,11,2) = 0.0; ARRAY2D(dfdu_,11,3) = 0.0; ARRAY2D(dfdu_,11,4) = 0.0; ARRAY2D(dfdu_,11,5) = 0.0; ARRAY2D(dfdu_,11,6) = 0.0; ARRAY2D(dfdu_,11,7) = 0.0; ARRAY2D(dfdu_,11,8) = 0.0; ARRAY2D(dfdu_,11,9) = 0.0; ARRAY2D(dfdu_,11,10) = 0.0; ARRAY2D(dfdu_,11,11) = 0.0; ARRAY2D(dfdu_,11,12) = 1.0; ARRAY2D(dfdu_,11,13) = 0.0; ARRAY2D(dfdu_,11,14) = 0.0; ARRAY2D(dfdu_,11,15) = 0.0; ARRAY2D(dfdu_,11,16) = 0.0; ARRAY2D(dfdu_,11,17) = 0.0; ARRAY2D(dfdu_,11,18) = 0.0; ARRAY2D(dfdu_,11,19) = 0.0; ARRAY2D(dfdu_,11,20) = 0.0; ARRAY2D(dfdu_,11,21) = 0.0; ARRAY2D(dfdu_,11,22) = 0.0; ARRAY2D(dfdu_,11,23) = 0.0; ARRAY2D(dfdu_,11,24) = 0.0; ARRAY2D(dfdu_,11,25) = 0.0; ARRAY2D(dfdu_,11,26) = 0.0; ARRAY2D(dfdu_,11,27) = 0.0; ARRAY2D(dfdu_,11,28) = 0.0; ARRAY2D(dfdu_,11,29) = 0.0; ARRAY2D(dfdu_,11,30) = 0.0; ARRAY2D(dfdu_,11,31) = 0.0; ARRAY2D(dfdu_,11,32) = 0.0; ARRAY2D(dfdu_,11,33) = 0.0; ARRAY2D(dfdu_,11,34) = 0.0; ARRAY2D(dfdu_,11,35) = 0.0; ARRAY2D(dfdu_,11,36) = 0.0; ARRAY2D(dfdu_,11,37) = 0.0; ARRAY2D(dfdu_,11,38) = 0.0; ARRAY2D(dfdu_,11,39) = 0.0; ARRAY2D(dfdu_,11,40) = 0.0; ARRAY2D(dfdu_,12,0) = 0.0; ARRAY2D(dfdu_,12,1) = 0.0; ARRAY2D(dfdu_,12,2) = 0.0; ARRAY2D(dfdu_,12,3) = 0.0; ARRAY2D(dfdu_,12,4) = 0.0; ARRAY2D(dfdu_,12,5) = 0.0; ARRAY2D(dfdu_,12,6) = 0.0; ARRAY2D(dfdu_,12,7) = 0.0; ARRAY2D(dfdu_,12,8) = 0.0; ARRAY2D(dfdu_,12,9) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,12,10) = 0.0; ARRAY2D(dfdu_,12,11) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,12,12) = -40.0*1.0/tau; ARRAY2D(dfdu_,12,13) = 0.0; ARRAY2D(dfdu_,12,14) = 0.0; ARRAY2D(dfdu_,12,15) = 0.0; ARRAY2D(dfdu_,12,16) = 0.0; ARRAY2D(dfdu_,12,17) = 0.0; ARRAY2D(dfdu_,12,18) = 0.0; ARRAY2D(dfdu_,12,19) = 0.0; ARRAY2D(dfdu_,12,20) = 0.0; ARRAY2D(dfdu_,12,21) = 0.0; ARRAY2D(dfdu_,12,22) = 0.0; ARRAY2D(dfdu_,12,23) = 0.0; ARRAY2D(dfdu_,12,24) = 0.0; ARRAY2D(dfdu_,12,25) = 0.0; ARRAY2D(dfdu_,12,26) = 0.0; ARRAY2D(dfdu_,12,27) = 0.0; ARRAY2D(dfdu_,12,28) = 0.0; ARRAY2D(dfdu_,12,29) = 0.0; ARRAY2D(dfdu_,12,30) = 0.0; ARRAY2D(dfdu_,12,31) = 0.0; ARRAY2D(dfdu_,12,32) = 0.0; ARRAY2D(dfdu_,12,33) = 0.0; ARRAY2D(dfdu_,12,34) = 0.0; ARRAY2D(dfdu_,12,35) = 0.0; ARRAY2D(dfdu_,12,36) = 0.0; ARRAY2D(dfdu_,12,37) = 0.0; ARRAY2D(dfdu_,12,38) = 0.0; ARRAY2D(dfdu_,12,39) = 0.0; ARRAY2D(dfdu_,12,40) = 0.0; ARRAY2D(dfdu_,13,0) = 0.0; ARRAY2D(dfdu_,13,1) = 0.0; ARRAY2D(dfdu_,13,2) = 0.0; ARRAY2D(dfdu_,13,3) = 0.0; ARRAY2D(dfdu_,13,4) = 0.0; ARRAY2D(dfdu_,13,5) = 0.0; ARRAY2D(dfdu_,13,6) = 0.0; ARRAY2D(dfdu_,13,7) = 0.0; ARRAY2D(dfdu_,13,8) = 0.0; ARRAY2D(dfdu_,13,9) = 0.0; ARRAY2D(dfdu_,13,10) = 0.0; ARRAY2D(dfdu_,13,11) = 0.0; ARRAY2D(dfdu_,13,12) = 0.0; ARRAY2D(dfdu_,13,13) = 0.0; ARRAY2D(dfdu_,13,14) = 1.0; ARRAY2D(dfdu_,13,15) = 0.0; ARRAY2D(dfdu_,13,16) = 0.0; ARRAY2D(dfdu_,13,17) = 0.0; ARRAY2D(dfdu_,13,18) = 0.0; ARRAY2D(dfdu_,13,19) = 0.0; ARRAY2D(dfdu_,13,20) = 0.0; ARRAY2D(dfdu_,13,21) = 0.0; ARRAY2D(dfdu_,13,22) = 0.0; ARRAY2D(dfdu_,13,23) = 0.0; ARRAY2D(dfdu_,13,24) = 0.0; ARRAY2D(dfdu_,13,25) = 0.0; ARRAY2D(dfdu_,13,26) = 0.0; ARRAY2D(dfdu_,13,27) = 0.0; ARRAY2D(dfdu_,13,28) = 0.0; ARRAY2D(dfdu_,13,29) = 0.0; ARRAY2D(dfdu_,13,30) = 0.0; ARRAY2D(dfdu_,13,31) = 0.0; ARRAY2D(dfdu_,13,32) = 0.0; ARRAY2D(dfdu_,13,33) = 0.0; ARRAY2D(dfdu_,13,34) = 0.0; ARRAY2D(dfdu_,13,35) = 0.0; ARRAY2D(dfdu_,13,36) = 0.0; ARRAY2D(dfdu_,13,37) = 0.0; ARRAY2D(dfdu_,13,38) = 0.0; ARRAY2D(dfdu_,13,39) = 0.0; ARRAY2D(dfdu_,13,40) = 0.0; ARRAY2D(dfdu_,14,0) = 0.0; ARRAY2D(dfdu_,14,1) = 0.0; ARRAY2D(dfdu_,14,2) = 0.0; ARRAY2D(dfdu_,14,3) = 0.0; ARRAY2D(dfdu_,14,4) = 0.0; ARRAY2D(dfdu_,14,5) = 0.0; ARRAY2D(dfdu_,14,6) = 0.0; ARRAY2D(dfdu_,14,7) = 0.0; ARRAY2D(dfdu_,14,8) = 0.0; ARRAY2D(dfdu_,14,9) = 0.0; ARRAY2D(dfdu_,14,10) = 0.0; ARRAY2D(dfdu_,14,11) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,14,12) = 0.0; ARRAY2D(dfdu_,14,13) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,14,14) = -40.0*1.0/tau; ARRAY2D(dfdu_,14,15) = 0.0; ARRAY2D(dfdu_,14,16) = 0.0; ARRAY2D(dfdu_,14,17) = 0.0; ARRAY2D(dfdu_,14,18) = 0.0; ARRAY2D(dfdu_,14,19) = 0.0; ARRAY2D(dfdu_,14,20) = 0.0; ARRAY2D(dfdu_,14,21) = 0.0; ARRAY2D(dfdu_,14,22) = 0.0; ARRAY2D(dfdu_,14,23) = 0.0; ARRAY2D(dfdu_,14,24) = 0.0; ARRAY2D(dfdu_,14,25) = 0.0; ARRAY2D(dfdu_,14,26) = 0.0; ARRAY2D(dfdu_,14,27) = 0.0; ARRAY2D(dfdu_,14,28) = 0.0; ARRAY2D(dfdu_,14,29) = 0.0; ARRAY2D(dfdu_,14,30) = 0.0; ARRAY2D(dfdu_,14,31) = 0.0; ARRAY2D(dfdu_,14,32) = 0.0; ARRAY2D(dfdu_,14,33) = 0.0; ARRAY2D(dfdu_,14,34) = 0.0; ARRAY2D(dfdu_,14,35) = 0.0; ARRAY2D(dfdu_,14,36) = 0.0; ARRAY2D(dfdu_,14,37) = 0.0; ARRAY2D(dfdu_,14,38) = 0.0; ARRAY2D(dfdu_,14,39) = 0.0; ARRAY2D(dfdu_,14,40) = 0.0; ARRAY2D(dfdu_,15,0) = 0.0; ARRAY2D(dfdu_,15,1) = 0.0; ARRAY2D(dfdu_,15,2) = 0.0; ARRAY2D(dfdu_,15,3) = 0.0; ARRAY2D(dfdu_,15,4) = 0.0; ARRAY2D(dfdu_,15,5) = 0.0; ARRAY2D(dfdu_,15,6) = 0.0; ARRAY2D(dfdu_,15,7) = 0.0; ARRAY2D(dfdu_,15,8) = 0.0; ARRAY2D(dfdu_,15,9) = 0.0; ARRAY2D(dfdu_,15,10) = 0.0; ARRAY2D(dfdu_,15,11) = 0.0; ARRAY2D(dfdu_,15,12) = 0.0; ARRAY2D(dfdu_,15,13) = 0.0; ARRAY2D(dfdu_,15,14) = 0.0; ARRAY2D(dfdu_,15,15) = 0.0; ARRAY2D(dfdu_,15,16) = 1.0; ARRAY2D(dfdu_,15,17) = 0.0; ARRAY2D(dfdu_,15,18) = 0.0; ARRAY2D(dfdu_,15,19) = 0.0; ARRAY2D(dfdu_,15,20) = 0.0; ARRAY2D(dfdu_,15,21) = 0.0; ARRAY2D(dfdu_,15,22) = 0.0; ARRAY2D(dfdu_,15,23) = 0.0; ARRAY2D(dfdu_,15,24) = 0.0; ARRAY2D(dfdu_,15,25) = 0.0; ARRAY2D(dfdu_,15,26) = 0.0; ARRAY2D(dfdu_,15,27) = 0.0; ARRAY2D(dfdu_,15,28) = 0.0; ARRAY2D(dfdu_,15,29) = 0.0; ARRAY2D(dfdu_,15,30) = 0.0; ARRAY2D(dfdu_,15,31) = 0.0; ARRAY2D(dfdu_,15,32) = 0.0; ARRAY2D(dfdu_,15,33) = 0.0; ARRAY2D(dfdu_,15,34) = 0.0; ARRAY2D(dfdu_,15,35) = 0.0; ARRAY2D(dfdu_,15,36) = 0.0; ARRAY2D(dfdu_,15,37) = 0.0; ARRAY2D(dfdu_,15,38) = 0.0; ARRAY2D(dfdu_,15,39) = 0.0; ARRAY2D(dfdu_,15,40) = 0.0; ARRAY2D(dfdu_,16,0) = 0.0; ARRAY2D(dfdu_,16,1) = 0.0; ARRAY2D(dfdu_,16,2) = 0.0; ARRAY2D(dfdu_,16,3) = 0.0; ARRAY2D(dfdu_,16,4) = 0.0; ARRAY2D(dfdu_,16,5) = 0.0; ARRAY2D(dfdu_,16,6) = 0.0; ARRAY2D(dfdu_,16,7) = 0.0; ARRAY2D(dfdu_,16,8) = 0.0; ARRAY2D(dfdu_,16,9) = 0.0; ARRAY2D(dfdu_,16,10) = 0.0; ARRAY2D(dfdu_,16,11) = 0.0; ARRAY2D(dfdu_,16,12) = 0.0; ARRAY2D(dfdu_,16,13) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,16,14) = 0.0; ARRAY2D(dfdu_,16,15) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,16,16) = -40.0*1.0/tau; ARRAY2D(dfdu_,16,17) = 0.0; ARRAY2D(dfdu_,16,18) = 0.0; ARRAY2D(dfdu_,16,19) = 0.0; ARRAY2D(dfdu_,16,20) = 0.0; ARRAY2D(dfdu_,16,21) = 0.0; ARRAY2D(dfdu_,16,22) = 0.0; ARRAY2D(dfdu_,16,23) = 0.0; ARRAY2D(dfdu_,16,24) = 0.0; ARRAY2D(dfdu_,16,25) = 0.0; ARRAY2D(dfdu_,16,26) = 0.0; ARRAY2D(dfdu_,16,27) = 0.0; ARRAY2D(dfdu_,16,28) = 0.0; ARRAY2D(dfdu_,16,29) = 0.0; ARRAY2D(dfdu_,16,30) = 0.0; ARRAY2D(dfdu_,16,31) = 0.0; ARRAY2D(dfdu_,16,32) = 0.0; ARRAY2D(dfdu_,16,33) = 0.0; ARRAY2D(dfdu_,16,34) = 0.0; ARRAY2D(dfdu_,16,35) = 0.0; ARRAY2D(dfdu_,16,36) = 0.0; ARRAY2D(dfdu_,16,37) = 0.0; ARRAY2D(dfdu_,16,38) = 0.0; ARRAY2D(dfdu_,16,39) = 0.0; ARRAY2D(dfdu_,16,40) = 0.0; ARRAY2D(dfdu_,17,0) = 0.0; ARRAY2D(dfdu_,17,1) = 0.0; ARRAY2D(dfdu_,17,2) = 0.0; ARRAY2D(dfdu_,17,3) = 0.0; ARRAY2D(dfdu_,17,4) = 0.0; ARRAY2D(dfdu_,17,5) = 0.0; ARRAY2D(dfdu_,17,6) = 0.0; ARRAY2D(dfdu_,17,7) = 0.0; ARRAY2D(dfdu_,17,8) = 0.0; ARRAY2D(dfdu_,17,9) = 0.0; ARRAY2D(dfdu_,17,10) = 0.0; ARRAY2D(dfdu_,17,11) = 0.0; ARRAY2D(dfdu_,17,12) = 0.0; ARRAY2D(dfdu_,17,13) = 0.0; ARRAY2D(dfdu_,17,14) = 0.0; ARRAY2D(dfdu_,17,15) = 0.0; ARRAY2D(dfdu_,17,16) = 0.0; ARRAY2D(dfdu_,17,17) = 0.0; ARRAY2D(dfdu_,17,18) = 1.0; ARRAY2D(dfdu_,17,19) = 0.0; ARRAY2D(dfdu_,17,20) = 0.0; ARRAY2D(dfdu_,17,21) = 0.0; ARRAY2D(dfdu_,17,22) = 0.0; ARRAY2D(dfdu_,17,23) = 0.0; ARRAY2D(dfdu_,17,24) = 0.0; ARRAY2D(dfdu_,17,25) = 0.0; ARRAY2D(dfdu_,17,26) = 0.0; ARRAY2D(dfdu_,17,27) = 0.0; ARRAY2D(dfdu_,17,28) = 0.0; ARRAY2D(dfdu_,17,29) = 0.0; ARRAY2D(dfdu_,17,30) = 0.0; ARRAY2D(dfdu_,17,31) = 0.0; ARRAY2D(dfdu_,17,32) = 0.0; ARRAY2D(dfdu_,17,33) = 0.0; ARRAY2D(dfdu_,17,34) = 0.0; ARRAY2D(dfdu_,17,35) = 0.0; ARRAY2D(dfdu_,17,36) = 0.0; ARRAY2D(dfdu_,17,37) = 0.0; ARRAY2D(dfdu_,17,38) = 0.0; ARRAY2D(dfdu_,17,39) = 0.0; ARRAY2D(dfdu_,17,40) = 0.0; ARRAY2D(dfdu_,18,0) = 0.0; ARRAY2D(dfdu_,18,1) = 0.0; ARRAY2D(dfdu_,18,2) = 0.0; ARRAY2D(dfdu_,18,3) = 0.0; ARRAY2D(dfdu_,18,4) = 0.0; ARRAY2D(dfdu_,18,5) = 0.0; ARRAY2D(dfdu_,18,6) = 0.0; ARRAY2D(dfdu_,18,7) = 0.0; ARRAY2D(dfdu_,18,8) = 0.0; ARRAY2D(dfdu_,18,9) = 0.0; ARRAY2D(dfdu_,18,10) = 0.0; ARRAY2D(dfdu_,18,11) = 0.0; ARRAY2D(dfdu_,18,12) = 0.0; ARRAY2D(dfdu_,18,13) = 0.0; ARRAY2D(dfdu_,18,14) = 0.0; ARRAY2D(dfdu_,18,15) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,18,16) = 0.0; ARRAY2D(dfdu_,18,17) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,18,18) = -40.0*1.0/tau; ARRAY2D(dfdu_,18,19) = 0.0; ARRAY2D(dfdu_,18,20) = 0.0; ARRAY2D(dfdu_,18,21) = 0.0; ARRAY2D(dfdu_,18,22) = 0.0; ARRAY2D(dfdu_,18,23) = 0.0; ARRAY2D(dfdu_,18,24) = 0.0; ARRAY2D(dfdu_,18,25) = 0.0; ARRAY2D(dfdu_,18,26) = 0.0; ARRAY2D(dfdu_,18,27) = 0.0; ARRAY2D(dfdu_,18,28) = 0.0; ARRAY2D(dfdu_,18,29) = 0.0; ARRAY2D(dfdu_,18,30) = 0.0; ARRAY2D(dfdu_,18,31) = 0.0; ARRAY2D(dfdu_,18,32) = 0.0; ARRAY2D(dfdu_,18,33) = 0.0; ARRAY2D(dfdu_,18,34) = 0.0; ARRAY2D(dfdu_,18,35) = 0.0; ARRAY2D(dfdu_,18,36) = 0.0; ARRAY2D(dfdu_,18,37) = 0.0; ARRAY2D(dfdu_,18,38) = 0.0; ARRAY2D(dfdu_,18,39) = 0.0; ARRAY2D(dfdu_,18,40) = 0.0; ARRAY2D(dfdu_,19,0) = 0.0; ARRAY2D(dfdu_,19,1) = 0.0; ARRAY2D(dfdu_,19,2) = 0.0; ARRAY2D(dfdu_,19,3) = 0.0; ARRAY2D(dfdu_,19,4) = 0.0; ARRAY2D(dfdu_,19,5) = 0.0; ARRAY2D(dfdu_,19,6) = 0.0; ARRAY2D(dfdu_,19,7) = 0.0; ARRAY2D(dfdu_,19,8) = 0.0; ARRAY2D(dfdu_,19,9) = 0.0; ARRAY2D(dfdu_,19,10) = 0.0; ARRAY2D(dfdu_,19,11) = 0.0; ARRAY2D(dfdu_,19,12) = 0.0; ARRAY2D(dfdu_,19,13) = 0.0; ARRAY2D(dfdu_,19,14) = 0.0; ARRAY2D(dfdu_,19,15) = 0.0; ARRAY2D(dfdu_,19,16) = 0.0; ARRAY2D(dfdu_,19,17) = 0.0; ARRAY2D(dfdu_,19,18) = 0.0; ARRAY2D(dfdu_,19,19) = 0.0; ARRAY2D(dfdu_,19,20) = 1.0; ARRAY2D(dfdu_,19,21) = 0.0; ARRAY2D(dfdu_,19,22) = 0.0; ARRAY2D(dfdu_,19,23) = 0.0; ARRAY2D(dfdu_,19,24) = 0.0; ARRAY2D(dfdu_,19,25) = 0.0; ARRAY2D(dfdu_,19,26) = 0.0; ARRAY2D(dfdu_,19,27) = 0.0; ARRAY2D(dfdu_,19,28) = 0.0; ARRAY2D(dfdu_,19,29) = 0.0; ARRAY2D(dfdu_,19,30) = 0.0; ARRAY2D(dfdu_,19,31) = 0.0; ARRAY2D(dfdu_,19,32) = 0.0; ARRAY2D(dfdu_,19,33) = 0.0; ARRAY2D(dfdu_,19,34) = 0.0; ARRAY2D(dfdu_,19,35) = 0.0; ARRAY2D(dfdu_,19,36) = 0.0; ARRAY2D(dfdu_,19,37) = 0.0; ARRAY2D(dfdu_,19,38) = 0.0; ARRAY2D(dfdu_,19,39) = 0.0; ARRAY2D(dfdu_,19,40) = 0.0; ARRAY2D(dfdu_,20,0) = 0.0; ARRAY2D(dfdu_,20,1) = 0.0; ARRAY2D(dfdu_,20,2) = 0.0; ARRAY2D(dfdu_,20,3) = 0.0; ARRAY2D(dfdu_,20,4) = 0.0; ARRAY2D(dfdu_,20,5) = 0.0; ARRAY2D(dfdu_,20,6) = 0.0; ARRAY2D(dfdu_,20,7) = 0.0; ARRAY2D(dfdu_,20,8) = 0.0; ARRAY2D(dfdu_,20,9) = 0.0; ARRAY2D(dfdu_,20,10) = 0.0; ARRAY2D(dfdu_,20,11) = 0.0; ARRAY2D(dfdu_,20,12) = 0.0; ARRAY2D(dfdu_,20,13) = 0.0; ARRAY2D(dfdu_,20,14) = 0.0; ARRAY2D(dfdu_,20,15) = 0.0; ARRAY2D(dfdu_,20,16) = 0.0; ARRAY2D(dfdu_,20,17) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,20,18) = 0.0; ARRAY2D(dfdu_,20,19) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,20,20) = -40.0*1.0/tau; ARRAY2D(dfdu_,20,21) = 0.0; ARRAY2D(dfdu_,20,22) = 0.0; ARRAY2D(dfdu_,20,23) = 0.0; ARRAY2D(dfdu_,20,24) = 0.0; ARRAY2D(dfdu_,20,25) = 0.0; ARRAY2D(dfdu_,20,26) = 0.0; ARRAY2D(dfdu_,20,27) = 0.0; ARRAY2D(dfdu_,20,28) = 0.0; ARRAY2D(dfdu_,20,29) = 0.0; ARRAY2D(dfdu_,20,30) = 0.0; ARRAY2D(dfdu_,20,31) = 0.0; ARRAY2D(dfdu_,20,32) = 0.0; ARRAY2D(dfdu_,20,33) = 0.0; ARRAY2D(dfdu_,20,34) = 0.0; ARRAY2D(dfdu_,20,35) = 0.0; ARRAY2D(dfdu_,20,36) = 0.0; ARRAY2D(dfdu_,20,37) = 0.0; ARRAY2D(dfdu_,20,38) = 0.0; ARRAY2D(dfdu_,20,39) = 0.0; ARRAY2D(dfdu_,20,40) = 0.0; ARRAY2D(dfdu_,21,0) = 0.0; ARRAY2D(dfdu_,21,1) = 0.0; ARRAY2D(dfdu_,21,2) = 0.0; ARRAY2D(dfdu_,21,3) = 0.0; ARRAY2D(dfdu_,21,4) = 0.0; ARRAY2D(dfdu_,21,5) = 0.0; ARRAY2D(dfdu_,21,6) = 0.0; ARRAY2D(dfdu_,21,7) = 0.0; ARRAY2D(dfdu_,21,8) = 0.0; ARRAY2D(dfdu_,21,9) = 0.0; ARRAY2D(dfdu_,21,10) = 0.0; ARRAY2D(dfdu_,21,11) = 0.0; ARRAY2D(dfdu_,21,12) = 0.0; ARRAY2D(dfdu_,21,13) = 0.0; ARRAY2D(dfdu_,21,14) = 0.0; ARRAY2D(dfdu_,21,15) = 0.0; ARRAY2D(dfdu_,21,16) = 0.0; ARRAY2D(dfdu_,21,17) = 0.0; ARRAY2D(dfdu_,21,18) = 0.0; ARRAY2D(dfdu_,21,19) = 0.0; ARRAY2D(dfdu_,21,20) = 0.0; ARRAY2D(dfdu_,21,21) = 0.0; ARRAY2D(dfdu_,21,22) = 1.0; ARRAY2D(dfdu_,21,23) = 0.0; ARRAY2D(dfdu_,21,24) = 0.0; ARRAY2D(dfdu_,21,25) = 0.0; ARRAY2D(dfdu_,21,26) = 0.0; ARRAY2D(dfdu_,21,27) = 0.0; ARRAY2D(dfdu_,21,28) = 0.0; ARRAY2D(dfdu_,21,29) = 0.0; ARRAY2D(dfdu_,21,30) = 0.0; ARRAY2D(dfdu_,21,31) = 0.0; ARRAY2D(dfdu_,21,32) = 0.0; ARRAY2D(dfdu_,21,33) = 0.0; ARRAY2D(dfdu_,21,34) = 0.0; ARRAY2D(dfdu_,21,35) = 0.0; ARRAY2D(dfdu_,21,36) = 0.0; ARRAY2D(dfdu_,21,37) = 0.0; ARRAY2D(dfdu_,21,38) = 0.0; ARRAY2D(dfdu_,21,39) = 0.0; ARRAY2D(dfdu_,21,40) = 0.0; ARRAY2D(dfdu_,22,0) = 0.0; ARRAY2D(dfdu_,22,1) = 0.0; ARRAY2D(dfdu_,22,2) = 0.0; ARRAY2D(dfdu_,22,3) = 0.0; ARRAY2D(dfdu_,22,4) = 0.0; ARRAY2D(dfdu_,22,5) = 0.0; ARRAY2D(dfdu_,22,6) = 0.0; ARRAY2D(dfdu_,22,7) = 0.0; ARRAY2D(dfdu_,22,8) = 0.0; ARRAY2D(dfdu_,22,9) = 0.0; ARRAY2D(dfdu_,22,10) = 0.0; ARRAY2D(dfdu_,22,11) = 0.0; ARRAY2D(dfdu_,22,12) = 0.0; ARRAY2D(dfdu_,22,13) = 0.0; ARRAY2D(dfdu_,22,14) = 0.0; ARRAY2D(dfdu_,22,15) = 0.0; ARRAY2D(dfdu_,22,16) = 0.0; ARRAY2D(dfdu_,22,17) = 0.0; ARRAY2D(dfdu_,22,18) = 0.0; ARRAY2D(dfdu_,22,19) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,22,20) = 0.0; ARRAY2D(dfdu_,22,21) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,22,22) = -40.0*1.0/tau; ARRAY2D(dfdu_,22,23) = 0.0; ARRAY2D(dfdu_,22,24) = 0.0; ARRAY2D(dfdu_,22,25) = 0.0; ARRAY2D(dfdu_,22,26) = 0.0; ARRAY2D(dfdu_,22,27) = 0.0; ARRAY2D(dfdu_,22,28) = 0.0; ARRAY2D(dfdu_,22,29) = 0.0; ARRAY2D(dfdu_,22,30) = 0.0; ARRAY2D(dfdu_,22,31) = 0.0; ARRAY2D(dfdu_,22,32) = 0.0; ARRAY2D(dfdu_,22,33) = 0.0; ARRAY2D(dfdu_,22,34) = 0.0; ARRAY2D(dfdu_,22,35) = 0.0; ARRAY2D(dfdu_,22,36) = 0.0; ARRAY2D(dfdu_,22,37) = 0.0; ARRAY2D(dfdu_,22,38) = 0.0; ARRAY2D(dfdu_,22,39) = 0.0; ARRAY2D(dfdu_,22,40) = 0.0; ARRAY2D(dfdu_,23,0) = 0.0; ARRAY2D(dfdu_,23,1) = 0.0; ARRAY2D(dfdu_,23,2) = 0.0; ARRAY2D(dfdu_,23,3) = 0.0; ARRAY2D(dfdu_,23,4) = 0.0; ARRAY2D(dfdu_,23,5) = 0.0; ARRAY2D(dfdu_,23,6) = 0.0; ARRAY2D(dfdu_,23,7) = 0.0; ARRAY2D(dfdu_,23,8) = 0.0; ARRAY2D(dfdu_,23,9) = 0.0; ARRAY2D(dfdu_,23,10) = 0.0; ARRAY2D(dfdu_,23,11) = 0.0; ARRAY2D(dfdu_,23,12) = 0.0; ARRAY2D(dfdu_,23,13) = 0.0; ARRAY2D(dfdu_,23,14) = 0.0; ARRAY2D(dfdu_,23,15) = 0.0; ARRAY2D(dfdu_,23,16) = 0.0; ARRAY2D(dfdu_,23,17) = 0.0; ARRAY2D(dfdu_,23,18) = 0.0; ARRAY2D(dfdu_,23,19) = 0.0; ARRAY2D(dfdu_,23,20) = 0.0; ARRAY2D(dfdu_,23,21) = 0.0; ARRAY2D(dfdu_,23,22) = 0.0; ARRAY2D(dfdu_,23,23) = 0.0; ARRAY2D(dfdu_,23,24) = 1.0; ARRAY2D(dfdu_,23,25) = 0.0; ARRAY2D(dfdu_,23,26) = 0.0; ARRAY2D(dfdu_,23,27) = 0.0; ARRAY2D(dfdu_,23,28) = 0.0; ARRAY2D(dfdu_,23,29) = 0.0; ARRAY2D(dfdu_,23,30) = 0.0; ARRAY2D(dfdu_,23,31) = 0.0; ARRAY2D(dfdu_,23,32) = 0.0; ARRAY2D(dfdu_,23,33) = 0.0; ARRAY2D(dfdu_,23,34) = 0.0; ARRAY2D(dfdu_,23,35) = 0.0; ARRAY2D(dfdu_,23,36) = 0.0; ARRAY2D(dfdu_,23,37) = 0.0; ARRAY2D(dfdu_,23,38) = 0.0; ARRAY2D(dfdu_,23,39) = 0.0; ARRAY2D(dfdu_,23,40) = 0.0; ARRAY2D(dfdu_,24,0) = 0.0; ARRAY2D(dfdu_,24,1) = 0.0; ARRAY2D(dfdu_,24,2) = 0.0; ARRAY2D(dfdu_,24,3) = 0.0; ARRAY2D(dfdu_,24,4) = 0.0; ARRAY2D(dfdu_,24,5) = 0.0; ARRAY2D(dfdu_,24,6) = 0.0; ARRAY2D(dfdu_,24,7) = 0.0; ARRAY2D(dfdu_,24,8) = 0.0; ARRAY2D(dfdu_,24,9) = 0.0; ARRAY2D(dfdu_,24,10) = 0.0; ARRAY2D(dfdu_,24,11) = 0.0; ARRAY2D(dfdu_,24,12) = 0.0; ARRAY2D(dfdu_,24,13) = 0.0; ARRAY2D(dfdu_,24,14) = 0.0; ARRAY2D(dfdu_,24,15) = 0.0; ARRAY2D(dfdu_,24,16) = 0.0; ARRAY2D(dfdu_,24,17) = 0.0; ARRAY2D(dfdu_,24,18) = 0.0; ARRAY2D(dfdu_,24,19) = 0.0; ARRAY2D(dfdu_,24,20) = 0.0; ARRAY2D(dfdu_,24,21) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,24,22) = 0.0; ARRAY2D(dfdu_,24,23) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,24,24) = -40.0*1.0/tau; ARRAY2D(dfdu_,24,25) = 0.0; ARRAY2D(dfdu_,24,26) = 0.0; ARRAY2D(dfdu_,24,27) = 0.0; ARRAY2D(dfdu_,24,28) = 0.0; ARRAY2D(dfdu_,24,29) = 0.0; ARRAY2D(dfdu_,24,30) = 0.0; ARRAY2D(dfdu_,24,31) = 0.0; ARRAY2D(dfdu_,24,32) = 0.0; ARRAY2D(dfdu_,24,33) = 0.0; ARRAY2D(dfdu_,24,34) = 0.0; ARRAY2D(dfdu_,24,35) = 0.0; ARRAY2D(dfdu_,24,36) = 0.0; ARRAY2D(dfdu_,24,37) = 0.0; ARRAY2D(dfdu_,24,38) = 0.0; ARRAY2D(dfdu_,24,39) = 0.0; ARRAY2D(dfdu_,24,40) = 0.0; ARRAY2D(dfdu_,25,0) = 0.0; ARRAY2D(dfdu_,25,1) = 0.0; ARRAY2D(dfdu_,25,2) = 0.0; ARRAY2D(dfdu_,25,3) = 0.0; ARRAY2D(dfdu_,25,4) = 0.0; ARRAY2D(dfdu_,25,5) = 0.0; ARRAY2D(dfdu_,25,6) = 0.0; ARRAY2D(dfdu_,25,7) = 0.0; ARRAY2D(dfdu_,25,8) = 0.0; ARRAY2D(dfdu_,25,9) = 0.0; ARRAY2D(dfdu_,25,10) = 0.0; ARRAY2D(dfdu_,25,11) = 0.0; ARRAY2D(dfdu_,25,12) = 0.0; ARRAY2D(dfdu_,25,13) = 0.0; ARRAY2D(dfdu_,25,14) = 0.0; ARRAY2D(dfdu_,25,15) = 0.0; ARRAY2D(dfdu_,25,16) = 0.0; ARRAY2D(dfdu_,25,17) = 0.0; ARRAY2D(dfdu_,25,18) = 0.0; ARRAY2D(dfdu_,25,19) = 0.0; ARRAY2D(dfdu_,25,20) = 0.0; ARRAY2D(dfdu_,25,21) = 0.0; ARRAY2D(dfdu_,25,22) = 0.0; ARRAY2D(dfdu_,25,23) = 0.0; ARRAY2D(dfdu_,25,24) = 0.0; ARRAY2D(dfdu_,25,25) = 0.0; ARRAY2D(dfdu_,25,26) = 1.0; ARRAY2D(dfdu_,25,27) = 0.0; ARRAY2D(dfdu_,25,28) = 0.0; ARRAY2D(dfdu_,25,29) = 0.0; ARRAY2D(dfdu_,25,30) = 0.0; ARRAY2D(dfdu_,25,31) = 0.0; ARRAY2D(dfdu_,25,32) = 0.0; ARRAY2D(dfdu_,25,33) = 0.0; ARRAY2D(dfdu_,25,34) = 0.0; ARRAY2D(dfdu_,25,35) = 0.0; ARRAY2D(dfdu_,25,36) = 0.0; ARRAY2D(dfdu_,25,37) = 0.0; ARRAY2D(dfdu_,25,38) = 0.0; ARRAY2D(dfdu_,25,39) = 0.0; ARRAY2D(dfdu_,25,40) = 0.0; ARRAY2D(dfdu_,26,0) = 0.0; ARRAY2D(dfdu_,26,1) = 0.0; ARRAY2D(dfdu_,26,2) = 0.0; ARRAY2D(dfdu_,26,3) = 0.0; ARRAY2D(dfdu_,26,4) = 0.0; ARRAY2D(dfdu_,26,5) = 0.0; ARRAY2D(dfdu_,26,6) = 0.0; ARRAY2D(dfdu_,26,7) = 0.0; ARRAY2D(dfdu_,26,8) = 0.0; ARRAY2D(dfdu_,26,9) = 0.0; ARRAY2D(dfdu_,26,10) = 0.0; ARRAY2D(dfdu_,26,11) = 0.0; ARRAY2D(dfdu_,26,12) = 0.0; ARRAY2D(dfdu_,26,13) = 0.0; ARRAY2D(dfdu_,26,14) = 0.0; ARRAY2D(dfdu_,26,15) = 0.0; ARRAY2D(dfdu_,26,16) = 0.0; ARRAY2D(dfdu_,26,17) = 0.0; ARRAY2D(dfdu_,26,18) = 0.0; ARRAY2D(dfdu_,26,19) = 0.0; ARRAY2D(dfdu_,26,20) = 0.0; ARRAY2D(dfdu_,26,21) = 0.0; ARRAY2D(dfdu_,26,22) = 0.0; ARRAY2D(dfdu_,26,23) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,26,24) = 0.0; ARRAY2D(dfdu_,26,25) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,26,26) = -40.0*1.0/tau; ARRAY2D(dfdu_,26,27) = 0.0; ARRAY2D(dfdu_,26,28) = 0.0; ARRAY2D(dfdu_,26,29) = 0.0; ARRAY2D(dfdu_,26,30) = 0.0; ARRAY2D(dfdu_,26,31) = 0.0; ARRAY2D(dfdu_,26,32) = 0.0; ARRAY2D(dfdu_,26,33) = 0.0; ARRAY2D(dfdu_,26,34) = 0.0; ARRAY2D(dfdu_,26,35) = 0.0; ARRAY2D(dfdu_,26,36) = 0.0; ARRAY2D(dfdu_,26,37) = 0.0; ARRAY2D(dfdu_,26,38) = 0.0; ARRAY2D(dfdu_,26,39) = 0.0; ARRAY2D(dfdu_,26,40) = 0.0; ARRAY2D(dfdu_,27,0) = 0.0; ARRAY2D(dfdu_,27,1) = 0.0; ARRAY2D(dfdu_,27,2) = 0.0; ARRAY2D(dfdu_,27,3) = 0.0; ARRAY2D(dfdu_,27,4) = 0.0; ARRAY2D(dfdu_,27,5) = 0.0; ARRAY2D(dfdu_,27,6) = 0.0; ARRAY2D(dfdu_,27,7) = 0.0; ARRAY2D(dfdu_,27,8) = 0.0; ARRAY2D(dfdu_,27,9) = 0.0; ARRAY2D(dfdu_,27,10) = 0.0; ARRAY2D(dfdu_,27,11) = 0.0; ARRAY2D(dfdu_,27,12) = 0.0; ARRAY2D(dfdu_,27,13) = 0.0; ARRAY2D(dfdu_,27,14) = 0.0; ARRAY2D(dfdu_,27,15) = 0.0; ARRAY2D(dfdu_,27,16) = 0.0; ARRAY2D(dfdu_,27,17) = 0.0; ARRAY2D(dfdu_,27,18) = 0.0; ARRAY2D(dfdu_,27,19) = 0.0; ARRAY2D(dfdu_,27,20) = 0.0; ARRAY2D(dfdu_,27,21) = 0.0; ARRAY2D(dfdu_,27,22) = 0.0; ARRAY2D(dfdu_,27,23) = 0.0; ARRAY2D(dfdu_,27,24) = 0.0; ARRAY2D(dfdu_,27,25) = 0.0; ARRAY2D(dfdu_,27,26) = 0.0; ARRAY2D(dfdu_,27,27) = 0.0; ARRAY2D(dfdu_,27,28) = 1.0; ARRAY2D(dfdu_,27,29) = 0.0; ARRAY2D(dfdu_,27,30) = 0.0; ARRAY2D(dfdu_,27,31) = 0.0; ARRAY2D(dfdu_,27,32) = 0.0; ARRAY2D(dfdu_,27,33) = 0.0; ARRAY2D(dfdu_,27,34) = 0.0; ARRAY2D(dfdu_,27,35) = 0.0; ARRAY2D(dfdu_,27,36) = 0.0; ARRAY2D(dfdu_,27,37) = 0.0; ARRAY2D(dfdu_,27,38) = 0.0; ARRAY2D(dfdu_,27,39) = 0.0; ARRAY2D(dfdu_,27,40) = 0.0; ARRAY2D(dfdu_,28,0) = 0.0; ARRAY2D(dfdu_,28,1) = 0.0; ARRAY2D(dfdu_,28,2) = 0.0; ARRAY2D(dfdu_,28,3) = 0.0; ARRAY2D(dfdu_,28,4) = 0.0; ARRAY2D(dfdu_,28,5) = 0.0; ARRAY2D(dfdu_,28,6) = 0.0; ARRAY2D(dfdu_,28,7) = 0.0; ARRAY2D(dfdu_,28,8) = 0.0; ARRAY2D(dfdu_,28,9) = 0.0; ARRAY2D(dfdu_,28,10) = 0.0; ARRAY2D(dfdu_,28,11) = 0.0; ARRAY2D(dfdu_,28,12) = 0.0; ARRAY2D(dfdu_,28,13) = 0.0; ARRAY2D(dfdu_,28,14) = 0.0; ARRAY2D(dfdu_,28,15) = 0.0; ARRAY2D(dfdu_,28,16) = 0.0; ARRAY2D(dfdu_,28,17) = 0.0; ARRAY2D(dfdu_,28,18) = 0.0; ARRAY2D(dfdu_,28,19) = 0.0; ARRAY2D(dfdu_,28,20) = 0.0; ARRAY2D(dfdu_,28,21) = 0.0; ARRAY2D(dfdu_,28,22) = 0.0; ARRAY2D(dfdu_,28,23) = 0.0; ARRAY2D(dfdu_,28,24) = 0.0; ARRAY2D(dfdu_,28,25) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,28,26) = 0.0; ARRAY2D(dfdu_,28,27) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,28,28) = -40.0*1.0/tau; ARRAY2D(dfdu_,28,29) = 0.0; ARRAY2D(dfdu_,28,30) = 0.0; ARRAY2D(dfdu_,28,31) = 0.0; ARRAY2D(dfdu_,28,32) = 0.0; ARRAY2D(dfdu_,28,33) = 0.0; ARRAY2D(dfdu_,28,34) = 0.0; ARRAY2D(dfdu_,28,35) = 0.0; ARRAY2D(dfdu_,28,36) = 0.0; ARRAY2D(dfdu_,28,37) = 0.0; ARRAY2D(dfdu_,28,38) = 0.0; ARRAY2D(dfdu_,28,39) = 0.0; ARRAY2D(dfdu_,28,40) = 0.0; ARRAY2D(dfdu_,29,0) = 0.0; ARRAY2D(dfdu_,29,1) = 0.0; ARRAY2D(dfdu_,29,2) = 0.0; ARRAY2D(dfdu_,29,3) = 0.0; ARRAY2D(dfdu_,29,4) = 0.0; ARRAY2D(dfdu_,29,5) = 0.0; ARRAY2D(dfdu_,29,6) = 0.0; ARRAY2D(dfdu_,29,7) = 0.0; ARRAY2D(dfdu_,29,8) = 0.0; ARRAY2D(dfdu_,29,9) = 0.0; ARRAY2D(dfdu_,29,10) = 0.0; ARRAY2D(dfdu_,29,11) = 0.0; ARRAY2D(dfdu_,29,12) = 0.0; ARRAY2D(dfdu_,29,13) = 0.0; ARRAY2D(dfdu_,29,14) = 0.0; ARRAY2D(dfdu_,29,15) = 0.0; ARRAY2D(dfdu_,29,16) = 0.0; ARRAY2D(dfdu_,29,17) = 0.0; ARRAY2D(dfdu_,29,18) = 0.0; ARRAY2D(dfdu_,29,19) = 0.0; ARRAY2D(dfdu_,29,20) = 0.0; ARRAY2D(dfdu_,29,21) = 0.0; ARRAY2D(dfdu_,29,22) = 0.0; ARRAY2D(dfdu_,29,23) = 0.0; ARRAY2D(dfdu_,29,24) = 0.0; ARRAY2D(dfdu_,29,25) = 0.0; ARRAY2D(dfdu_,29,26) = 0.0; ARRAY2D(dfdu_,29,27) = 0.0; ARRAY2D(dfdu_,29,28) = 0.0; ARRAY2D(dfdu_,29,29) = 0.0; ARRAY2D(dfdu_,29,30) = 1.0; ARRAY2D(dfdu_,29,31) = 0.0; ARRAY2D(dfdu_,29,32) = 0.0; ARRAY2D(dfdu_,29,33) = 0.0; ARRAY2D(dfdu_,29,34) = 0.0; ARRAY2D(dfdu_,29,35) = 0.0; ARRAY2D(dfdu_,29,36) = 0.0; ARRAY2D(dfdu_,29,37) = 0.0; ARRAY2D(dfdu_,29,38) = 0.0; ARRAY2D(dfdu_,29,39) = 0.0; ARRAY2D(dfdu_,29,40) = 0.0; ARRAY2D(dfdu_,30,0) = 0.0; ARRAY2D(dfdu_,30,1) = 0.0; ARRAY2D(dfdu_,30,2) = 0.0; ARRAY2D(dfdu_,30,3) = 0.0; ARRAY2D(dfdu_,30,4) = 0.0; ARRAY2D(dfdu_,30,5) = 0.0; ARRAY2D(dfdu_,30,6) = 0.0; ARRAY2D(dfdu_,30,7) = 0.0; ARRAY2D(dfdu_,30,8) = 0.0; ARRAY2D(dfdu_,30,9) = 0.0; ARRAY2D(dfdu_,30,10) = 0.0; ARRAY2D(dfdu_,30,11) = 0.0; ARRAY2D(dfdu_,30,12) = 0.0; ARRAY2D(dfdu_,30,13) = 0.0; ARRAY2D(dfdu_,30,14) = 0.0; ARRAY2D(dfdu_,30,15) = 0.0; ARRAY2D(dfdu_,30,16) = 0.0; ARRAY2D(dfdu_,30,17) = 0.0; ARRAY2D(dfdu_,30,18) = 0.0; ARRAY2D(dfdu_,30,19) = 0.0; ARRAY2D(dfdu_,30,20) = 0.0; ARRAY2D(dfdu_,30,21) = 0.0; ARRAY2D(dfdu_,30,22) = 0.0; ARRAY2D(dfdu_,30,23) = 0.0; ARRAY2D(dfdu_,30,24) = 0.0; ARRAY2D(dfdu_,30,25) = 0.0; ARRAY2D(dfdu_,30,26) = 0.0; ARRAY2D(dfdu_,30,27) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,30,28) = 0.0; ARRAY2D(dfdu_,30,29) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,30,30) = -40.0*1.0/tau; ARRAY2D(dfdu_,30,31) = 0.0; ARRAY2D(dfdu_,30,32) = 0.0; ARRAY2D(dfdu_,30,33) = 0.0; ARRAY2D(dfdu_,30,34) = 0.0; ARRAY2D(dfdu_,30,35) = 0.0; ARRAY2D(dfdu_,30,36) = 0.0; ARRAY2D(dfdu_,30,37) = 0.0; ARRAY2D(dfdu_,30,38) = 0.0; ARRAY2D(dfdu_,30,39) = 0.0; ARRAY2D(dfdu_,30,40) = 0.0; ARRAY2D(dfdu_,31,0) = 0.0; ARRAY2D(dfdu_,31,1) = 0.0; ARRAY2D(dfdu_,31,2) = 0.0; ARRAY2D(dfdu_,31,3) = 0.0; ARRAY2D(dfdu_,31,4) = 0.0; ARRAY2D(dfdu_,31,5) = 0.0; ARRAY2D(dfdu_,31,6) = 0.0; ARRAY2D(dfdu_,31,7) = 0.0; ARRAY2D(dfdu_,31,8) = 0.0; ARRAY2D(dfdu_,31,9) = 0.0; ARRAY2D(dfdu_,31,10) = 0.0; ARRAY2D(dfdu_,31,11) = 0.0; ARRAY2D(dfdu_,31,12) = 0.0; ARRAY2D(dfdu_,31,13) = 0.0; ARRAY2D(dfdu_,31,14) = 0.0; ARRAY2D(dfdu_,31,15) = 0.0; ARRAY2D(dfdu_,31,16) = 0.0; ARRAY2D(dfdu_,31,17) = 0.0; ARRAY2D(dfdu_,31,18) = 0.0; ARRAY2D(dfdu_,31,19) = 0.0; ARRAY2D(dfdu_,31,20) = 0.0; ARRAY2D(dfdu_,31,21) = 0.0; ARRAY2D(dfdu_,31,22) = 0.0; ARRAY2D(dfdu_,31,23) = 0.0; ARRAY2D(dfdu_,31,24) = 0.0; ARRAY2D(dfdu_,31,25) = 0.0; ARRAY2D(dfdu_,31,26) = 0.0; ARRAY2D(dfdu_,31,27) = 0.0; ARRAY2D(dfdu_,31,28) = 0.0; ARRAY2D(dfdu_,31,29) = 0.0; ARRAY2D(dfdu_,31,30) = 0.0; ARRAY2D(dfdu_,31,31) = 0.0; ARRAY2D(dfdu_,31,32) = 1.0; ARRAY2D(dfdu_,31,33) = 0.0; ARRAY2D(dfdu_,31,34) = 0.0; ARRAY2D(dfdu_,31,35) = 0.0; ARRAY2D(dfdu_,31,36) = 0.0; ARRAY2D(dfdu_,31,37) = 0.0; ARRAY2D(dfdu_,31,38) = 0.0; ARRAY2D(dfdu_,31,39) = 0.0; ARRAY2D(dfdu_,31,40) = 0.0; ARRAY2D(dfdu_,32,0) = 0.0; ARRAY2D(dfdu_,32,1) = 0.0; ARRAY2D(dfdu_,32,2) = 0.0; ARRAY2D(dfdu_,32,3) = 0.0; ARRAY2D(dfdu_,32,4) = 0.0; ARRAY2D(dfdu_,32,5) = 0.0; ARRAY2D(dfdu_,32,6) = 0.0; ARRAY2D(dfdu_,32,7) = 0.0; ARRAY2D(dfdu_,32,8) = 0.0; ARRAY2D(dfdu_,32,9) = 0.0; ARRAY2D(dfdu_,32,10) = 0.0; ARRAY2D(dfdu_,32,11) = 0.0; ARRAY2D(dfdu_,32,12) = 0.0; ARRAY2D(dfdu_,32,13) = 0.0; ARRAY2D(dfdu_,32,14) = 0.0; ARRAY2D(dfdu_,32,15) = 0.0; ARRAY2D(dfdu_,32,16) = 0.0; ARRAY2D(dfdu_,32,17) = 0.0; ARRAY2D(dfdu_,32,18) = 0.0; ARRAY2D(dfdu_,32,19) = 0.0; ARRAY2D(dfdu_,32,20) = 0.0; ARRAY2D(dfdu_,32,21) = 0.0; ARRAY2D(dfdu_,32,22) = 0.0; ARRAY2D(dfdu_,32,23) = 0.0; ARRAY2D(dfdu_,32,24) = 0.0; ARRAY2D(dfdu_,32,25) = 0.0; ARRAY2D(dfdu_,32,26) = 0.0; ARRAY2D(dfdu_,32,27) = 0.0; ARRAY2D(dfdu_,32,28) = 0.0; ARRAY2D(dfdu_,32,29) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,32,30) = 0.0; ARRAY2D(dfdu_,32,31) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,32,32) = -40.0*1.0/tau; ARRAY2D(dfdu_,32,33) = 0.0; ARRAY2D(dfdu_,32,34) = 0.0; ARRAY2D(dfdu_,32,35) = 0.0; ARRAY2D(dfdu_,32,36) = 0.0; ARRAY2D(dfdu_,32,37) = 0.0; ARRAY2D(dfdu_,32,38) = 0.0; ARRAY2D(dfdu_,32,39) = 0.0; ARRAY2D(dfdu_,32,40) = 0.0; ARRAY2D(dfdu_,33,0) = 0.0; ARRAY2D(dfdu_,33,1) = 0.0; ARRAY2D(dfdu_,33,2) = 0.0; ARRAY2D(dfdu_,33,3) = 0.0; ARRAY2D(dfdu_,33,4) = 0.0; ARRAY2D(dfdu_,33,5) = 0.0; ARRAY2D(dfdu_,33,6) = 0.0; ARRAY2D(dfdu_,33,7) = 0.0; ARRAY2D(dfdu_,33,8) = 0.0; ARRAY2D(dfdu_,33,9) = 0.0; ARRAY2D(dfdu_,33,10) = 0.0; ARRAY2D(dfdu_,33,11) = 0.0; ARRAY2D(dfdu_,33,12) = 0.0; ARRAY2D(dfdu_,33,13) = 0.0; ARRAY2D(dfdu_,33,14) = 0.0; ARRAY2D(dfdu_,33,15) = 0.0; ARRAY2D(dfdu_,33,16) = 0.0; ARRAY2D(dfdu_,33,17) = 0.0; ARRAY2D(dfdu_,33,18) = 0.0; ARRAY2D(dfdu_,33,19) = 0.0; ARRAY2D(dfdu_,33,20) = 0.0; ARRAY2D(dfdu_,33,21) = 0.0; ARRAY2D(dfdu_,33,22) = 0.0; ARRAY2D(dfdu_,33,23) = 0.0; ARRAY2D(dfdu_,33,24) = 0.0; ARRAY2D(dfdu_,33,25) = 0.0; ARRAY2D(dfdu_,33,26) = 0.0; ARRAY2D(dfdu_,33,27) = 0.0; ARRAY2D(dfdu_,33,28) = 0.0; ARRAY2D(dfdu_,33,29) = 0.0; ARRAY2D(dfdu_,33,30) = 0.0; ARRAY2D(dfdu_,33,31) = 0.0; ARRAY2D(dfdu_,33,32) = 0.0; ARRAY2D(dfdu_,33,33) = 0.0; ARRAY2D(dfdu_,33,34) = 1.0; ARRAY2D(dfdu_,33,35) = 0.0; ARRAY2D(dfdu_,33,36) = 0.0; ARRAY2D(dfdu_,33,37) = 0.0; ARRAY2D(dfdu_,33,38) = 0.0; ARRAY2D(dfdu_,33,39) = 0.0; ARRAY2D(dfdu_,33,40) = 0.0; ARRAY2D(dfdu_,34,0) = 0.0; ARRAY2D(dfdu_,34,1) = 0.0; ARRAY2D(dfdu_,34,2) = 0.0; ARRAY2D(dfdu_,34,3) = 0.0; ARRAY2D(dfdu_,34,4) = 0.0; ARRAY2D(dfdu_,34,5) = 0.0; ARRAY2D(dfdu_,34,6) = 0.0; ARRAY2D(dfdu_,34,7) = 0.0; ARRAY2D(dfdu_,34,8) = 0.0; ARRAY2D(dfdu_,34,9) = 0.0; ARRAY2D(dfdu_,34,10) = 0.0; ARRAY2D(dfdu_,34,11) = 0.0; ARRAY2D(dfdu_,34,12) = 0.0; ARRAY2D(dfdu_,34,13) = 0.0; ARRAY2D(dfdu_,34,14) = 0.0; ARRAY2D(dfdu_,34,15) = 0.0; ARRAY2D(dfdu_,34,16) = 0.0; ARRAY2D(dfdu_,34,17) = 0.0; ARRAY2D(dfdu_,34,18) = 0.0; ARRAY2D(dfdu_,34,19) = 0.0; ARRAY2D(dfdu_,34,20) = 0.0; ARRAY2D(dfdu_,34,21) = 0.0; ARRAY2D(dfdu_,34,22) = 0.0; ARRAY2D(dfdu_,34,23) = 0.0; ARRAY2D(dfdu_,34,24) = 0.0; ARRAY2D(dfdu_,34,25) = 0.0; ARRAY2D(dfdu_,34,26) = 0.0; ARRAY2D(dfdu_,34,27) = 0.0; ARRAY2D(dfdu_,34,28) = 0.0; ARRAY2D(dfdu_,34,29) = 0.0; ARRAY2D(dfdu_,34,30) = 0.0; ARRAY2D(dfdu_,34,31) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,34,32) = 0.0; ARRAY2D(dfdu_,34,33) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,34,34) = -40.0*1.0/tau; ARRAY2D(dfdu_,34,35) = 0.0; ARRAY2D(dfdu_,34,36) = 0.0; ARRAY2D(dfdu_,34,37) = 0.0; ARRAY2D(dfdu_,34,38) = 0.0; ARRAY2D(dfdu_,34,39) = 0.0; ARRAY2D(dfdu_,34,40) = 0.0; ARRAY2D(dfdu_,35,0) = 0.0; ARRAY2D(dfdu_,35,1) = 0.0; ARRAY2D(dfdu_,35,2) = 0.0; ARRAY2D(dfdu_,35,3) = 0.0; ARRAY2D(dfdu_,35,4) = 0.0; ARRAY2D(dfdu_,35,5) = 0.0; ARRAY2D(dfdu_,35,6) = 0.0; ARRAY2D(dfdu_,35,7) = 0.0; ARRAY2D(dfdu_,35,8) = 0.0; ARRAY2D(dfdu_,35,9) = 0.0; ARRAY2D(dfdu_,35,10) = 0.0; ARRAY2D(dfdu_,35,11) = 0.0; ARRAY2D(dfdu_,35,12) = 0.0; ARRAY2D(dfdu_,35,13) = 0.0; ARRAY2D(dfdu_,35,14) = 0.0; ARRAY2D(dfdu_,35,15) = 0.0; ARRAY2D(dfdu_,35,16) = 0.0; ARRAY2D(dfdu_,35,17) = 0.0; ARRAY2D(dfdu_,35,18) = 0.0; ARRAY2D(dfdu_,35,19) = 0.0; ARRAY2D(dfdu_,35,20) = 0.0; ARRAY2D(dfdu_,35,21) = 0.0; ARRAY2D(dfdu_,35,22) = 0.0; ARRAY2D(dfdu_,35,23) = 0.0; ARRAY2D(dfdu_,35,24) = 0.0; ARRAY2D(dfdu_,35,25) = 0.0; ARRAY2D(dfdu_,35,26) = 0.0; ARRAY2D(dfdu_,35,27) = 0.0; ARRAY2D(dfdu_,35,28) = 0.0; ARRAY2D(dfdu_,35,29) = 0.0; ARRAY2D(dfdu_,35,30) = 0.0; ARRAY2D(dfdu_,35,31) = 0.0; ARRAY2D(dfdu_,35,32) = 0.0; ARRAY2D(dfdu_,35,33) = 0.0; ARRAY2D(dfdu_,35,34) = 0.0; ARRAY2D(dfdu_,35,35) = 0.0; ARRAY2D(dfdu_,35,36) = 1.0; ARRAY2D(dfdu_,35,37) = 0.0; ARRAY2D(dfdu_,35,38) = 0.0; ARRAY2D(dfdu_,35,39) = 0.0; ARRAY2D(dfdu_,35,40) = 0.0; ARRAY2D(dfdu_,36,0) = 0.0; ARRAY2D(dfdu_,36,1) = 0.0; ARRAY2D(dfdu_,36,2) = 0.0; ARRAY2D(dfdu_,36,3) = 0.0; ARRAY2D(dfdu_,36,4) = 0.0; ARRAY2D(dfdu_,36,5) = 0.0; ARRAY2D(dfdu_,36,6) = 0.0; ARRAY2D(dfdu_,36,7) = 0.0; ARRAY2D(dfdu_,36,8) = 0.0; ARRAY2D(dfdu_,36,9) = 0.0; ARRAY2D(dfdu_,36,10) = 0.0; ARRAY2D(dfdu_,36,11) = 0.0; ARRAY2D(dfdu_,36,12) = 0.0; ARRAY2D(dfdu_,36,13) = 0.0; ARRAY2D(dfdu_,36,14) = 0.0; ARRAY2D(dfdu_,36,15) = 0.0; ARRAY2D(dfdu_,36,16) = 0.0; ARRAY2D(dfdu_,36,17) = 0.0; ARRAY2D(dfdu_,36,18) = 0.0; ARRAY2D(dfdu_,36,19) = 0.0; ARRAY2D(dfdu_,36,20) = 0.0; ARRAY2D(dfdu_,36,21) = 0.0; ARRAY2D(dfdu_,36,22) = 0.0; ARRAY2D(dfdu_,36,23) = 0.0; ARRAY2D(dfdu_,36,24) = 0.0; ARRAY2D(dfdu_,36,25) = 0.0; ARRAY2D(dfdu_,36,26) = 0.0; ARRAY2D(dfdu_,36,27) = 0.0; ARRAY2D(dfdu_,36,28) = 0.0; ARRAY2D(dfdu_,36,29) = 0.0; ARRAY2D(dfdu_,36,30) = 0.0; ARRAY2D(dfdu_,36,31) = 0.0; ARRAY2D(dfdu_,36,32) = 0.0; ARRAY2D(dfdu_,36,33) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,36,34) = 0.0; ARRAY2D(dfdu_,36,35) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,36,36) = -40.0*1.0/tau; ARRAY2D(dfdu_,36,37) = 0.0; ARRAY2D(dfdu_,36,38) = 0.0; ARRAY2D(dfdu_,36,39) = 0.0; ARRAY2D(dfdu_,36,40) = 0.0; ARRAY2D(dfdu_,37,0) = 0.0; ARRAY2D(dfdu_,37,1) = 0.0; ARRAY2D(dfdu_,37,2) = 0.0; ARRAY2D(dfdu_,37,3) = 0.0; ARRAY2D(dfdu_,37,4) = 0.0; ARRAY2D(dfdu_,37,5) = 0.0; ARRAY2D(dfdu_,37,6) = 0.0; ARRAY2D(dfdu_,37,7) = 0.0; ARRAY2D(dfdu_,37,8) = 0.0; ARRAY2D(dfdu_,37,9) = 0.0; ARRAY2D(dfdu_,37,10) = 0.0; ARRAY2D(dfdu_,37,11) = 0.0; ARRAY2D(dfdu_,37,12) = 0.0; ARRAY2D(dfdu_,37,13) = 0.0; ARRAY2D(dfdu_,37,14) = 0.0; ARRAY2D(dfdu_,37,15) = 0.0; ARRAY2D(dfdu_,37,16) = 0.0; ARRAY2D(dfdu_,37,17) = 0.0; ARRAY2D(dfdu_,37,18) = 0.0; ARRAY2D(dfdu_,37,19) = 0.0; ARRAY2D(dfdu_,37,20) = 0.0; ARRAY2D(dfdu_,37,21) = 0.0; ARRAY2D(dfdu_,37,22) = 0.0; ARRAY2D(dfdu_,37,23) = 0.0; ARRAY2D(dfdu_,37,24) = 0.0; ARRAY2D(dfdu_,37,25) = 0.0; ARRAY2D(dfdu_,37,26) = 0.0; ARRAY2D(dfdu_,37,27) = 0.0; ARRAY2D(dfdu_,37,28) = 0.0; ARRAY2D(dfdu_,37,29) = 0.0; ARRAY2D(dfdu_,37,30) = 0.0; ARRAY2D(dfdu_,37,31) = 0.0; ARRAY2D(dfdu_,37,32) = 0.0; ARRAY2D(dfdu_,37,33) = 0.0; ARRAY2D(dfdu_,37,34) = 0.0; ARRAY2D(dfdu_,37,35) = 0.0; ARRAY2D(dfdu_,37,36) = 0.0; ARRAY2D(dfdu_,37,37) = 0.0; ARRAY2D(dfdu_,37,38) = 1.0; ARRAY2D(dfdu_,37,39) = 0.0; ARRAY2D(dfdu_,37,40) = 0.0; ARRAY2D(dfdu_,38,0) = 0.0; ARRAY2D(dfdu_,38,1) = 0.0; ARRAY2D(dfdu_,38,2) = 0.0; ARRAY2D(dfdu_,38,3) = 0.0; ARRAY2D(dfdu_,38,4) = 0.0; ARRAY2D(dfdu_,38,5) = 0.0; ARRAY2D(dfdu_,38,6) = 0.0; ARRAY2D(dfdu_,38,7) = 0.0; ARRAY2D(dfdu_,38,8) = 0.0; ARRAY2D(dfdu_,38,9) = 0.0; ARRAY2D(dfdu_,38,10) = 0.0; ARRAY2D(dfdu_,38,11) = 0.0; ARRAY2D(dfdu_,38,12) = 0.0; ARRAY2D(dfdu_,38,13) = 0.0; ARRAY2D(dfdu_,38,14) = 0.0; ARRAY2D(dfdu_,38,15) = 0.0; ARRAY2D(dfdu_,38,16) = 0.0; ARRAY2D(dfdu_,38,17) = 0.0; ARRAY2D(dfdu_,38,18) = 0.0; ARRAY2D(dfdu_,38,19) = 0.0; ARRAY2D(dfdu_,38,20) = 0.0; ARRAY2D(dfdu_,38,21) = 0.0; ARRAY2D(dfdu_,38,22) = 0.0; ARRAY2D(dfdu_,38,23) = 0.0; ARRAY2D(dfdu_,38,24) = 0.0; ARRAY2D(dfdu_,38,25) = 0.0; ARRAY2D(dfdu_,38,26) = 0.0; ARRAY2D(dfdu_,38,27) = 0.0; ARRAY2D(dfdu_,38,28) = 0.0; ARRAY2D(dfdu_,38,29) = 0.0; ARRAY2D(dfdu_,38,30) = 0.0; ARRAY2D(dfdu_,38,31) = 0.0; ARRAY2D(dfdu_,38,32) = 0.0; ARRAY2D(dfdu_,38,33) = 0.0; ARRAY2D(dfdu_,38,34) = 0.0; ARRAY2D(dfdu_,38,35) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,38,36) = 0.0; ARRAY2D(dfdu_,38,37) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,38,38) = -40.0*1.0/tau; ARRAY2D(dfdu_,38,39) = 0.0; ARRAY2D(dfdu_,38,40) = 0.0; ARRAY2D(dfdu_,39,0) = 0.0; ARRAY2D(dfdu_,39,1) = 0.0; ARRAY2D(dfdu_,39,2) = 0.0; ARRAY2D(dfdu_,39,3) = 0.0; ARRAY2D(dfdu_,39,4) = 0.0; ARRAY2D(dfdu_,39,5) = 0.0; ARRAY2D(dfdu_,39,6) = 0.0; ARRAY2D(dfdu_,39,7) = 0.0; ARRAY2D(dfdu_,39,8) = 0.0; ARRAY2D(dfdu_,39,9) = 0.0; ARRAY2D(dfdu_,39,10) = 0.0; ARRAY2D(dfdu_,39,11) = 0.0; ARRAY2D(dfdu_,39,12) = 0.0; ARRAY2D(dfdu_,39,13) = 0.0; ARRAY2D(dfdu_,39,14) = 0.0; ARRAY2D(dfdu_,39,15) = 0.0; ARRAY2D(dfdu_,39,16) = 0.0; ARRAY2D(dfdu_,39,17) = 0.0; ARRAY2D(dfdu_,39,18) = 0.0; ARRAY2D(dfdu_,39,19) = 0.0; ARRAY2D(dfdu_,39,20) = 0.0; ARRAY2D(dfdu_,39,21) = 0.0; ARRAY2D(dfdu_,39,22) = 0.0; ARRAY2D(dfdu_,39,23) = 0.0; ARRAY2D(dfdu_,39,24) = 0.0; ARRAY2D(dfdu_,39,25) = 0.0; ARRAY2D(dfdu_,39,26) = 0.0; ARRAY2D(dfdu_,39,27) = 0.0; ARRAY2D(dfdu_,39,28) = 0.0; ARRAY2D(dfdu_,39,29) = 0.0; ARRAY2D(dfdu_,39,30) = 0.0; ARRAY2D(dfdu_,39,31) = 0.0; ARRAY2D(dfdu_,39,32) = 0.0; ARRAY2D(dfdu_,39,33) = 0.0; ARRAY2D(dfdu_,39,34) = 0.0; ARRAY2D(dfdu_,39,35) = 0.0; ARRAY2D(dfdu_,39,36) = 0.0; ARRAY2D(dfdu_,39,37) = 0.0; ARRAY2D(dfdu_,39,38) = 0.0; ARRAY2D(dfdu_,39,39) = 0.0; ARRAY2D(dfdu_,39,40) = 1.0; ARRAY2D(dfdu_,40,0) = 0.0; ARRAY2D(dfdu_,40,1) = 0.0; ARRAY2D(dfdu_,40,2) = 0.0; ARRAY2D(dfdu_,40,3) = 0.0; ARRAY2D(dfdu_,40,4) = 0.0; ARRAY2D(dfdu_,40,5) = 0.0; ARRAY2D(dfdu_,40,6) = 0.0; ARRAY2D(dfdu_,40,7) = 0.0; ARRAY2D(dfdu_,40,8) = 0.0; ARRAY2D(dfdu_,40,9) = 0.0; ARRAY2D(dfdu_,40,10) = 0.0; ARRAY2D(dfdu_,40,11) = 0.0; ARRAY2D(dfdu_,40,12) = 0.0; ARRAY2D(dfdu_,40,13) = 0.0; ARRAY2D(dfdu_,40,14) = 0.0; ARRAY2D(dfdu_,40,15) = 0.0; ARRAY2D(dfdu_,40,16) = 0.0; ARRAY2D(dfdu_,40,17) = 0.0; ARRAY2D(dfdu_,40,18) = 0.0; ARRAY2D(dfdu_,40,19) = 0.0; ARRAY2D(dfdu_,40,20) = 0.0; ARRAY2D(dfdu_,40,21) = 0.0; ARRAY2D(dfdu_,40,22) = 0.0; ARRAY2D(dfdu_,40,23) = 0.0; ARRAY2D(dfdu_,40,24) = 0.0; ARRAY2D(dfdu_,40,25) = 0.0; ARRAY2D(dfdu_,40,26) = 0.0; ARRAY2D(dfdu_,40,27) = 0.0; ARRAY2D(dfdu_,40,28) = 0.0; ARRAY2D(dfdu_,40,29) = 0.0; ARRAY2D(dfdu_,40,30) = 0.0; ARRAY2D(dfdu_,40,31) = 0.0; ARRAY2D(dfdu_,40,32) = 0.0; ARRAY2D(dfdu_,40,33) = 0.0; ARRAY2D(dfdu_,40,34) = 0.0; ARRAY2D(dfdu_,40,35) = 0.0; ARRAY2D(dfdu_,40,36) = 0.0; ARRAY2D(dfdu_,40,37) = 800.0*1.0/(tau*tau); ARRAY2D(dfdu_,40,38) = 0.0; ARRAY2D(dfdu_,40,39) = -800.0*1.0/(tau*tau); ARRAY2D(dfdu_,40,40) = -40.0*1.0/tau; if (ijac_ == 1) return 0; /* * The Jacobian with respect to the parameters. */ ARRAY2D(dfdp_,0,0) = symbol9_1_20/( ((symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20))*(symbol9_1_20*(symbol9_1_20*symbol9_1_20)*(symbol9_1_20*symbol9_1_20)))+1.0); ARRAY2D(dfdp_,0,1) = -x; ARRAY2D(dfdp_,0,2) = 0.0; ARRAY2D(dfdp_,1,0) = 0.0; ARRAY2D(dfdp_,1,1) = 0.0; ARRAY2D(dfdp_,1,2) = 0.0; ARRAY2D(dfdp_,2,0) = 0.0; ARRAY2D(dfdp_,2,1) = 0.0; ARRAY2D(dfdp_,2,2) = 800.0*( symbol9_1_1-x)/(tau*tau*tau)+40.0*1.0/(tau*tau)*( symbol9_2_1+20.0*( symbol9_1_1-x)/tau); ARRAY2D(dfdp_,3,0) = 0.0; ARRAY2D(dfdp_,3,1) = 0.0; ARRAY2D(dfdp_,3,2) = 0.0; ARRAY2D(dfdp_,4,0) = 0.0; ARRAY2D(dfdp_,4,1) = 0.0; ARRAY2D(dfdp_,4,2) = -800.0*1.0/(tau*tau*tau)*( symbol9_1_1-symbol9_1_2)+40.0*( symbol9_2_2+-20.0*1.0/tau*( symbol9_1_1-symbol9_1_2))/(tau*tau); ARRAY2D(dfdp_,5,0) = 0.0; ARRAY2D(dfdp_,5,1) = 0.0; ARRAY2D(dfdp_,5,2) = 0.0; ARRAY2D(dfdp_,6,0) = 0.0; ARRAY2D(dfdp_,6,1) = 0.0; ARRAY2D(dfdp_,6,2) = 800.0*( symbol9_1_3-symbol9_1_2)/(tau*tau*tau)+40.0*1.0/(tau*tau)*( 20.0*( symbol9_1_3-symbol9_1_2)/tau+symbol9_2_3); ARRAY2D(dfdp_,7,0) = 0.0; ARRAY2D(dfdp_,7,1) = 0.0; ARRAY2D(dfdp_,7,2) = 0.0; ARRAY2D(dfdp_,8,0) = 0.0; ARRAY2D(dfdp_,8,1) = 0.0; ARRAY2D(dfdp_,8,2) = -800.0*1.0/(tau*tau*tau)*( symbol9_1_3-symbol9_1_4)+40.0*1.0/(tau*tau)*( symbol9_2_4+-20.0*1.0/tau*( symbol9_1_3-symbol9_1_4)); ARRAY2D(dfdp_,9,0) = 0.0; ARRAY2D(dfdp_,9,1) = 0.0; ARRAY2D(dfdp_,9,2) = 0.0; ARRAY2D(dfdp_,10,0) = 0.0; ARRAY2D(dfdp_,10,1) = 0.0; ARRAY2D(dfdp_,10,2) = 800.0*1.0/(tau*tau*tau)*( symbol9_1_5-symbol9_1_4)+40.0*( 20.0*1.0/tau*( symbol9_1_5-symbol9_1_4)+symbol9_2_5)/(tau*tau); ARRAY2D(dfdp_,11,0) = 0.0; ARRAY2D(dfdp_,11,1) = 0.0; ARRAY2D(dfdp_,11,2) = 0.0; ARRAY2D(dfdp_,12,0) = 0.0; ARRAY2D(dfdp_,12,1) = 0.0; ARRAY2D(dfdp_,12,2) = -40.0*1.0/(tau*tau)*( 20.0*( symbol9_1_5-symbol9_1_6)/tau-symbol9_2_6)+-800.0*( symbol9_1_5-symbol9_1_6)/(tau*tau*tau); ARRAY2D(dfdp_,13,0) = 0.0; ARRAY2D(dfdp_,13,1) = 0.0; ARRAY2D(dfdp_,13,2) = 0.0; ARRAY2D(dfdp_,14,0) = 0.0; ARRAY2D(dfdp_,14,1) = 0.0; ARRAY2D(dfdp_,14,2) = -40.0*1.0/(tau*tau)*( 20.0*1.0/tau*( symbol9_1_6-symbol9_1_7)-symbol9_2_7)+-800.0*1.0/(tau*tau*tau)*( symbol9_1_6-symbol9_1_7); ARRAY2D(dfdp_,15,0) = 0.0; ARRAY2D(dfdp_,15,1) = 0.0; ARRAY2D(dfdp_,15,2) = 0.0; ARRAY2D(dfdp_,16,0) = 0.0; ARRAY2D(dfdp_,16,1) = 0.0; ARRAY2D(dfdp_,16,2) = 800.0*1.0/(tau*tau*tau)*( symbol9_1_8-symbol9_1_7)+40.0*1.0/(tau*tau)*( symbol9_2_8+20.0*1.0/tau*( symbol9_1_8-symbol9_1_7)); ARRAY2D(dfdp_,17,0) = 0.0; ARRAY2D(dfdp_,17,1) = 0.0; ARRAY2D(dfdp_,17,2) = 0.0; ARRAY2D(dfdp_,18,0) = 0.0; ARRAY2D(dfdp_,18,1) = 0.0; ARRAY2D(dfdp_,18,2) = -800.0*( symbol9_1_8-symbol9_1_9)/(tau*tau*tau)+40.0*1.0/(tau*tau)*( symbol9_2_9+-20.0*( symbol9_1_8-symbol9_1_9)/tau); ARRAY2D(dfdp_,19,0) = 0.0; ARRAY2D(dfdp_,19,1) = 0.0; ARRAY2D(dfdp_,19,2) = 0.0; ARRAY2D(dfdp_,20,0) = 0.0; ARRAY2D(dfdp_,20,1) = 0.0; ARRAY2D(dfdp_,20,2) = -40.0*( 20.0*( symbol9_1_9-symbol9_1_10)/tau-symbol9_2_10)/(tau*tau)+-800.0*( symbol9_1_9-symbol9_1_10)/(tau*tau*tau); ARRAY2D(dfdp_,21,0) = 0.0; ARRAY2D(dfdp_,21,1) = 0.0; ARRAY2D(dfdp_,21,2) = 0.0; ARRAY2D(dfdp_,22,0) = 0.0; ARRAY2D(dfdp_,22,1) = 0.0; ARRAY2D(dfdp_,22,2) = 40.0*( 20.0*( symbol9_1_11-symbol9_1_10)/tau+symbol9_2_11)/(tau*tau)+800.0*( symbol9_1_11-symbol9_1_10)/(tau*tau*tau); ARRAY2D(dfdp_,23,0) = 0.0; ARRAY2D(dfdp_,23,1) = 0.0; ARRAY2D(dfdp_,23,2) = 0.0; ARRAY2D(dfdp_,24,0) = 0.0; ARRAY2D(dfdp_,24,1) = 0.0; ARRAY2D(dfdp_,24,2) = -40.0*( 20.0*1.0/tau*( symbol9_1_11-symbol9_1_12)-symbol9_2_12)/(tau*tau)+-800.0*1.0/(tau*tau*tau)*( symbol9_1_11-symbol9_1_12); ARRAY2D(dfdp_,25,0) = 0.0; ARRAY2D(dfdp_,25,1) = 0.0; ARRAY2D(dfdp_,25,2) = 0.0; ARRAY2D(dfdp_,26,0) = 0.0; ARRAY2D(dfdp_,26,1) = 0.0; ARRAY2D(dfdp_,26,2) = 40.0*1.0/(tau*tau)*( 20.0*1.0/tau*( symbol9_1_13-symbol9_1_12)+symbol9_2_13)+800.0*1.0/(tau*tau*tau)*( symbol9_1_13-symbol9_1_12); ARRAY2D(dfdp_,27,0) = 0.0; ARRAY2D(dfdp_,27,1) = 0.0; ARRAY2D(dfdp_,27,2) = 0.0; ARRAY2D(dfdp_,28,0) = 0.0; ARRAY2D(dfdp_,28,1) = 0.0; ARRAY2D(dfdp_,28,2) = 800.0*1.0/(tau*tau*tau)*( symbol9_1_14-symbol9_1_13)+40.0*( symbol9_2_14+20.0*1.0/tau*( symbol9_1_14-symbol9_1_13))/(tau*tau); ARRAY2D(dfdp_,29,0) = 0.0; ARRAY2D(dfdp_,29,1) = 0.0; ARRAY2D(dfdp_,29,2) = 0.0; ARRAY2D(dfdp_,30,0) = 0.0; ARRAY2D(dfdp_,30,1) = 0.0; ARRAY2D(dfdp_,30,2) = 40.0*1.0/(tau*tau)*( symbol9_2_15+-20.0*( symbol9_1_14-symbol9_1_15)/tau)+-800.0*( symbol9_1_14-symbol9_1_15)/(tau*tau*tau); ARRAY2D(dfdp_,31,0) = 0.0; ARRAY2D(dfdp_,31,1) = 0.0; ARRAY2D(dfdp_,31,2) = 0.0; ARRAY2D(dfdp_,32,0) = 0.0; ARRAY2D(dfdp_,32,1) = 0.0; ARRAY2D(dfdp_,32,2) = -800.0*1.0/(tau*tau*tau)*( symbol9_1_15-symbol9_1_16)+40.0*( symbol9_2_16+-20.0*1.0/tau*( symbol9_1_15-symbol9_1_16))/(tau*tau); ARRAY2D(dfdp_,33,0) = 0.0; ARRAY2D(dfdp_,33,1) = 0.0; ARRAY2D(dfdp_,33,2) = 0.0; ARRAY2D(dfdp_,34,0) = 0.0; ARRAY2D(dfdp_,34,1) = 0.0; ARRAY2D(dfdp_,34,2) = 40.0*( 20.0*1.0/tau*( symbol9_1_17-symbol9_1_16)+symbol9_2_17)/(tau*tau)+800.0*1.0/(tau*tau*tau)*( symbol9_1_17-symbol9_1_16); ARRAY2D(dfdp_,35,0) = 0.0; ARRAY2D(dfdp_,35,1) = 0.0; ARRAY2D(dfdp_,35,2) = 0.0; ARRAY2D(dfdp_,36,0) = 0.0; ARRAY2D(dfdp_,36,1) = 0.0; ARRAY2D(dfdp_,36,2) = 40.0*( symbol9_2_18+-20.0*1.0/tau*( symbol9_1_17-symbol9_1_18))/(tau*tau)+-800.0*1.0/(tau*tau*tau)*( symbol9_1_17-symbol9_1_18); ARRAY2D(dfdp_,37,0) = 0.0; ARRAY2D(dfdp_,37,1) = 0.0; ARRAY2D(dfdp_,37,2) = 0.0; ARRAY2D(dfdp_,38,0) = 0.0; ARRAY2D(dfdp_,38,1) = 0.0; ARRAY2D(dfdp_,38,2) = 800.0*( symbol9_1_19-symbol9_1_18)/(tau*tau*tau)+40.0*1.0/(tau*tau)*( 20.0*( symbol9_1_19-symbol9_1_18)/tau+symbol9_2_19); ARRAY2D(dfdp_,39,0) = 0.0; ARRAY2D(dfdp_,39,1) = 0.0; ARRAY2D(dfdp_,39,2) = 0.0; ARRAY2D(dfdp_,40,0) = 0.0; ARRAY2D(dfdp_,40,1) = 0.0; ARRAY2D(dfdp_,40,2) = -800.0*1.0/(tau*tau*tau)*( symbol9_1_19-symbol9_1_20)+40.0*1.0/(tau*tau)*( symbol9_2_20+-20.0*1.0/tau*( symbol9_1_19-symbol9_1_20)); return 0; } /* * STPNT Gives a starting point */ int stpnt(integer ndim_, doublereal t_, doublereal *u_, doublereal *par_) { double x, symbol9_1_1, symbol9_2_1, symbol9_1_2, symbol9_2_2, symbol9_1_3, symbol9_2_3, symbol9_1_4, symbol9_2_4, symbol9_1_5, symbol9_2_5, symbol9_1_6, symbol9_2_6, symbol9_1_7, symbol9_2_7, symbol9_1_8, symbol9_2_8, symbol9_1_9, symbol9_2_9, symbol9_1_10, symbol9_2_10, symbol9_1_11, symbol9_2_11, symbol9_1_12, symbol9_2_12, symbol9_1_13, symbol9_2_13, symbol9_1_14, symbol9_2_14, symbol9_1_15, symbol9_2_15, symbol9_1_16, symbol9_2_16, symbol9_1_17, symbol9_2_17, symbol9_1_18, symbol9_2_18, symbol9_1_19, symbol9_2_19, symbol9_1_20, symbol9_2_20; double a, b, tau; /* Change the parameter values and the starting point to correct values! */ a = 2.0000000000000001e-01; b = 1.0000000000000001e-01; tau = 1.0000000000000001e-01; x = 1.0; symbol9_1_1 = 1.0; symbol9_2_1 = 0.0; symbol9_1_2 = 1.0; symbol9_2_2 = 0.0; symbol9_1_3 = 1.0; symbol9_2_3 = 0.0; symbol9_1_4 = 1.0; symbol9_2_4 = 0.0; symbol9_1_5 = 1.0; symbol9_2_5 = 0.0; symbol9_1_6 = 1.0; symbol9_2_6 = 0.0; symbol9_1_7 = 1.0; symbol9_2_7 = 0.0; symbol9_1_8 = 1.0; symbol9_2_8 = 0.0; symbol9_1_9 = 1.0; symbol9_2_9 = 0.0; symbol9_1_10 = 1.0; symbol9_2_10 = 0.0; symbol9_1_11 = 1.0; symbol9_2_11 = 0.0; symbol9_1_12 = 1.0; symbol9_2_12 = 0.0; symbol9_1_13 = 1.0; symbol9_2_13 = 0.0; symbol9_1_14 = 1.0; symbol9_2_14 = 0.0; symbol9_1_15 = 1.0; symbol9_2_15 = 0.0; symbol9_1_16 = 1.0; symbol9_2_16 = 0.0; symbol9_1_17 = 1.0; symbol9_2_17 = 0.0; symbol9_1_18 = 1.0; symbol9_2_18 = 0.0; symbol9_1_19 = 1.0; symbol9_2_19 = 0.0; symbol9_1_20 = 1.0; symbol9_2_20 = 0.0; par_[0] = a; par_[1] = b; par_[2] = tau; u_[0] = x; u_[1] = symbol9_1_1; u_[2] = symbol9_2_1; u_[3] = symbol9_1_2; u_[4] = symbol9_2_2; u_[5] = symbol9_1_3; u_[6] = symbol9_2_3; u_[7] = symbol9_1_4; u_[8] = symbol9_2_4; u_[9] = symbol9_1_5; u_[10] = symbol9_2_5; u_[11] = symbol9_1_6; u_[12] = symbol9_2_6; u_[13] = symbol9_1_7; u_[14] = symbol9_2_7; u_[15] = symbol9_1_8; u_[16] = symbol9_2_8; u_[17] = symbol9_1_9; u_[18] = symbol9_2_9; u_[19] = symbol9_1_10; u_[20] = symbol9_2_10; u_[21] = symbol9_1_11; u_[22] = symbol9_2_11; u_[23] = symbol9_1_12; u_[24] = symbol9_2_12; u_[25] = symbol9_1_13; u_[26] = symbol9_2_13; u_[27] = symbol9_1_14; u_[28] = symbol9_2_14; u_[29] = symbol9_1_15; u_[30] = symbol9_2_15; u_[31] = symbol9_1_16; u_[32] = symbol9_2_16; u_[33] = symbol9_1_17; u_[34] = symbol9_2_17; u_[35] = symbol9_1_18; u_[36] = symbol9_2_18; u_[37] = symbol9_1_19; u_[38] = symbol9_2_19; u_[39] = symbol9_1_20; u_[40] = symbol9_2_20; return 0; } /* * The remaining functions are just stubs. * You will have to edit these by hand if you need these functions. */ /* * BCND Defines the boundary conditions */ int bcnd(integer ndim_, const doublereal *par_, const integer *icp_, integer nbc_, const doublereal *u0_, const doublereal *u1_, integer ijac_, doublereal *fb_, doublereal *dbc_) { return 0; } /* * ICND Defines the integral conditions */ int icnd(integer ndim_, const doublereal *par_, const integer *icp_, integer nint_, const doublereal *u_, const doublereal *uold_, const doublereal *udot_, const doublereal *upold_, integer ijac_, doublereal *fi_, doublereal *dint_) { return 0; } /* * FOPT */ int fopt(integer ndim_, const doublereal *u_, const integer *icp_, const doublereal *par_, integer ijac_, doublereal *fs_, doublereal *dfdu_, doublereal *dfdp_) { return 0; } /* * PVLS */ int pvls(integer ndim_, const doublereal *u_, doublereal *par_) { return 0; }