| 查看: 477 | 回复: 3 | ||
言论压力
|
[求助]
请教一个问题
|
|
请问下面这个定步长四阶龙格-库塔法程序为什么 不能在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; } |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有65人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个吸附问题
已经有5人回复
- 请问一个关于DIBALH的问题
已经有5人回复
- 请教接一个接收函的问题
已经有12人回复
- 请教一个HWE反应的问题
已经有7人回复
- 请教关于假发票的一个问题！
已经有10人回复
- 请教一个积分用留数定理计算的问题
已经有15人回复
- 请教一个关于卡方检验样本规模的问题
已经有5人回复
- FFT算法的一个问题请教
已经有18人回复
- 请教关于不饱和聚酯树脂的一个问题
已经有12人回复
- 请教一个免疫共沉淀技术的问题
已经有4人回复
- 请教一个手性物质的旋光性问题
已经有12人回复
- [求助]外行请教一个简单的生物教科书题目！！！！
已经有3人回复
- 公派出国德国签证一个小问题请教
已经有14人回复
- 请教大家一个问题--提质粒
已经有8人回复
- 请教大家一个审稿意见问题~
已经有3人回复
- 【求助】请教一个美国签证问题：）
已经有9人回复
- 请教文章修回的一个问题
已经有7人回复
- 【求助】请教一个材料力学问题
已经有30人回复
- 【请教】请教一个关于图片精度的问题
已经有6人回复
- 【求助/交流】请教一个农杆菌转导的问题
已经有4人回复
- 【求助】再请教大家一个问题
已经有16人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
2楼2013-09-23 11:02:07
言论压力 
至尊木虫 (职业作家)
我是个好人
- 应助: 19 (小学生)
- 贵宾: 4.585
- 金币: 13152.6
- 散金: 1190
- 红花: 5
- 沙发: 9
- 帖子: 3303
- 在线: 410.7小时
- 虫号: 678275
- 注册: 2008-12-20
- 性别: GG
- 专业: 高分子合成化学
3楼2013-09-23 11:05:54
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
4楼2013-09-23 11:36:39
