求助：Mathematica 求解微分方程组

Magic111

萌新求助如何求解一个非线性方程组：
3 V& T) c8 C5 Z( q# z: N
求助大神，我要解的非线性方程组不是规范的，无法直接套用NDSolve,也没有在Help中查到处理方法。原因是微分方程组中有一项表达式需要再解一个方程才能得到，具体的：
* x: [- X4 I' l
/ K5 C6 F. ]# h- d% A3 PNDSolve[{ EIS'[t] == (segma[t] + f[EIS[t]]/(str*omega))*str*
  k' X4 ~* z" [% m    omega,
* D$ ^& w3 A- w9 R1 y  segma'[t] ==
7 U6 h  U% C1 ^( }/ m) }   shear*str -
    1/(str*omega)*Df[EIS[t]]*EIS'[t]，EIS[0] == -4.842, segma[0] == 0,}, {EIS, segma}, {t, 0,
  10}]
0 f& `/ A" _4 ]; R. Q7 k  n以上是我想要求解的微分方程，形式很简单。其中各个函数除了f[EIS]以外均有直接的表达式。而f[EIS]需要由以下表达式得到：
# ~, p( U( e& U0 G/ x/ g
7 Y. p8 c4 J5 Q% df[EIS_] :=
  Re[Coef /.
' F$ j' U& O7 x- q% v    FindRoot[
  q- |3 G$ G. e9 [, |3 y     Integrate[P[EA, EIS, Coef], {EA, 0, Infinity}] == 1, {Coef, 1}]];
' S; J8 G3 V3 v9 B3 l& U* m也就是f[EIS]每一点的值都需要解一个上述的方程。


' b9 G  ]4 e7 C3 ?
2 T- V; @# R4 ^  H+ |. D# H程序报错是求解f[EIS]时得不到数值无法进行。请问这样的情况应该怎么处理呢？我也尝试了将f[EIS]也作为一个方程放入NDSolve中求解微分代数方程，但是也无法得到结果。



十分感谢大神帮助

有生之莲

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