exercise4:
Solve the stiff ODE system using back_euler_lam().
First 4 values of Y are  [[0.        ]
 [0.15633494]
 [0.30891986]
 [0.4538982 ]]
Graphics saved as "exercise4.jpg"