dxdt[0]=-((a + del)*x[0]) + 16*b*x[16];
dxdt[1]=-(a*(-1 + eps)*x[1]) - 16.*(b + 0.00625*del + 0.0625*a*eps)*x[1] - 2*a*(-1 + eps)*(x[0] + x[2] + x[3] + x[4] + x[5] + x[6] + x[7] + x[8] + x[9] + x[10] + x[11] + x[12] + x[13] + x[14] + x[15] + x[16]) - 2*a*(-1 + eps)*(x[33] + x[34] + x[35] + x[36] + x[37] + x[38] + x[39] + x[40] + x[41] + x[42] + x[43] + x[44] + x[45] + x[46] + x[47] + x[48] + x[50]);
dxdt[2]=16.*b*x[1] - 1.*a*x[2] - 16.*b*x[2] - 0.1*del*x[2];
dxdt[3]=16.*b*x[2] - 1.*a*x[3] - 16.*b*x[3] - 0.1*del*x[3];
dxdt[4]=16.*b*x[3] - 1.*a*x[4] - 16.*b*x[4] - 0.1*del*x[4];
dxdt[5]=16.*b*x[4] - 1.*a*x[5] - 16.*b*x[5] - 0.1*del*x[5];
dxdt[6]=16.*b*x[5] - 1.*a*x[6] - 16.*b*x[6] - 0.1*del*x[6];
dxdt[7]=16.*b*x[6] - 1.*a*x[7] - 16.*b*x[7] - 0.1*del*x[7];
dxdt[8]=16.*b*x[7] - 1.*a*x[8] - 16.*b*x[8] - 0.1*del*x[8];
dxdt[9]=16.*b*x[8] - 1.*a*x[9] - 16.*b*x[9] - 0.1*del*x[9];
dxdt[10]=16.*b*x[9] - 1.*a*x[10] - 16.*b*x[10] - 0.1*del*x[10];
dxdt[11]=16.*b*x[10] - 1.*a*x[11] - 16.*b*x[11] - 0.1*del*x[11];
dxdt[12]=16.*b*x[11] - 1.*a*x[12] - 16.*b*x[12] - 0.1*del*x[12];
dxdt[13]=16.*b*x[12] - 1.*a*x[13] - 16.*b*x[13] - 0.1*del*x[13];
dxdt[14]=16.*b*x[13] - 1.*a*x[14] - 16.*b*x[14] - 0.1*del*x[14];
dxdt[15]=16.*b*x[14] - 1.*a*x[15] - 16.*b*x[15] - 0.1*del*x[15];
dxdt[16]=16.*b*x[15] - 1.*a*x[16] - 16.*b*x[16] - 0.1*del*x[16];
dxdt[17]=N0*sig + (-1.*a - 16.*b - 0.1*del)*x[17] + 2.*a*eps*(x[0] + x[1] + x[2] + x[3] + x[4] + x[5] + x[6] + x[7] + x[8] + x[9] + x[10] + x[11] + x[12] + x[13] + x[14] + x[15] + x[16] + 1.*x[17] + x[18] + x[19] + x[20] + x[21] + x[22] + x[23] + x[24] + x[25] + x[26] + x[27] + x[28] + x[29] + x[30] + x[31] + x[32] + x[33] + x[34] + x[35] + x[36] + x[37] + x[38] + x[39] + x[40] + x[41] + x[42] + x[43] + x[44] + x[45] + x[46] + x[47] + x[48] + x[49] + x[50]);
dxdt[18]=16.*b*x[17] - 1.*a*x[18] - 16.*b*x[18] - 0.1*del*x[18];
dxdt[19]=16.*b*x[18] - 1.*a*x[19] - 16.*b*x[19] - 0.1*del*x[19];
dxdt[20]=16.*b*x[19] - 1.*a*x[20] - 16.*b*x[20] - 0.1*del*x[20];
dxdt[21]=16.*b*x[20] - 1.*a*x[21] - 16.*b*x[21] - 0.1*del*x[21];
dxdt[22]=16.*b*x[21] - 1.*a*x[22] - 16.*b*x[22] - 0.1*del*x[22];
dxdt[23]=16.*b*x[22] - 1.*a*x[23] - 16.*b*x[23] - 0.1*del*x[23];
dxdt[24]=16.*b*x[23] - 1.*a*x[24] - 16.*b*x[24] - 0.1*del*x[24];
dxdt[25]=16.*b*x[24] - 1.*a*x[25] - 16.*b*x[25] - 0.1*del*x[25];
dxdt[26]=16.*b*x[25] - 1.*a*x[26] - 16.*b*x[26] - 0.1*del*x[26];
dxdt[27]=16.*b*x[26] - 1.*a*x[27] - 16.*b*x[27] - 0.1*del*x[27];
dxdt[28]=16.*b*x[27] - 1.*a*x[28] - 16.*b*x[28] - 0.1*del*x[28];
dxdt[29]=16.*b*x[28] - 1.*a*x[29] - 16.*b*x[29] - 0.1*del*x[29];
dxdt[30]=16.*b*x[29] - 1.*a*x[30] - 16.*b*x[30] - 0.1*del*x[30];
dxdt[31]=16.*b*x[30] - 1.*a*x[31] - 16.*b*x[31] - 0.1*del*x[31];
dxdt[32]=16.*b*x[31] - 1.*a*x[32] - 16.*b*x[32] - 0.1*del*x[32];
dxdt[33]=-1.*(a + 16.*b + 0.1*del)*x[33] - 2*a*(-1 + eps)*(x[17] + x[18] + x[19] + x[20] + x[21] + x[22] + x[23] + x[24] + x[25] + x[26] + x[27] + x[28] + x[29] + x[30] + x[31] + x[32] + x[49]);
dxdt[34]=16.*b*x[33] - 1.*a*x[34] - 16.*b*x[34] - 0.1*del*x[34];
dxdt[35]=16.*b*x[34] - 1.*a*x[35] - 16.*b*x[35] - 0.1*del*x[35];
dxdt[36]=16.*b*x[35] - 1.*a*x[36] - 16.*b*x[36] - 0.1*del*x[36];
dxdt[37]=16.*b*x[36] - 1.*a*x[37] - 16.*b*x[37] - 0.1*del*x[37];
dxdt[38]=16.*b*x[37] - 1.*a*x[38] - 16.*b*x[38] - 0.1*del*x[38];
dxdt[39]=16.*b*x[38] - 1.*a*x[39] - 16.*b*x[39] - 0.1*del*x[39];
dxdt[40]=16.*b*x[39] - 1.*a*x[40] - 16.*b*x[40] - 0.1*del*x[40];
dxdt[41]=16.*b*x[40] - 1.*a*x[41] - 16.*b*x[41] - 0.1*del*x[41];
dxdt[42]=16.*b*x[41] - 1.*a*x[42] - 16.*b*x[42] - 0.1*del*x[42];
dxdt[43]=16.*b*x[42] - 1.*a*x[43] - 16.*b*x[43] - 0.1*del*x[43];
dxdt[44]=16.*b*x[43] - 1.*a*x[44] - 16.*b*x[44] - 0.1*del*x[44];
dxdt[45]=16.*b*x[44] - 1.*a*x[45] - 16.*b*x[45] - 0.1*del*x[45];
dxdt[46]=16.*b*x[45] - 1.*a*x[46] - 16.*b*x[46] - 0.1*del*x[46];
dxdt[47]=16.*b*x[46] - 1.*a*x[47] - 16.*b*x[47] - 0.1*del*x[47];
dxdt[48]=16.*b*x[47] - 1.*a*x[48] - 16.*b*x[48] - 0.1*del*x[48];
dxdt[49]=16*b*x[32] - (a + del)*x[49];
dxdt[50]=16*b*x[48] - (a + del)*x[50];