求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
8 u$ c2 @$ U) d# [
" D; E2 ^2 D9 Y" `& }- K5 @7 S8 y& S) D求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

8 Y; _- b4 h  R3 E+ S0 DNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
  segma'[t] ==
: E( h0 x( F0 v6 _- p   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
5 V7 S  Z( d$ M! Y, z! A) D! A  10}]
- \% K* v, p6 g4 W, d以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

9 s/ n) B$ B2 y# n5 j" c( of[EIS_] :=
  Re[Coef /.
% I$ U5 j) `! }% Z    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
# O3 h( `7 ]; G8 ?6 q; y也就是f[EIS]每一点的值都需要解一个上述的方程。
' x8 ]; A" E2 Y2 `! [+ s. M1 ?
* g9 x& f. l+ b' a: c

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

6 m7 T. K* a: _/ r. t
9 m, g* C  u2 q! p: Z- E
十分感谢大神帮助

* f8 {# ?/ c& E/ B: [: x

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
( Y* f+ N+ e8 N! k