GSL 集成失败

问题描述

我有一个程序可以对以下形式的常微分方程组进行积分：

dp0/dt = a - (1+gamma) * p0
dp1/dt = p0 - (1+gamma) * p1
dp2/dt = p1 - (1+gamma) * p2
dp3/dt = p2 - (1+gamma) * p3

我以通常的方式使用 gsl_odeiv_evolve_apply 作为集成引擎：

while ((*t) < t1)
{      
    status = gsl_odeiv_evolve_apply(e,c,s,&sys,t,t1,&h,y);
    if (status != GSL_SUCCESS)
    {
        printf("\n ************ status: %d\n",status);
        break; 
    }
}

while 循环失败，错误代码为 status = 32526。 p[i]和dp[i]/dt的值我已经一步步打印出来了，没有NAN也没有inf。

我如何知道导致此错误的原因？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

c gsl gsl