求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
  Y- K7 ?5 T% X! o9 {/ u
2 `3 e6 }7 W6 Y5 s. o! A; g! {求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
; ?4 ^5 n7 ?1 `6 J) G6 Q4 x6 S1 v  segma'[t] ==
& R( V( B% V, |- y/ I: ?; V& C) e1 ]   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
1 o: \" L3 s0 x3 T以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

f[EIS_] :=
9 b8 ^( R7 T4 ?+ e' Y. O: E* r* S  Re[Coef /.
! `6 K# k. K4 g' D    FindRoot[
9 f" L  ~& |: @$ @     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。

) }# l& |, o7 v  T5 c7 E% L1 W" w
3 p+ M+ o3 i' \% w1 i/ x6 J
- H( r3 {9 I& F& w6 s* a程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
* s+ G$ r. W/ }- [
: k3 ~0 M( W, D* }: v& ]. l

$ _$ S' Z$ ^* X$ _. M十分感谢大神帮助
4 M* k$ g. L% v5 ^7 \$ l

: d% E- V' _  E2 k' d  o

有生之莲

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