求助：Mathematica 求解微分方程组

Magic111

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

+ m- K. a- ]. a' _9 A$ G求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：

NDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
; X- `' V7 Y9 u: h( {* s    omega,
  segma'[t] ==
3 `5 ?; X' H& a( ?# h" Q   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
, c! p( ?1 e0 z' O1 f8 M  10}]
% |. x# e) e7 n: t以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
* C" Y. L" s4 c. ?; A
6 \+ T7 _7 ~7 V- |f[EIS_] :=
  Re[Coef /.
- _8 V# B' t: b9 T. O    FindRoot[
     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
也就是f[EIS]每一点的值都需要解一个上述的方程。
& r7 H6 C8 R0 {! M' M. X2 B+ }- `
, v% P9 p3 g' {1 P& o) N1 {6 g
1 j3 E" W( e9 x2 x0 P# j% h
程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。
" Q/ d1 H- B$ m$ v
" b+ N3 p) p/ w8 c0 z; Z
; s7 z! G4 n# T5 n, T
. h& A' ~3 S- |. c! x十分感谢大神帮助

6 q! G+ P0 S, V
! Z8 {9 q+ q9 I/ \1 c" z

有生之莲

泛函就是没办法。。。尝试重新转化处理问题吧
7 z+ U9 i5 |5 b7 ~7 a2 W4 n