In this section, we will consider the mathematical task of numerically solving a system of ordinary equations for given initial values:
y'(t) = f(t, y) y(t0) = y0∈ ℝn
The solution of this problem is a function y. A numerical method aims at computing good approximations, yi≈ y(ti) at discrete points, the communications points ti, within the interval of interest [t0, te]. We collect the data that describes the problem in a class, as follows:
class IV_Problem: """ Initial value problem (IVP) class """ def __init__(self, rhs, y0, interval, name='IVP'): """ rhs 'right hand side' function of the ordinary differential equation f(t,y) y0 array with initial values interval start and end value of the interval of independent variables often initial and end time name descriptive name of the problem...