求助：Mathematica 求解微分方程组

Magic111

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

求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
& _# K3 z! ^6 r
NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
9 R5 n& I/ ?) p( o: x8 h. h/ Z& l    omega,
) g% u: v. w. k* w  segma'[t] ==
   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
8 M  Z+ ?- r4 S) W9 J0 N- b8 M  10}]
. ]* d: }  l! b5 @4 D以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
& C) Z$ q: \8 M: h/ s" I
) @# e! P) a  a5 N. ~: cf[EIS_] :=
% }5 `# ^9 x( E$ K  Re[Coef /.
    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
' s# C+ o6 h( a$ p9 b

) T: z/ W. u9 I+ S0 _
, u0 w2 s- @. T5 Q+ n  A1 b程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
( \( P: t8 X& H8 Q8 K1 @7 ?, ?; y


& Y, M* T! W% O7 n# [十分感谢大神帮助
1 p3 D8 {+ j: `, u* A

有生之莲

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