求助：Mathematica 求解微分方程组

Magic111

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

8 E4 V9 X+ O& _; A/ M$ _4 G求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
0 M3 I7 H" N- Q. x
3 T% u& A' F' c7 n0 ]. LNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
7 [  q7 _1 R& e. J    omega,
% \3 [  ]/ ~% Y7 n  }1 _( K  segma'[t] ==
8 s/ J; p% e) D/ e, q5 k8 P" M   shear*str -
# k$ [" ^( {; M" a: q9 K# _    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
1 S7 P6 Y/ {% J% z  10}]
. ^8 ~+ i" |9 Q! j以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

f[EIS_] :=
6 p3 e5 y' Q4 ]$ T0 J4 i6 ]. j  Re[Coef /.
    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。

0 \$ e$ `0 g' e
/ K: C% {6 H" G. r# K
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。

( X: Q: I, L& w! R# T! J# N

# I* B. ?7 ~2 Y5 |  r. m  j' n十分感谢大神帮助
) S$ t# j8 X+ c9 ^

有生之莲

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