求助：Mathematica 求解微分方程组

Magic111

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

. J7 c* ]% j( k/ ?, o求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
    omega,
% c" Z" w; ~' R5 ?3 i9 \  segma'[t] ==
   shear*str -
# S7 I7 J5 [+ m; V    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
. M0 J# j5 W" G; G# f  10}]
/ w, [! R, w3 f; v3 l以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：

% t/ n. d9 O* u. d8 [! vf[EIS_] :=
1 J/ w: R( L8 t8 E5 w  Re[Coef /.
; E$ ?* }' Q4 `" |5 I2 ~  s0 {    FindRoot[
  {& N* O$ u; q& d     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
$ q7 c1 s1 }0 p
) H# r% r  }* v% K

2 |" V, j6 f; G9 C+ E程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。


3 F% X- B8 ]  K( t6 J6 X
6 J9 [  {9 G5 o) a) G& L; Q十分感谢大神帮助
/ I/ E5 M4 J' ^, v/ O  R

有生之莲

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