| 查看: 955 | 回复: 13 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
数值分析的c语言程序，帮忙看下，改改错
|
|||
|
http://d.namipan.com/d/79567bb4f ... dda82dcaf2d009e0100 http://d.namipan.com/d/79567bb4f ... dda82dcaf2d009e0100 我的题目已经放到米人里了，请朋友们帮忙看看。 先谢过了。 其中之一是下面的，大约有6个题，其他都传到米人了。 用定步长四阶Runge-Kutta法求解 dy1/dt=1 dy2/dt=y3 dy3/dt=1000-1000y2-100y3 y1(0)=0 y2(0)=0 y3(0)=0 方程组是上面六个式子 h=0.0005,打印yi(0.025), yi(0.045), yi(0.085), yi(0.1),(i=1,2,3) 我是这么编的 #include main() { double t,h=0.0005,y1=0,y2=0,y3=0,ky1[5],ky2[5],ky3[5]; for(t=h;t<=0.1001;t+=h) /*Runge-Kutta算法具体过程*/ { ky1[1]=h*1; ky2[1]=h*y3; ky3[1]=h*(1000-1000*y2-100*y3); ky1[2]=h*1; ky2[2]=h*(y3+ky3[1]/2); ky3[2]=h*(1000-1000*(y2+ky2[1]/2)-100*(y3+ky3[1]/2)); ky1[3]=h*1; ky2[3]=h*(y3+ky3[2]/2); ky3[3]=h*(1000-1000*(y2+ky2[2]/2)-100*(y3+ky3[2]/2)); ky1[4]=h*1; ky2[4]=h*(y3+ky3[3]); ky3[4]=h*(1000-1000*(y2+ky2[3])-100*(y3+ky3[3])); y1+=(ky1[1]+2*ky1[2]+2*ky1[3]+ky1[4])/6; y2+=(ky2[1]+2*ky2[2]+2*ky2[3]+ky2[4])/6; y3+=(ky3[1]+2*ky3[2]+2*ky3[3]+ky3[4])/6; if((t>0.0249&&t<0.0251)||(t>0.0449&&t<0.0451)||(t>0.0849&&t<0.0851)||(t>0.0999&&t<0.1001)) printf("y1(%f)=%f,y2(%f)=%f,y3(%f)=%f\n",t,y1,t,y2,t,y3); /*输出函数值*/ } } |
回复此楼» 猜你喜欢
267求调剂
已经有13人回复
-
292求调剂
已经有3人回复
-
322求调剂
已经有8人回复
-
348求调剂
已经有7人回复
-
材料调剂
已经有5人回复
-
349求调剂
已经有4人回复
-
环境工程 085701，267求调剂
已经有16人回复
-
一志愿北交材料工程总分358
已经有5人回复
-
材料0856 英一数二 323 求调剂
已经有8人回复
-
318求调剂
已经有4人回复
14楼2009-12-17 19:52:32
haixing2008
新虫
- 应助: 22 (小学生)
- 贵宾: 2.252
- 金币: 32603.8
- 散金: 2020
- 红花: 44
- 沙发: 1
- 帖子: 25068
- 在线: 429.1小时
- 虫号: 540264
- 注册: 2008-04-06
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 信息科学
2楼2009-11-30 22:46:29
殛殛殛残酷
管理员
- 应助: 3 (幼儿园)
- 金币: 3095.4
- 散金: 416
- 红花: 2
- 帖子: 1421
- 在线: 84.5小时
- 虫号: 853364
- 注册: 2009-09-21
- 性别: GG
- 专业: 计算数学与科学工程计算
3楼2009-12-01 13:08:17
bluesine
专家顾问
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9853.9
- 散金: 89
- 红花: 19
- 帖子: 3623
- 在线: 377.6小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
★ ★ ★ ★
formleaf(金币+2,VIP+0):谢谢参与！ 12-1 16:03
adjoinff(金币+2,VIP+0):谢谢帮助 12-1 21:24
formleaf(金币+2,VIP+0):谢谢参与！ 12-1 16:03
adjoinff(金币+2,VIP+0):谢谢帮助 12-1 21:24
你的这个程序可以运行啊。。。 y1(0.025000)=0.025000,y2(0.025000)=0.151599,y3(0.025000)=8.335426 y1(0.045000)=0.045000,y2(0.045000)=0.312860,y3(0.045000)=7.536280 y1(0.085000)=0.085000,y2(0.085000)=0.560581,y3(0.085000)=4.946353 y1(0.100000)=0.100000,y2(0.100000)=0.628881,y3(0.100000)=4.180993 Press any key to continue |
4楼2009-12-01 15:49:30
