Н. Д. Решить задачу в системе уравнений в математике

Question

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

   `bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb + 
    f\[Alpha]*\[Eta]f*wf - whe, (3/2)*ni[t]*k*
   Ti'[t] == (1 - \[Eta]d)*wd - wie + f\[Alpha]*(1 - \[Eta]f)*wf - 
   whi, nd'[t] = -nd[t]*nT[t]*\[Sigma], 
   nT'[t] = -nd[t]*nT[t]*\[Sigma], 
   n\[Alpha]'[t] = nd[t]*nT[t]*\[Sigma], Te[0] = 1, Ti[0] = 1, 
  nd[0] == nT[0] == \!\(TraditionalForm\`
   \*FractionBox[\(1.4447999999999998`*^26\), \(2\)]\), n\[Alpha][0] = 0}

    sol = NDSolve[bal, {Te, Ti, nd, nT, n\[Alpha]}, t]

вот ошибка.

NDSolve::deqn: Equation or list of equations expected instead of - 
((2.85474*10^-12 E^(-19.983 ((1<<1>>Plus[<<3>>]<<1>>Power<<1>> 
<<1>>])/Ti[t])^(1/3)) nd[t] nT[t])/(Ti[t]^(2/3) (1-(15.136 Ti[t]+4.6064 
Ti[<<1>>]^2-0.10675 Ti[<<1>>]^3)/(1000+75.189 Ti[<<1>>]+13.5 
Power[<<2>>]+0.01366 Power[<<2>>]))^(5/6))) in the first argument 
{4.8*10^-9 nd[t] (Te^\[Prime])[t]==900000000000000000-8.70051*10^-25 
(nd[t]+nT[t])^2 Sqrt[Te[t]]-(5.2266*10^46 <<1>>^<<1>> (11.92 
+1.69505*10^-9 <<1>>^3))/((nd[t]+nT[t]) (3.77 +<<21>> 
<<1>>+1.32084*10^-19 Power[<<2>>]))+(8.7331*10^17 (24-Log[Times[<<4>>]]) 
nd[t]^2 (-Te[t]+Ti[t]))/(1.09626*10^24 Te[<<1>>]+5.97059*10^20 
 <<1>>)^(3/2)+(5.152*10^-16 E^(-19.983 Times[<<2>>]^(1/3)) nd[t] nT[t] 
 (8/3 (4.32916*10^-7+Times[<<3>>])+64/9 Plus[<<2>>]^2))/((1+104/27 
    Plus[<<2>>]+64/9 Power[<<2>>]) (32+Te[t]) Ti[t]^(2/3) (1-Plus[<<3>>] 
 Power[<<2>>])^(5/6)),<<7>>,0}.`enter code here`

Answer 1

Я изменяю значение = на == и присваиваю значения всем неизвестным переменным и избавляюсь от оболочки TraditionalForm, чтобы превратить это в обычную дробь

k=1;\[Eta]d=1;wd=1;wie=1;wb=1;f\[Alpha]=1;\[Eta]f=1;wf=1;whe=1;whi=1;\[Sigma]=1;
ne[t_]:=2t+1;ni[t_]:=3t+2;
bal = {(3/2)*ne[t]*k*Te'[t] == \[Eta]d*wd + wie - wb + f\[Alpha]*\[Eta]f*wf - whe,
  (3/2)*ni[t]*k*Ti'[t] == (1 - \[Eta]d)*wd - wie + f\[Alpha]*(1 - \[Eta]f)*wf - whi,
  nd'[t] == -nd[t]*nT[t]*\[Sigma], 
  nT'[t] == -nd[t]*nT[t]*\[Sigma], 
  n\[Alpha]'[t] == nd[t]*nT[t]*\[Sigma],
  Te[0] == 1,
  Ti[0] == 1, 
  nd[0] == nT[0] == 1.4447999999999998`*^26/2,
  n\[Alpha][0] == 0};
sol = NDSolve[bal, {Te,Ti,nd,nT,n\[Alpha]}, {t,0,1}];
Plot[{Te[t],Ti[t]}/.sol[[1]],{t,0,1}]

Теперь подставимфактические значения для всех этих переменных и ваши фактические функции для ne [t] и ni [t] и посмотрите, что вы получите.