Question

我还是很新的与数学的工作，我要在这个有很多，我真的不明白我是什么（很清楚）做错了。

我基本上是想重新获得一个获得三个函数的结果，我称之为C1，C2和f：

我将约束用于C1和C2的导数以及f的初始值。我需要照顾其他的我的事情是确保当f * 22 = 1，创建零一个部门，我会做出C1' 变为零，整个方程式。（技术上的理论是，你定义LAMBDA / P是东西是对称于方程x的特定值，其中此奇点出现的第二项的其余部分）。

我为这一点，但由于非常遗憾提前。

c = 2; \[Lambda] = 1; k = 1; Dp = 0.5; pm = 1;

NDSolve[{C1'[x] == pm*((1 - f[x]*C2[x]^2)^(-1/2) )* (C2[x]*((C2[x] - c)^2 - C1[x]^2) + 2*(\[Lambda]/k)*C2[x] + (Dp/k) - 2*(1 - f[x]*C2[x]^2)*(C1[x] - C2[x])/f[x]), 
C2'[x] == pm*2*((1 - f[x]*C2[x]^2)^(1/2))*(C1[x] - C2[x])/f[x], 
f'[x] == pm*4*(1 - f[x]*C2[x]^2)^(1/2),
C2'[0.] == (1/3)*C1'[0.], C2'[1] == (1/3)*C1'[1], f[0.] == 0.1, 
WhenEvent[C2[x]^2 == 1/f[x], {C1'[x] == 0.}]},
{C1[x], C2[x], f[x]}, {x, 0., 1}, Method -> {"ExplicitRungeKutta", "DifferenceOrder" -> 4}]

数学中三个微分方程的NDSolve问题

0 个答案: