24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zouji1983

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 567.4
散金: 2
红花: 3
帖子: 323
在线: 550.6小时
虫号: 445773
注册: 2007-10-28
性别: GG
专业: 材料科学

[求助] 求助一个小程序！

我是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。非常感谢各位的帮组！

回复此楼

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助一个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人回复

1楼 2012-12-04 17:17:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zouji1983

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 567.4
散金: 2
红花: 3
帖子: 323
在线: 550.6小时
虫号: 445773
注册: 2007-10-28
性别: GG
专业: 材料科学

引用回帖:

2楼: Originally posted by chembetsey at 2012-12-04 23:00:54
call subroutine的时候把数组元素值传进去就是了，call calculation (x(i,j,k),A(i,j,k))
或者先把数组元素值赋给一个临时变量，再传递进去

我按照你说的做了，编译都过了，但是运行的时候，显示：log domian error. 不知道怎么解决呢？谢谢！