| 查看: 738 | 回复: 9 | ||
[求助]
求助一个小程序！
|
|
我是fortran的初学者，现在用别人的程序。但我现在需要用newton-raphson解一个非线性方程，例如：lnx+x^3+A=0, 这里我要解的未知数x和常量A都是一个数组，像x是这样一个数组：(1,1,100,100), A是（1,100,100）. 我只能从网上找到关于解x不是数组的code,如下： subroutine calculation (x,A) integer ::k real(8) :: x,x0,A real(8):: tol f(x)=lnx+x**3+A df(x)=1/x+3*x**2 tol=1e-6 x0=0.2d0 k=0 10 k=k+1 x=x0-f(x0)/df(x0) if(abs(x-x0).lt.tol)goto 20 x0= x goto 10 20 write(*,*) x end do end do end subroutine calculation 但是对于数组来说应该怎么改呢？我想把这个作为一个subroutine。非常感谢各位的帮组！ |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助一个matlab关于分类讨论的问题
已经有9人回复
- 求助一个编程问题！
已经有10人回复
- c语言求助，一个小程序详细讲解
已经有13人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 【求助】如何让一个程序文件不能被复制
已经有4人回复
- 求助高手帮忙编一个Fortran的小程序，重金悬赏......
已经有15人回复
- 【求助】求助帖：请教大家一个程序
已经有10人回复
- 【求助】一个求数组绝对值的小问题
已经有7人回复
- 【求助】请教一个关于Mathematica中计算小量的问题
已经有4人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】求助帮助编一个读取TEXT文件画图的程序
已经有8人回复
- 【求助】急需一个程序
已经有8人回复
- 【求助】运行ms的时候跳出“由于另一个程序正在运行中，此操作无法完成。请选择切换…
已经有28人回复
- 【求助】新手求助。。。。。。。。。。。。
已经有6人回复
- 【求助】一个小程序--出现错误
已经有3人回复
- 【求助】帮忙看一段程序
已经有9人回复
- 【求助】100金币求助高人一个程序解方程
已经有15人回复
- 【求助】求翻译一个句子【已完成】
已经有4人回复
- 【求助】帮忙编个用于VASP小程序
已经有22人回复
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
2楼2012-12-04 23:00:54
3楼2012-12-04 23:49:32
4楼2012-12-04 23:54:32
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
5楼2012-12-05 01:47:49
6楼2012-12-05 16:54:09
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
7楼2012-12-05 22:12:41
8楼2012-12-06 18:44:33
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
9楼2012-12-06 22:16:46
10楼2012-12-07 23:53:17
