일부 매개 변수에 대해 숫자 값을 제공하지 않았기 때문에 일부를 구성했습니다. 일부의 경우 솔루션을 얻지 못하거나 복잡한 솔루션을 얻을 수 있습니다. 약 s.
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]
제공
(b t varx - m Log[m/(v0x - varx)] + m Log[b t - m/(-v0x + varx)])/b
그런 다음 함수를 사용할 수 있습니다. 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)
기타 …
Nasser 팔로우하기 ” 의 답변, 다음은 약간의 변형입니다.
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}] ]
그런 다음 다음과 같이 평가할 수 있습니다.
parms = {1, 2, 1, 1, 0}; x[4, parms] Plot[x[t, parms], {t, 0, 1}]
등. ODE는 각 매개 변수 벡터에 대해 한 번만 계산됩니다.