Hello everyone,

I'm a PhD student and I have to take some post-grad courses. One is about numerical methods to solve differential equations.

I understood methods such as the Euler, Backward Euler, etc... fine, BUT this kind of sytem I can't solve numerically:

(A, B, C, D some constants) We know X(0), Y(0) and Z(0)

X' = -A*X + B*Y*Z

Y' = A*X - B*Y*Z - DY'

Z' = DY'

My main problem comes from the multiplications between the function (Y*Z here).

I tried to do X' + Y', this kind of trick, to simplify the system but it led to no solution, and anyways I have other exercises that do not seem to simplify either, yet it should be possible to solve them with Euler method...

I know how to solve an equation such as X' = AX + B in a numerical way, but for this system I have no clue how to do it...

If anyone can help, thank you very much!