24小时热门版块排行榜    

查看: 3321  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 8/400 2026-03-16 13:42 by Dream007008
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
信息提示
请填处理意见