- 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一元方程数值解法
第8章 一元方程数值解法 * * * * program Read(*,*)x,m Do 10 i=1,m X1=(-x**3-2*x*x-1)/2 If(abs(x-x1).gt.1e-6) then X=x1 Else Stop End if 10 continue Write(*,*)x end 程序运行情况如下: 1,20 i=1 x1=-1.0000000 i=2 x1=-1.5000000 i=3 x1=-1.5625000 . . . i=19 x1=-1.5436890 Stop(循环19次后,已满足所给条件) program Read(*,*)x N=1 10 x1=x F=x1**3-2*x1**2+4*x1+1.0 F1=3*x1**2-4*x1+4 X=x1-f/f1 Write(*,*)n,x1,x N=n+1 If(abs(x-x1).gt.1e-6)goto 10 end 运行过程记录如下: 1.0 N=1 x1=1.0000000 x=-0.3333333 N=2 x1=-0.3333333 x=-0.2287582 N=3 x1=-0.2287582 x=-0.2225152 … … … … N=5 x1=-0.2224945 x=-0.2224945 program Read(*,*)x1,x2 F1=x1**3-6*x1-1 F2=x2**3-6*x2-1 10 x=(x1+x2)/2 F=x**3-6*x-1 If(sign(f,f1).eq.f) then X1=x F1=f Else X2=x F2=f End if If(abs(x1-x2).gt.1e-5).and. $ (abs(f).gt.1e-6)) goto 10 If(abs(f).gt.1e-6) x=(x1+x2)/2 Write(*,*) x end 程序运行结果: 1.0, 5.0 x=2.5289170 2.0, 4.0 x=2.5289170 10 Read(*,*)x1,x2 F1=x1**3-2*x1**2+7*x1+4 F2=x2**3-2*x2**2+7*x2+4 If(sign(f1,f2).eq.f1) goto 10 F=1.0 20 if((abs(x1-x2).gt.1e-5).and.(abs(f).gt.1e-6)) then x=x2-(x2-x1)/(f2-f1)*f2 f=x**3-2*x**2+7*x+4 If(sign(f,f1).eq.f) then X1=x F1=f Else X2=x F2=f End if Goto 20 End if If(abs(f).gt.1e-6) x=(x1+x2)/2 Write(*,*) x end 运行结果如下: -1,1 X=-0.4871201 -5 5 X=-0.4871201
文档评论(0)