Qurak

Nd solve event locator limitations

Підручник 3 кроків, що виконуються по порядку в одній сесії. Кожен крок було виконано на рушії, і вивід нижче — те, що він дав.

Крок 1
Block[{n = 0}, NDSolve[{y'[t] == y[t], y[-1] == E^-1}, y, {t, 5}, Method -> {"EventLocator", "Event" -> Sin[π y[t]], "EventAction" :> n++}];n]
Вивід
1
Крок 2
Block[{n = 0}, NDSolve[{y'[t] == y[t], y[-1] == E^-1}, y, {t, 5}, Method -> {"EventLocator", "Event" -> Sin[π y[t]], "EventAction" :> n++}, MaxStepSize -> 0.001];n]
Вивід
1
Крок 3
Block[{n = 0}, NDSolve[{y'[t] == y[t], y[-1] == E^-1, z'[t] == D[Sin[π y[t]], t], z[-1] == Sin[π E^-1]}, {y, z}, {t, 10}, Method -> {"EventLocator", "Event" -> z[t], "EventAction" :> n++}, MaxSteps -> ∞];n]
Вивід
1

Використані функції

Пов'язані рецепти

Усі рецепти · Довідник функцій · Використати з MCP-клієнта