求助：Mathematica 求解微分方程组

Magic111

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

: i+ ^6 p0 J7 s求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
1 Y/ z. i6 J7 T) [: B% ?% Q
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
7 j, q& F5 d4 H, _! a8 N  ~2 {  segma'[t] ==
& ?( o/ D( v9 v9 c, X" N& F   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
, J& H3 v3 i" X6 ^; I' J( r以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
4 u0 c+ K& W7 V5 h9 J) X7 _& z) H8 _
f[EIS_] :=
0 Q$ v. P2 h- Y  g  Re[Coef /.
    FindRoot[
# j, F8 i. Y( a& b' O     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。

  x. {+ n, y9 n' {: ?
7 P5 ^4 s( k, D/ P- z
: |: Y" t+ o$ b- I程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。

9 r; m7 X4 p* D- u' m% w

十分感谢大神帮助
; x5 ?$ e/ r% j, [8 E) C4 i

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
) \9 p# R9 A8 G% j7 J