求助：Mathematica 求解微分方程组

Magic111

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

7 u! P3 l& i& `$ z求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

8 W- e# ^+ X& v3 d- bNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
5 W: `2 f# o* Y6 [9 b/ [  segma'[t] ==
- j& J4 @/ h* }* S/ @   shear*str -
( `$ [/ ]. ?! L% p2 i6 x    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
+ R% x, U) d  a/ U. H2 |1 G以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

f[EIS_] :=
3 {. a$ f) n9 M2 }, C) f( A  Re[Coef /.
4 H% F! @8 g4 o- j7 `    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
) V1 ~9 _% J* d4 g3 h4 G3 `也就是f[EIS]每一点的值都需要解一个上述的方程。

! H/ U- A+ L' ]  L6 {$ R

程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
6 p# j$ y! }5 u# ~
9 V1 ?9 v6 G* t8 w7 C: d8 V
) C0 v# e9 D1 s) S3 q' x
% r1 i* b: `! h7 h. k  |  F十分感谢大神帮助


2 g0 B5 M* U1 f) K) O* C! M; f

有生之莲

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