求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
: l1 {/ w0 ^4 P8 r8 X, Q
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

" S" ^& r* ]. B' x, yNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
% R- w" f0 H& Y- S3 D& r    omega,
* u: g. a( b) k; |5 d9 }" G7 ?) o  segma'[t] ==
/ i! \+ ]( o! v. c5 A; V/ L+ w2 {   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
2 _6 k% L$ B" N* m/ N
f[EIS_] :=
& n: {; l. p9 H: p5 T+ G6 n+ k  Re[Coef /.
    FindRoot[
2 A, R$ f" `" s0 K, H     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
- |: Q) k, m7 h) H也就是f[EIS]每一点的值都需要解一个上述的方程。


4 S+ Q% j6 P# V' n6 L
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。

( v  s- j- K1 D

/ e2 g+ B! L! k5 x- m& B十分感谢大神帮助
  O) J4 ]# X6 p4 ?7 \

" C: d9 T2 Q2 _" u# _4 ?5 B

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧