Leierkastenmann
Прошу прощения но что-то опять не так. График у меня не строится. Получается прямая, вместо закручивающейся окружности.
Вот код:
Код:
myeps = 0.1;
mymu = 2;
myM = 1;
sol = NDSolve[
{
x'[t] == -y[t] +
Which[x[t] + z[t] < -myeps, -myeps^3, x[t] + z[t] < myeps,
mymu*(x[t] + z[t]), x[t] + z[t] > myeps, myeps^3],
y'[t] ==
x[t] + Which[x[t] + z[t] < -myeps, -myeps^3, x[t] + z[t] < myeps,
mymu*(x[t] + z[t]), x[t] + z[t] > myeps, myeps^3],
z'[t] == -z[t] +
Which[x[t] + z[t] < -myeps, -myeps^3, x[t] + z[t] < myeps,
mymu*(x[t] + z[t]), x[t] + z[t] > myeps, myeps^3],
x[0] == 0, y[0] == -0.8165, z[0] == 0
},
{x, y, z},
{t, 0, 1256.6371},
Method -> {"ExplicitRungeKutta", "DifferenceOrder" -> 4},
InterpolationOrder -> All,
StartingStepSize -> 0.0001
]
ParametricPlot[{x[t], y[t]} /. sol, {t, 0, 1256.6371},
PlotRange -> All]
еще предупреждение кидает:
Код:
NDSolve::mxst: Maximum number of 10000 steps reached at the point t == 166.6857351451608`. >>
Но t нужно больше!
