求助：Mathematica 求解微分方程组

Magic111

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

NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
/ {2 S# `5 ]$ Z, s. {4 f0 H8 G* ^0 n    omega,
. ~3 a; i7 v% K3 \  m# h3 p; L  [1 N  segma'[t] ==
   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
+ \3 L# X& O/ _& f( ^& p$ X  10}]
7 h" G4 r1 P0 S, p, l以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

f[EIS_] :=
  Re[Coef /.
    FindRoot[
! o8 U6 Y4 L" `6 ~  f     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
, @* t) k5 g# n; B& A


1 m  |# @8 K* ~( q8 _/ N; w# p5 \6 N: n程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。


$ @$ b" A- g6 G2 x& Q
% C" \$ D$ E3 G3 ~7 t* j8 p5 _: \十分感谢大神帮助

$ l. P" T% M' o3 G
# g# T+ d% w2 b* S( C

有生之莲

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