求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
* ^+ `8 p( @9 d0 z
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
  segma'[t] ==
7 q; z" U0 J5 q  I4 n/ ?% S   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
8 C8 ~5 A5 @' j4 o以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
0 V% z  C, ?- y, F: x
f[EIS_] :=
  Re[Coef /.
    FindRoot[
9 q  v, V7 C6 e$ v     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
9 e& l2 B4 W9 X. V也就是f[EIS]每一点的值都需要解一个上述的方程。
1 x' I% e4 q+ V# _8 X; F" S

; r8 j3 H6 V7 l( Q
0 u, O) S6 f3 P" d$ [程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
* I6 @& d. i2 _
$ L% V5 W  Z7 z' E" T: d; O

& ]8 G7 n) o- B% _  f2 |4 c- T十分感谢大神帮助

! ?6 k9 r' {" f, u

有生之莲

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