求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
/ H3 I, C) j5 H
0 x* [: Y4 I. S; Q. L2 A求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
: H! k! p+ d( v5 ?5 h7 b; k+ G; x% t
8 _- t% K' k! n9 y3 H* R1 yNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
  segma'[t] ==
8 `3 w7 g. e' N, L. ^# d   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
' X' R# c8 N# m- B  10}]
  l2 @* y2 g0 ]5 ~; Y以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
5 h' m' N; s' t3 I0 k" [9 A+ h
f[EIS_] :=
  Re[Coef /.
. m! N9 M$ H; t4 b' g2 |    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。

: |  v" N! W. ^6 M. w3 K# D( E) D
6 @: b4 i/ R# M* }0 k: D
: f( k# Q, Y) E7 i程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
8 C& ?6 ]% t9 [) y% S# q3 ^& j& Z

, O! H& t; y9 @* V* f
# K. T/ i. K! t; i9 r. E5 [0 s十分感谢大神帮助


6 T" @" M. i4 C- L3 ]! v

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
5 A# w, X2 w- ^