| 查看: 3535 | 回复: 7 | ||
sddhywj2008木虫 (小有名气)
|
[求助]
一元高次方程组解法
|
| 本人刚接触程序，想求助各位达人用FORTRAN如何求解一元高次方程组，比如X的15次方或者30次方等。金币重谢。 |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有12人回复
-
2026国自然函评费到账
已经有3人回复
-
科研孤儿太难了
已经有17人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有19人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一道方程组的求解
已经有16人回复
- 用mathematics解方程组
已经有19人回复
- 三元二次的三角方程如何用MATLAB求解？？
已经有10人回复
- 请教大家一个MATLAB的方程求解问题
已经有6人回复
- 一元高次方程怎么解
已经有5人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求助Peng-Robinson状态方程！
已经有5人回复
- 【求解】 一元四次方程
已经有5人回复
- 微分方程当做一元二次方程来解
已经有5人回复
- 含偏微分之绝对值项的一阶偏微分方程组的数值解法
已经有7人回复
- 非线性方程组的解法
已经有5人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- 一元五次方程怎么求解？
已经有8人回复
- Matlab如何找一个一元二次方程式的eigenvalues
已经有3人回复
- 求助！matlab求解一元高次方程的问题！
已经有8人回复
- 如何求解含阻尼的双自由度系统强迫振动的固有频率
已经有8人回复
- 求助解一道高中数学方程题
已经有8人回复
- 请问如何求解二元一阶微分方程组
已经有24人回复
- 【求助】帮忙求解一个一元四次方程
已经有9人回复
- 一元高次方程因式分解
已经有14人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
li2008jie
铁杆木虫 (著名写手)
逍遥孤鹤
- 应助: 10 (幼儿园)
- 金币: 8670.4
- 红花: 7
- 帖子: 1090
- 在线: 531.9小时
- 虫号: 590822
- 注册: 2008-09-01
- 性别: GG
- 专业: 金属功能材料
2楼2013-11-25 18:44:42
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11208.1
- 红花: 91
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
3楼2013-11-25 19:07:13
sddhywj2008
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3027.7
- 散金: 120
- 帖子: 290
- 在线: 156.4小时
- 虫号: 895195
- 注册: 2009-11-06
- 性别: GG
- 专业: 计算机应用技术
4楼2013-11-25 19:45:57
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11208.1
- 红花: 91
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
5楼2013-11-25 20:37:55
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
sddhywj2008: 金币+20, ★★★很有帮助 2013-12-08 10:44:16
感谢参与，应助指数 +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
sddhywj2008
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3027.7
- 散金: 120
- 帖子: 290
- 在线: 156.4小时
- 虫号: 895195
- 注册: 2009-11-06
- 性别: GG
- 专业: 计算机应用技术
7楼2013-11-26 14:19:29
8楼2013-11-27 09:30:01
40