Gratis flashcards

Man definerer høyresiden f(t, y) og kaller $solve_ivp$:

from scipy.integrate import $solve_ivp$

import numpy as np

def f(t, y):

return -2*y

sol $= solve_ivp(f, t_span=(0, 5), y0=$[1]$, t_eval=np.linspace(0, 5, 100))$

Første argument er funksjonen $f(t, y), t_span$ er $(t_start, t_slutt), y0$ er en liste/array med startverdier (her [1]), og $t_eval$ angir hvilke t-verdier løsningen skal lagres i. Returverdien sol er et objekt der sol.t er tidspunktene og sol.y er løsningen (en 2D-array med én rad per komponent, her sol.y[0]). Standardmetoden er RK45 (adaptiv Runge-Kutta av orden 5(4)). Den eksakte løsningen er y(t) = e^(-2t).