Bonjours à tous
serait il possible d'avoir un script RK4 pour résoudre le système suivant:
dx/ds=cos(teta)
dz/ds=sin(teta)
d(teta)/ds=2*5,838-13,45z-sin(teta)/x

cond initiales:
pour s=0 sin(teta)/x=5,838
x=0
z=0
teta=0

ceci servant a obtenir le profil d'une goutte pendante dont s est l'abscisse curviligne
merci