When I execute this it gives no result with no error.
NDSolve[{x1'[t] == z1[t] + (y1[t] - 0.1) x1[t] ,
y1'[t] == 1 - 0.3 y1[t] - (x1[t] x1[t]),
z1'[t] == -x1[t] - 1.5 z1[t],
x2'[t] == (a2[t] - 36) x2[t] + z1[t] + (y1[t] x1[t]) -
a1[t] x1[t] - 2 (x2[t] - x1[t]),
y2'[t] == (14 - b2[t]) y2[t] +
1 + (0.3 - b1[t]) y1[t] - (x1[t] x1[t]) - 2 (y2[t] - y1[t]),
z2'[t] == (36 - a2[t]) x2[t] y2[t] + (c2[t] - 5) z2[
t] + (c1[t] - 1.5) z1[t] - 2 (z2[t] - z1[t]),
a1'[t] == x1[t] (x2[t] - x1[t]) + 2 (a1[t] - 0.1),
b1'[t] == y1[t] (y2[t] - y1[t]) + 2 (b1[t] - 0.3),
c1'[t] == -z1[t] (z2[t] - z1[t]) - 2 (c1[t] - 1.5),
a2'[t] ==
x2[t] y2[t] (z2[t] - z1[t]) - x2[t] (x2[t] - x1[t]) -
2 (a2[t] - 36), b2'[t] == y2[t] (y2[t] - y1[t]) - 2 (b2[t] - 14),
c2'[t] == -z2[t] (z2[t] - z1[t]) - 2 (c2[t] - 5),
x1[0] == 1, y1[0] == 1, z1[0] == 1, x2[0] == 50, y2[0] == 10,
z2[0] == 10, a1[0] == 1, b1[0] == 1, c1[0] == 1, a2[0] == 30,
b2[0] == 10, c2[0] == 2}, {x1, y1, z1, x2, y2, z2, a1, b1, c1, a2,
b2, c2}, {t, 0, 20}]
Attachments:
