24小时热门版块排行榜    

查看: 810  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

利维坦6021

新虫 (初入文坛)

[求助] 【求助】使用fortran95牛顿法求根

最近被一个超越方程困扰了很久,因为对fortran不太熟悉,因此求教各位如何使用fortran语言求这个方程的根呢?
公示在图中,其中p为常数。
在网上搜到的程序包有点看不明白,希望各位不吝赐教

【求助】使用fortran95牛顿法求根
捕获.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylsxz2012

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by lijf517 at 2017-09-16 06:25:12
这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。
...

的确是最简单的。一般微积分书上都有。给您更正一个小小的笔误。
x_(k+1)=x_k-f‘_k*dx,
这里dx是步长。
4楼2017-09-16 08:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

scoot

新虫 (初入文坛)

2楼2017-09-15 18:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijf517

木虫 (小有名气)

这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。

发自小木虫Android客户端
3楼2017-09-16 06:25:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylsxz2012

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by lijf517 at 2017-09-16 06:25:12
这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。
...

我刚才应该是dt不是dx,你用x做因变量了,没有仔细看。呵呵。
x_(k+1)=x_k-x‘_k*dt
dt步长。
5楼2017-09-16 09:01:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见