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

我打算在微带线的一端加入二极管二极管的电场方程是个超越方程，我单独用牛顿法写了个方程，调用后结果一直不对，不知道为什么求指导啊！ 'jv[Gcss3L  
     module diode UL0%oJ#  
     USE common_data_1 5B+>28G%  
    implicit none } $:uN  
tz(\|0WDQ  
      integer :: i,j,k ,X Zo0!  
    REAL,PARAMETER :: zero=1.e-5 u(ep$>[F#_  
    REAL,PARAMETER :: eps=13.,q=1.6e-19,kerw=1.38e-23,temperature=300. _*b1]<  
    contains
x,
W)qv  
    !割线法的参数 $Iuf(J-5[  
    real function newton(a,f,df) [da,SM  
    implicit none $t' .  
%j],6wW5J  
    real,intent(in) :: a    !起始的猜想 P.bBu  
    real, external :: f         !输入的求值函数 B~HA 32  
    real, external :: df            !f'(x)的函数 #NZ\UmA  
    real :: b=0.                 !逼近得到的解 +dWx?$n  
    real :: fb <{rRcFR  
    real :: temp1 *2O4*Q1  
            temp1=0.0; Pkx*1.uo  
temp1=a   J)O1)fR  
   b=temp1-f(temp1)/df(temp1) nZ?BCO  
     fb=f(b) e|:#Y^  
   do while(abs(fb) > zero) vywd&7gK  
     temp1=b F,0@z/8a  
     b=temp1-f(temp1)/df(temp1)   O[O`4de9  
     fb=f(b) w3lR8R]  
    enddo ZKKz?reM'  
     newton=b AJ-p|[wPz  
    return }V.Wp6"S   
    end function newton 4&r+K`C0  
     HJpkR<h  
    real function func(x) Kp") %p#  
    implicit none s^KUe%am0  
    real :: x BXgAohg!  
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.) [;M31b3  
return x2B~1
edf  
    end function func V$u~}]z  
                 uLV@D r   
    real function dfunc(x) {~DYf*RZ  
    implicit none fu|N{$h%X  
 &n .. h4CDZ  
2XJn3wPi  

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