版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

adjoinff

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 158
帖子: 132
在线: 2.9小时
虫号: 853277
注册: 2009-09-21
性别: GG
专业: 食品加工技术

[交流] 数值分析的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);    /*输出函数值*/
}
}

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2009-11-30 20:16:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

haixing2008

荣誉版主 (文坛精英)

应助: 22 (小学生)
贵宾: 2.252
金币: 32603.8
散金: 2020
红花: 44
沙发: 1
帖子: 25068
在线: 429.1小时
虫号: 540264
注册: 2008-04-06
性别: GG
专业: 计算数学与科学工程计算
管辖: 信息科学

★ ★ ★
小木虫(金币+0.5):恭喜抢沙发，给个红包
adjoinff(金币+2,VIP+0):谢谢帮助 12-1 21:24

main()这个前面要加个 void 即： void main()
t<=0.1001;这个好像也有问题
还有就是编程似乎不是很美观。
我也学得不好，仅供参考，期待专家解答

赞一下(2人)

回复此楼

平平淡淡才是真！

2楼2009-11-30 22:46:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

殛殛殛残酷

木虫 (著名写手)

应助: 3 (幼儿园)
金币: 3095.4
散金: 416
红花: 2
帖子: 1421
在线: 84.5小时
虫号: 853364
注册: 2009-09-21
性别: GG
专业: 计算数学与科学工程计算

★ ★
小木虫(金币+0.2):抢了个小板凳，给个红包
adjoinff(金币+1,VIP+0): 12-1 21:24

，我试一下！1

回复此楼

学术问题，交朋友

3楼2009-12-01 13:08:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9496.3
散金: 89
红花: 19
帖子: 3534
在线: 363.8小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

★ ★ ★ ★
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

赞一下(12人)

回复此楼

板凳要做十年冷文章不发一个字

4楼2009-12-01 15:49:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖