24小时热门版块排行榜

>论坛更新日志 (2788)
>虫友互识 (587)
>文献求助 (135)
>导师招生 (125)
>休闲灌水 (123)
>考博 (116)
>教师之家 (71)
>找工作 (68)
>博后之家 (67)
>硕博家园 (65)
>论文投稿 (63)
>论文道贺祈福 (61)
>基金申请 (57)
>考研 (50)
>公派出国 (47)
>招聘信息布告栏 (30)

返回列表

woxiaolong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4
帖子: 7
在线: 1.7小时
虫号: 2526560
注册: 2013-06-30
专业: 金属材料的凝固与结晶学

[求助] fortran程序总是出现“root must be bracketed for zbrent”求指点

-13.34=2logA+3logB+2(-6.6B+0.045A-0.001A2)+3(-0.2B-3.9A+1.7A2)
这个是我想解的方程，想输入一个A就能对应得到一个B
A取值范围是10e-5-5e-2，B取值范围是5e-6--1e-4
希望高手能帮忙解答一下，小弟感激不尽，下面是程序

PROGRAM main
   INCLUDE 'DEFIN.DIM'
EXTERNAL func

x1=1.D-5
x2=5.D-2
tol=1.D-10
999 BB=zbrent(x1,x2,tol)
GOTO 999
WRITE(*,*)BB
END

   FUNCTION zbrent(x1,x2,tol)
   INCLUDE 'DEFIN.DIM'

c    INTEGER ITMAX
c    REAL zbrent,tol,x1,x2,func,EPS
EXTERNAL func
   PARAMETER (ITMAX=100,EPS=3.e-8)
c    INTEGER iter
c    REAL a,b,c,d,e,fa,fb,fc,p,q,r,s,tol1,xm
   a=x1
   b=x2
   fa=func(a)
   fb=func(b)
   if((fa.gt.0..and.fb.gt.0.).or.(fa.lt.0..and.fb.lt.0.))pause
   *'root must be bracketed for zbrent'
   c=b
   fc=fb
   do 11 iter=1,ITMAX
      if((fb.gt.0..and.fc.gt.0.).or.(fb.lt.0..and.fc.lt.0.))then
      c=a
      fc=fa
      d=b-a
      e=d
      endif
      if(abs(fc).lt.abs(fb)) then
      a=b
      b=c
      c=a
      fa=fb
      fb=fc
      fc=fa
      endif
      tol1=2.*EPS*abs(b)+0.5*tol
      xm=.5*(c-b)
      if(abs(xm).le.tol1 .or. fb.eq.0.)then
      zbrent=b
goto 99
c       return
      endif
      if(abs(e).ge.tol1 .and. abs(fa).gt.abs(fb)) then
      s=fb/fa
      if(a.eq.c) then
         p=2.*xm*s
         q=1.-s
      else
         q=fa/fc
         r=fb/fc
         p=s*(2.*xm*q*(q-r)-(b-a)*(r-1.))
         q=(q-1.)*(r-1.)*(s-1.)
      endif
      if(p.gt.0.) q=-q
      p=abs(p)
      if(2.*p .lt. min(3.*xm*q-abs(tol1*q),abs(e*q))) then
         e=d
         d=p/q
      else
         d=xm
         e=d
      endif
      else
      d=xm
      e=d
      endif
      a=b
      fa=fb
      if(abs(d) .gt. tol1) then
      b=b+d
      else
      b=b+sign(tol1,xm)
      endif
      fb=func(b)
11 continue
   pause 'zbrent exceeding maximum iterations'
   zbrent=b
99 return
   end

FUNCTION func(BB)
INCLUDE 'DEFIN.DIM'
AA=1.D-5
   func=2*log(AA)+3*log(BB)+2*(-6.6*BB+0.045*AA-0.001*AA**2)
   +    +3*(-0.2*BB-3.9*AA+1.7*AA**2)+13.34

c END IF
RETURN
END

回复此楼