# Résolution de l'équation différentielle définissant, en polaires dans # le plan de développement, la chaînette conique. # a: rho, b: rho' # ---------------------------------------------------------------------- # Usage: ode < chainette.ode > chainette.dta # ---------------------------------------------------------------------- # JMS - 25/11/07 lambda = 3 a' = b b' = (a*a*(2*a-lambda)+b*b*(3*a-2*lambda))/(a*(a-lambda)) a = 2.414 b = 0 print t, a step 0, 12.57