Méthode Runge Kutta d'ordre 4 en MATLAB

Bonsoir à tous !
Depuis quelques temps je me suis mis aux méthodes de résolutions d'equo diff de façon numérique et evidemment je suis tombé sur la méthode de Runge kutta d'ordre 4 qui permet une approximation assez précise de la solution a l'aide de la formule :

yn+1=yn+ h/6(k1+2k2+2k3+4)

Je cherche actuellement un algorithme en MATLAB qui permettrait d'avoir en entrée une equo diff du premier ordre

Connaissez vous donc un site ou je pourrais trouver cet algorithme?


Merci d'avance.

Salut,
tu peux jeter un oeil sur
http://www.mathworks.fr/matlabcentral/fileexchange/
Sinon, les équations ne sont pas très compliquées et c'est un très bon exercice que de les implémenter soit même ;)