- #1

- 23

- 0

I hope that any body could help me

For [a = 0, a <= 30, a = a + 10,

For[z = 3.7, z <= 5, z = z + 0.2,

For[w = 1000, w <= 20000, w = w + 1000,

{m = NDSolve[{Derivative[1][y1][t] == y2[t],

Derivative[1][y2][t] == y3[t], Derivative[1][y3][t] == y4[t],

Derivative[1][y4][t] == -z^4 y1[t] + 2 z^2 y3[t] +

z^2 w (-t + t^2) y5[t] + a y3[t],

Derivative[1][y5][t] == y6[t],

Derivative[1][y6][t] == (-1 + 2 t) y1[t] + z^2 y5[t] -

a y8[t], Derivative[1][y8][t] == y7[t],

Derivative[1][y7][t] == -y6[t] + z^2 y8[t],

y1[0] == y2[0] == y5[0] == y7[0] == 0, y3[0] == 1,

y4[0] == y6[0] == y8[0] == 0}, {y1, y2, y3, y4, y5, y6, y7,

y8}, {t, 0, 1}]}

{n = NDSolve[{Derivative[1][y9][t] == y10[t],

Derivative[1][y10][t] == y11[t],

Derivative[1][y11][t] == y12[t],

Derivative[1][y12][t] == -z^4 y9[t] + 2 z^2 y11[t] +

z^2 w (-t + t^2) y13[t] + a y11[t],

Derivative[1][y13][t] == y14[t],

Derivative[1][y14][t] == (-1 + 2 t) y9[t] + z^2 y13[t] -

a y16[t], Derivative[1][y16][t] == y15[t],

Derivative[1][y15][t] == -y14[t] + z^2 y16[t],

y9[0] == y10[0] == y13[0] == y15[0] == 0, y12[0] == 1,

y11[0] == y14[0] == y16[0] == 0}, {y9, y10, y11, y12, y13,

y14, y15, y16}, {t, 0, 1}]}

{l = NDSolve[{Derivative[1][y17][t] == y18[t],

Derivative[1][y18][t] == y19[t],

Derivative[1][y19][t] == y20[t],

Derivative[1][y20][t] == -z^4 y17[t] + 2 z^2 y19[t] +

z^2 w (-t + t^2) y21[t] + a y19[t],

Derivative[1][y21][t] == y22[t],

Derivative[1][y22][t] == (-1 + 2 t) y17[t] + z^2 y20[t] -

a y24[t], Derivative[1][y24][t] == y23[t],

Derivative[1][y23][t] == -y22[t] + z^2 y24[t],

y17[0] == y18[0] == y21[0] == y23[0] == 0, y22[0] == 1,

y20[0] == y19[0] == y24[0] == 0}, {y17, y18, y19, y20, y21,

y22, y23, y24}, {t, 0, 1}]}

{c = NDSolve[{Derivative[1][y25][t] == y26[t],

Derivative[1][y26][t] == y27[t],

Derivative[1][y27][t] == y28[t],

Derivative[1][y28][t] == -z^4 y25[t] + 2 z^2 y27[t] +

z^2 w (-t + t^2) y29[t] + a y27[t],

Derivative[1][y29][t] == y30[t],

Derivative[1][y30][t] == (-1 + 2 t) y25[t] + z^2 y28[t] -

a y32[t], Derivative[1][y32][t] == y31[t],

Derivative[1][y31][t] == -y30[t] + z^2 y32[t],

y25[0] == y26[0] == y29[0] == y31[0] == 0, y30[0] == 1,

y28[0] == y27[0] == y32[0] == 0}, {y25, y26, y27, y28, y29,

y30, y31, y32}, {t, 0, 1}]}

{d = Array[f, {4, 4}]}

{{f[1, 1]} = y1[1] /. m}

{{f[1, 2]} = y9[1] /. n}

{{f[1, 3]} = y17[1] /. l}

{{f[1, 4]} = y25[1] /. c}

{{f[2, 1]} = y2[1] /. m}

{{f[2, 2]} = y10[1] /. n}

{{f[2, 3]} = y18[1] /. l}

{{f[2, 4]} = y26[1] /. c}

{{f[3, 1]} = y5[1] /. m}

{{f[3, 2]} = y13[1] /. n}

{{f[3, 3]} = y21[1] /. l}

{{f[3, 4]} = y29[1] /. c}

{{f[4, 1]} = y7[1] /. m}

{{f[4, 2]} = y15[1] /. n}

{{f[4, 3]} = y23[1] /. l}

{{f[4, 4]} = y31[1] /. c}

{d = Array[f, {4, 4}]}

{b = Det[d]}

{Print[ a, " " , z, " ", w, " ", b]}]]]