24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3207  |  回复: 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 的主题更新
信息提示
请填处理意见