请教求解刚性微分方程组的问题
现有一刚性比为2.5*10^4的线性微分方程组,请问有什么好的算法对其进行求解?即求解刚性微分方程组的方法有哪些?望高手赐教,谢谢! 不知你是要自己编程序求解呢还是只想求出结果,如果是后者,直接利用Matlab中的ode类函数即可求解,方便的很。回复
求解刚性微分方程常用吉尔方法,他有标准的源程序,你只要编写自己的微分方程组就行了,他可以求出数值解 Malab的 ode15s ode23s0de23tb这三个可以试下,专门求解刚性ODE的如果需要过程,可以选择DASSL求解器试下,它是用来求解刚性常微分方程和微分—代数方程的,
Fortran语言的源代码 http://www.netlib.org/ 有下载
页:
[1]