求助：Mathematica 求解微分方程组

Magic111

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

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

# R- t3 e) Q& [4 ]: X8 z( @NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
2 Q6 n/ K6 @) g2 d, K  w; i; C' |    omega,
$ Q; E6 {3 a6 @5 Z7 z& ^  segma'[t] ==
   shear*str -
, a; Y  o( \1 r! g9 y" N    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
9 @! C4 ]. W. w$ t: W  10}]
  D% ~( }/ u9 X4 N2 O& h以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
: w% e. J8 E; f, Y$ d# K
f[EIS_] :=
/ f  F! J7 T$ ?# G. p$ d3 b  Re[Coef /.
4 w: x' ?9 z8 b2 M( h    FindRoot[
* [- ^5 p5 V: u     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。

* M! M2 l0 k) \3 z- C& Z5 l

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



$ C; W* _0 x! V  J5 W十分感谢大神帮助

有生之莲

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