Mathematica: задача с решением дифференциального уравнения - PullRequest
0 голосов
/ 28 января 2020

У меня есть проблема, я работал над решением двух различных дифференциальных уравнений, и в данный момент я начинаю подозревать, что данный ответ был неправильным, я думаю, что это моя вина, но я не вижу, где находится Моя ошибка.

Сначала давайте определим некоторые значения:

c = -0.884

H0 = 68

G = 6.67*10^{-11}

Om = 0.3

Rhom0 = (3 H0^2/(8 Pi G))*Om

Это первое дифференциальное уравнение

co22 = NDSolve[{(1 + x)^3D[p[x]/(1 + x)^3, x]+ 6*c(H0* ( (1 + c/1.5) (1 + x)^{3/2} - (c/1.5) ) )^2 == 0, p[0] == Rhom0}, p, {x, 0, 10}]

И второе дифференциальное уравнение:

sol25 = DSolve[{(1 + x)^3*(-((3 p4[x])/(1 + x)^4) + p4'[x]/(1 + x)^3) == 0, p4[0] == Rhom0}, p4, {x, 0, 10}]

И когда я строю решение обоих, они дают мне один и тот же график, когда я впервые думаю, что это может быть совпадением, но затем я ввожу третье дифференциальное уравнение

sol24 = DSolve[{(1 + x)^3*(-((3 p3[x])/(1 + x)^4) + p3'[x]/(1 + x)^3) == -6*c*(H0*(0.3*(1 + x)^3 + 0.7)^(1/2))^2, p3[0] == Rhom0}, p3,{x, 0, 10} ]

И снова тот же ответ, вы знаете, в чем проблема?

...