I've successfully (well I think so) made 2 different programs that can numerically solve an ODE using Euler and Runge-Kutta's methods.
Here they are:
program test
implicit none
real(8)::a,b,h,y_0,t
write(*,*)"Enter the interval a,b, the value of the step-size h and the value of y_0"...