Methode de Runge-Kutta

Bonjour tout le monde! J'aurais besoin d'aide pour programmer (fortan) la methode de résolution des equations différentielles par Runge_Kutta ordre 4. Le problème est le suivant: on a dx/dt=P*(y(t)-x(t))
dy/dt=R*x(t)-x(t)*z(t)+y(t)
dz/dt=x(t)*y(t)-b*z(t) avec P=10, R=28, b=8/3
Je connais la methode générale, le souci est que t n'apparaît jamais directement dans les equations, du coup je ne sais pas trop comment calculer les coefficients... toute aide est la bienvenue...