Cette question a déjà des réponses ici :
Commentaires
- @belisarius …! # @ * (& @ # ^%
- @ rm-rf Avez-vous besoin daide?
- Plus précisément, vous avez besoin de la réponse En utilisant le résultat des fonctions qui renvoient des règles de remplacement dans la question liée.
- Voir également: mathematica.stackexchange.com/q/9035/5
Réponse
Vous navez pas donné de valeurs numériques pour certains paramètres, alors jen ai inventé quelques-uns. Pour certains, vous pourriez ne pas obtenir de solutions, ou obtenir une solution complexe, cest donc quelque chose que vous pouvez examiner car je ne connais pas le physique s du problème.
Clear[x, t, b, varx, m, v0x]; y[t_] = x[t] /. First@DSolve[{b varx^2 - 2 b varx x"[t] + b x"[t]^2 + m x""[t] == 0, x[0] == x0, x"[0] == v0x}, x[t], t] donne
(b t varx - m Log[m/(v0x - varx)] + m Log[b t - m/(-v0x + varx)])/b Ensuite, vous pouvez utiliser la fonction y[t]
parms = {b -> 1, varx -> 2, m -> 1, x0 -> 1, v0x -> 0}; Plot[y[t] /. parms, {t, 0, 1}] 
D[y[t] /. parms, t] Out[48]= 2 + 1/(-(1/2) + t) etc …
Réponse
Suivre Nasser » réponse, voici une variante mineure:
x[tt_, {b_, varx_, m_, x0_, v0x_}] := Module[{}, x[t_, {b, varx, m, x0, v0x}] = Block[{x, t}, x[t] /. First@ DSolve[{b varx^2 - 2 b varx x"[t] + b x"[t]^2 + m x""[t] == 0, x[0] == x0, x"[0] == v0x}, x[t], t] ]; x[tt, {b, varx, m, x0, v0x}] ] Ensuite, vous pouvez lévaluer par:
parms = {1, 2, 1, 1, 0}; x[4, parms] Plot[x[t, parms], {t, 0, 1}] etc. Notez que lODE est calculé une seule fois pour chaque vecteur de paramètres.