Python reference

NumPy

import numpy as np
x = np.array([1.0, 2.0, 3.0])
t = np.linspace(0, 10, 101)

Plotting

import matplotlib.pyplot as plt
plt.plot(t, y)
plt.xlabel("Time")
plt.ylabel("State")
plt.show()

ODE solution

from scipy.integrate import solve_ivp
sol = solve_ivp(model, [0, 100], y0)

Random generator

rng = np.random.default_rng(123)
u = rng.uniform()
z = rng.normal()

Linear algebra

y = A @ x
values, vectors = np.linalg.eig(A)
solution = np.linalg.solve(A, b)