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

我打算在微带线的一端加入二极管二极管的电场方程是个超越方程，我单独用牛顿法写了个方程，调用后结果一直不对，不知道为什么求指导啊！ o7t{?|  
     module diode Lw2VdFi>E&  
     USE common_data_1 |]?zH~L  
    implicit none 0` .5gxm  
Ss%1{s~ok  
      integer :: i,j,k |kP utB  
    REAL,PARAMETER :: zero=1.e-5 mRa\ wEg%  
    REAL,PARAMETER :: eps=13.,q=1.6e-19,kerw=1.38e-23,temperature=300. YS_3Cq  
    contains G[1\5dK*uR  
    !割线法的参数 n*
Uk<_WA  
    real function newton(a,f,df) 3{fg3?  
    implicit none ;tSAQ  
j+@3.^vK
 
    real,intent(in) :: a    !起始的猜想 c_8<N7 C  
    real, external :: f         !输入的求值函数 [dAQrou6P  
    real, external :: df            !f'(x)的函数 =r*Ykd;W|E  
    real :: b=0.                 !逼近得到的解 !I.}[9N  
    real :: fb H8B$#.  
    real :: temp1 z;!"i~fFK  
            temp1=0.0; K_@?Q@#YhR  
temp1=a   G;$;$gM  
   b=temp1-f(temp1)/df(temp1) $xWUzg1<U  
     fb=f(b) n}I?.r@e  
   do while(abs(fb) > zero) BM'!odR
v  
     temp1=b Q;J( 5;  
     b=temp1-f(temp1)/df(temp1)   K{{_qFj@<y  
     fb=f(b) X8">DR&>Y  
    enddo vxHFNGI  
     newton=b J>5rkR@/  
    return =J4|"z:  
    end function newton xJ2I@*DN  
     4hODpIF  
    real function func(x) :eSsqt9]9  
    implicit none 0Ni{UV? k  
    real :: x ] |nW  
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.) GX0zir
z  
return Qn,6s%n  
    end function func 3H`{ A/r  
                 -rn6ZSD)  
    real function dfunc(x) Q2D!Agq=D  
    implicit none X fqhD&g  
 &n .. ?*%_:fB  
!3Pbu=(cte  

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