应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 请教一个问题
41/1返回列表
查看: 395  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

言论压力

至尊木虫 (职业作家)

我是个好人

[求助] 请教一个问题

请问下面这个定步长四阶龙格-库塔法程序为什么 不能在matlab运行，该怎么运行。
#include "stdlib.h"
  void rkt1(t,y,n,h,k,z,f)
  void  (*f)();
  int n,k;
  double t,h,y[],z[];
  { int i,j,l;
    double a[4],tt,*b,*d;
    b=malloc(n*sizeof(double));
    d=malloc(n*sizeof(double));
    a[0]=h/2.0; a[1]=a[0];
    a[2]=h; a[3]=h;
    for (i=0; i<=n-1; i++) z[i*k]=y;
    for (l=1; l<=k-1; l++)
      { (*f)(t,y,n,d);
        for (i=0; i<=n-1; i++) b=y;
        for (j=0; j<=2; j++)
          { for (i=0; i<=n-1; i++)
              { y=z[i*k+l-1]+a[j]*d;
                b=b+a[j+1]*d/3.0;
              }
            tt=t+a[j];
            (*f)(tt,y,n,d);
          }
        for (i=0; i<=n-1; i++)
          y=b+h*d/6.0;
        for (i=0; i<=n-1; i++)
          z[i*k+l]=y;
        t=t+h;
      }
    free(b); free(d);
    return;
  }
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-23 10:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

感谢参与，应助指数 +1
这个是c语言程序，当然不能在matlab底下运行。
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
2楼2013-09-23 11:02:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

言论压力

至尊木虫 (职业作家)

我是个好人
引用回帖:
2楼: Originally posted by baobiao007 at 2013-09-23 11:02:07
这个是c语言程序，当然不能在matlab底下运行。

怎么让C语言的程序在matlab下运行
赞一下 回复此楼
3楼2013-09-23 11:05:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

引用回帖:
3楼: Originally posted by 言论压力 at 2013-09-23 11:05:54
怎么让C语言的程序在matlab下运行...

看懂这个程序，然后翻译成matlab程序，或找个龙哥库塔的matlab程序
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
4楼2013-09-23 11:36:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 言论压力 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭