求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
" I7 d; E' l5 E& x# E7 G( X/ E
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

/ g% r6 F6 B/ q+ k* c" n# sNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
" v# z. N. O1 n6 d( [8 R( H    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]需要由以下表达式得到：
6 u% d' r5 \& s
& s# q( l4 p) J* V3 z+ E3 {f[EIS_] :=
  Re[Coef /.
    FindRoot[
2 v# F$ k+ c( S* Y     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
* {6 e/ \4 ^. S7 i& T! o
) R+ h7 @4 w& r5 j3 a/ ?
% e, p7 W3 Z" Y- q& e
# m! C7 w2 k6 I6 @+ b" l4 {% `程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
8 H1 S/ r% r5 k6 M( d  c" G

3 P% I$ O6 v7 V5 B+ F
十分感谢大神帮助
% l- C" k' h+ l) J/ m+ l

, v; \8 I, e4 F; T6 @) m0 S

有生之莲

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