j'ai besoin de générer des conditions initiales modélisant des bras spiraux galactiques.

je pars de l'équation polaire suivante :

rho=a./(log(b*tanh(theta/(2*n))) avec a, b et n des paramètres...