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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 数值分析的c语言程序，帮忙看下，改改错
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 955  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

adjoinff

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

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

» 猜你喜欢
1楼 2009-11-30 20:16:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haixing2008

无虫
★ ★ ★
小木虫(金币+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的回帖

殛殛殛残酷

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★
小木虫(金币+0.2):抢了个小板凳，给个红包
adjoinff(金币+1,VIP+0): 12-1 21:24
，我试一下！1
回复此楼
学术问题，交朋友
3楼2009-12-01 13:08:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

超级版主

科苑小木虫

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★ ★ ★
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的回帖

adjoinff

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
有都改好的帮忙发到  这个邮箱里，谢谢
1986join@163.com
赞一下 回复此楼
5楼2009-12-01 21:23:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
本帖仅楼主可见
赞一下
6楼2009-12-02 14:48:31
已阅   申请数学EPI   回复此楼   编辑   查看我的主页

bluesine

超级版主

科苑小木虫

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by wfgwfgwfg at 2009-12-2 14:48:

不会吧.怎么弄的?

copy 的结果啊，不然还能怎么弄。。。
赞一下 回复此楼
板凳要做十年冷文章不发一个字
7楼2009-12-03 12:53:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

adjoinff

主管区长

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
没有人了吗。。。
回复此楼
8楼2009-12-09 21:51:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

formleaf

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★
Doctorcbw(金币+1,VIP+0):谢谢 12-10 08:17
adjoinff(金币+1,VIP+0): 12-11 08:55
引用回帖:
Originally posted by adjoinff at 2009-12-9 21:51:
没有人了吗。。。

可能最近大家都很忙吧，毕竟期末了，年终岁尾事情多啊！如果你不急，等忙过这阵子我们帮你看看！
赞一下(7人) 回复此楼
9楼2009-12-09 22:05:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

adjoinff

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
谢谢了，其实我也是用来交作业的
赞一下 回复此楼
10楼2009-12-11 08:55:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 adjoinff 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂0703 +4 周嘉尧 2026-03-31 5/250 2026-04-01 09:25 by fangnagu
[考研] 262求调剂 +7 励志一定发文章 2026-03-31 7/350 2026-04-01 09:14 by jp9609
[考研] 294分080500材料科学与工程求调剂 +12 柳溪边 2026-03-26 13/650 2026-04-01 01:37 by 1018329917
[考研] 调剂推荐 +7 清酒714 2026-03-26 8/400 2026-03-31 22:12 by 544594351
[考研] 一志愿南师大0703化学 275求调剂 +5 Ripcord上岸 2026-03-27 5/250 2026-03-31 19:52 by mg1014
[考研] 346求调剂 一志愿070303有机化学 +11 萝卜炖青菜 2026-03-28 12/600 2026-03-31 17:54 by xhai2011
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 281求调剂 +5 亚克西good 2026-03-26 7/350 2026-03-30 20:42 by dophin1985
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 324求调剂 +9 hanamiko 2026-03-26 11/550 2026-03-30 14:27 by JourneyLucky
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 考研调剂 +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[考研] 材料与化工304求B区调剂 +4 邱gl 2026-03-26 7/350 2026-03-30 08:39 by 探123
[考研] 本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 打过很多竞赛，085406控制工程300分，求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭