求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
& `5 L, Q) E; r9 l  c
* Q% J& P+ y0 y8 K求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
; p8 G$ Y/ w7 o0 q0 A5 k! M
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
& V; `0 v; E( a6 y/ V  I    omega,
% m8 ~$ ?4 K( N3 X3 O  Q( `  segma'[t] ==
2 l, i) E6 ?6 Q, s$ N  M2 F/ h! o   shear*str -
' {3 r( s( h3 O: U    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

7 r/ e  ?" A2 v/ }0 w* z: h7 J: ~f[EIS_] :=
  Re[Coef /.
- L' J7 T+ z! |1 p- k* ^2 u    FindRoot[
& v) C& z8 f2 S     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
; t( y1 w& T. T, I也就是f[EIS]每一点的值都需要解一个上述的方程。
# r; G5 d/ }" a+ M! p


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

7 a, F" \8 R) e( X  ?) E# _- {
9 G* [0 E4 L4 g* V* }
十分感谢大神帮助

. c1 F2 k  e$ r, R8 h- |# ?6 G

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
. ~4 \. n1 o6 q3 h' t; l