求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：

8 j7 S; v( u8 k5 \7 J求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

) Z# K) M% |: v1 |, YNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
  \0 V: e* ?! P. d1 d, ^, o    omega,
" M3 p. c' I. G. f+ X/ l1 d2 M  segma'[t] ==
% e) [8 s! X- S5 P0 D   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
+ n0 y. U2 I8 m+ \* y$ \以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

/ m7 r  L2 e; Af[EIS_] :=
  Re[Coef /.
( ^7 E6 \% s; J/ e  F$ W    FindRoot[
1 B5 b7 K7 |/ u+ \4 {' b  u5 h2 p* w     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
, D8 `' i% Q& j


; G9 j' A6 U& h, F2 D" J( B  {程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。



十分感谢大神帮助


% ?8 H8 C, s! e% z/ u8 z  H

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
7 T  f/ |; r$ W6 g3 ]) ~5 ~