How can I program this algorithm with Python

[ichabod801]

What have you tried? We're not big on writing code for people here, but we would be happy to help you fix your code when you run into problems. When you do run into problems, please post your code in Python tags, and clearly explain the problem you are having, including the full text of any errors.

That said, this is just a couple for loops and some basic math. x sub i just becomes x[i]. So, x[i + 1] = x[i] + h * (s + r * y[i]) / (...) ...