24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3359  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 求调剂,一志愿厦门大学,生物与医药,总分272,本科211 +3 Electron1cc 2026-04-01 4/200 2026-04-05 20:24 by lys0704
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 材料调剂 +13 一样YWY 2026-04-03 14/700 2026-04-05 18:20 by 蓝云思雨
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 17:59 by 猪会飞
[考研] 调剂 +3 好好读书。 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +4 一名26考研生 2026-04-04 4/200 2026-04-05 10:01 by barlinike
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 301求调剂 +18 骆驼男人 2026-04-02 18/900 2026-04-04 20:33 by 蓝云思雨
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 291求调剂 +4 迷蒙木木 2026-04-01 5/250 2026-04-04 15:59 by sihailian3
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 一志愿郑大材料工程290求调剂 +20 Youth_ 2026-03-30 20/1000 2026-04-02 14:48 by 5896
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
信息提示
请填处理意见