24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3360  |  回复: 7

sddhywj2008

木虫 (小有名气)

[求助] 一元高次方程组解法

本人刚接触程序,想求助各位达人用FORTRAN如何求解一元高次方程组,比如X的15次方或者30次方等。金币重谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

li2008jie

铁杆木虫 (著名写手)

逍遥孤鹤

【答案】应助回帖

感谢参与,应助指数 +1
用matlab应该好解,可以尝试一下。 fortran语言不会。
阿甘与徐三多,执着+执着+执着=成功!
2楼2013-11-25 18:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
用Fortran写个牛顿法的程序求解。
3楼2013-11-25 19:07:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sddhywj2008

木虫 (小有名气)

引用回帖:
3楼: Originally posted by sciencejoy at 2013-11-25 19:07:13
用Fortran写个牛顿法的程序求解。

谢谢你,不会写呢。有现成的吗?呵呵
4楼2013-11-25 19:45:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)

引用回帖:
4楼: Originally posted by sddhywj2008 at 2013-11-25 19:45:57
谢谢你,不会写呢。有现成的吗?呵呵...

数值方法的书上都有。网上搜也能搜得到。
5楼2013-11-25 20:37:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sddhywj2008: 金币+20, ★★★很有帮助 2013-12-08 10:44:16
给你一个例子,使用Newton法求解方程,比如求解 f(x)=x^N-x-a=0, N 可以选大于2的 任意整数,a>0, 选取初值 x= 1.1


!flag
subroutine Newton7
  implicit none
  integer :: iter,N
  real*8 :: a,x,dx,f,f1,eps=1.d-10

  ! specify parameters
  N=20
  a=6d0

  ! initial value, carefully chosen, otherwise not convergent
  x=1.1d0

  iter=0
9 iter=iter+1

  ! f is the function, f1 is the derivative
  f=x**N-x-a
  f1=N*x**(N-1)-1
  dx=-f/f1
  x=x+dx
  print"(a,i3,9(a,f20.12))",' iter=',iter,' dx=',dx,' x=',x
  if(abs(dx) > eps) goto 9
  print*,''
  stop

end subroutine Newton7
6楼2013-11-26 10:22:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sddhywj2008

木虫 (小有名气)

引用回帖:
6楼: Originally posted by pippi6 at 2013-11-26 10:22:36
给你一个例子,使用Newton法求解方程,比如求解 f(x)=x^N-x-a=0, N 可以选大于2的 任意整数,a>0, 选取初值 x= 1.1


!flag
subroutine Newton7
  implicit none
  integer :: iter,N
  real*8 :: a,x,dx ...

谢谢您。因为的方程比这个复杂,可能不能直接把导数写出来。
7楼2013-11-26 14:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

isityou

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
7楼: Originally posted by sddhywj2008 at 2013-11-26 14:19:29
谢谢您。因为的方程比这个复杂,可能不能直接把导数写出来。...

不能求导数,那就用二分法吧,只是效率慢点,或者其他解法。
8楼2013-11-27 09:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sddhywj2008 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 316求调剂 +5 yyx想调剂 2026-04-05 5/250 2026-04-05 22:22 by 咔咔咔咔9
[考研] 一志愿河北工业大学材料工程,初试344求专硕调剂 +3 15933906766 2026-04-05 3/150 2026-04-05 22:17 by dongzh2009
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 化学0703-一志愿211-338分求调剂 +7 vants 2026-04-05 7/350 2026-04-05 18:17 by cql1109
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +16 披星河 2026-04-03 16/800 2026-04-05 09:00 by dick_runner
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 考研调剂 +4 美丽的youth_ 2026-04-04 5/250 2026-04-04 17:16 by imissbao
[考研] 0703求调剂 +6 zizimo 2026-03-31 6/300 2026-04-04 14:16 by 无际的草原
[考研] 化工调剂303分,过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 22408 266求调剂 +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by 笔落锦州
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 302求调剂一志愿北航070300,本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 生物学327,求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
信息提示
请填处理意见