求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
3 D8 k+ g+ G" {% f% ^
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
. d/ }3 E) c% ^' o2 L% \! o
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
" X3 H% ^3 C! B* q9 K- k6 a    omega,
  segma'[t] ==
   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

( k: G/ f* _/ n- @5 K* Zf[EIS_] :=
! s" C' Y) f8 G' B' u" @: F$ c8 ~  Re[Coef /.
+ x( O7 @0 j- c+ ], J8 |    FindRoot[
$ X% L8 u# z8 o4 o3 Q     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
  {2 k0 `  |8 y! \6 p- J也就是f[EIS]每一点的值都需要解一个上述的方程。
$ O8 C( B: y( \
4 i& z% Q+ A6 ?: s$ s/ f- ]! U

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



9 a) r: C4 w% W2 l十分感谢大神帮助
% o( ^3 y* ~8 Z3 [/ H
  \; ^" l6 b% K- e

有生之莲

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