应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 一元高次方程组解法
81/1返回列表
查看: 3388  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

sddhywj2008

木虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-25 17:55:59
已阅   回复此楼   关注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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 297，工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
信息提示
关闭
请填处理意见
关闭