[求助]含集总元件的FDTD仿真

我打算在微带线的一端加入二极管二极管的电场方程是个超越方程，我单独用牛顿法写了个方程，调用后结果一直不对，不知道为什么求指导啊！ qu;$I'Ul%  
     module diode nL]eGC  
     USE common_data_1 HNLr} Yj  
    implicit none >Udq{<]#r  
w8`B}Dr23  
      integer :: i,j,k RW<4",  
    REAL,PARAMETER :: zero=1.e-5 CF:!  
    REAL,PARAMETER :: eps=13.,q=1.6e-19,kerw=1.38e-23,temperature=300. ({cWb:+r  
    contains -U|Z9sia  
    !割线法的参数 w,OPM}) il  
    real function newton(a,f,df) bicbCC6kC  
    implicit none TQQh:y  
i*T -9IP  
    real,intent(in) :: a    !起始的猜想 FRsp?i K)  
    real, external :: f         !输入的求值函数 <00=bZzX  
    real, external :: df            !f'(x)的函数 fk\]wFj  
    real :: b=0.                 !逼近得到的解 Rt@O@oDI  
    real :: fb ~^fb`f+%  
    real :: temp1 equi26jhr  
            temp1=0.0; #S[Y}-]T  
temp1=a   27}0
  
   b=temp1-f(temp1)/df(temp1) N7_(,Gu*R  
     fb=f(b) G\~?.s|^  
   do while(abs(fb) > zero) !iK{q0  
     temp1=b .dwbJT  
     b=temp1-f(temp1)/df(temp1)   dt/-0~U  
     fb=f(b) jI9#OEH_g  
    enddo t~":'le`zr  
     newton=b XQ8q)B=  
    return ``>z8t[ks  
    end function newton & s:\tL  
     _\{/#J;lN  
    real function func(x) f6{.Uq%SGp  
    implicit none KWM}VZY:Z  
    real :: x }%D${.R]  
func=x-e(i,j,k,1)-(dt/(e0*eps))*((h(i,j+1,k,3)-h(i,j,k,3))/dy-(h(i,j,k+1,2)-h(i,j,k,2))/dz)+(canshu*dt/(e0*eps*dy*dz))*(exp(q*dx*(x+e(i,j,k,1))/(2*kerw*Temperature))-1.) Q
m_;o(  
return ^4et; F%  
    end function func \~l"  
                 q(v|@l|)yO  
    real function dfunc(x) CV2#G*  
    implicit none v<mSd2B*  
 &n .. moZeP#Q%  
8EBd`kiq  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册