Vous n'avez besoin que de l'EDO d'ordre 1 pour tracer le champ de directions. Avec y' = f(x,y) vous avez la valeur de y' pour tous les points (x,y) que vous voulez vérifier. À la main, c'est un peu long, donc il est utile d'utiliser la fonction quiver de Matlab.
https://www.mathworks.com/help/matlab/ref/quiver.html
https://www.mathworks.com/help/matlab/ref/quiver.html