function[imim1tmtm1tftt]=cjdlfsq_1（RLC）??????%P262沖擊電流發生器FORTRAN改編程序
dtf=0.05;
dtt=0.1;
u=122.5;
nm=280;
i=zeros（999）;
t=zeros（999）;
for?k=1:nm
????t（k+nm）=k*dtt;
????t（k）=k*dtf;
end
for?k=1:nm
????t（k+2*nm）=（nm+k）*dtt;
end
af=R/（2.0*L）;
w0=1.0/sqrt（L*C）;
afd=sqrt（af*af-w0*w0）;
w=sqrt（w0*w0-af*af）;
if?（（af-w0）<0）
????for?j=1:3*nm
????????i（j）=u*exp（-af*t（j））*sin（w*t（j））/（w*L）;
????end
????bt=atan（w/af）;
????tm1=bt/w;
????disp（[‘計算峰值時間Tm1為:?‘num2str（tm1）‘?μs‘]）;
????im1=u*exp（-af*bt/w）/sqrt（L/C）;
????disp（[‘計算電流峰值Im1為:?‘num2str（im1）‘?kA‘]）;
elseif?（（af-w0）==0）
????for?j=1:2*nm;
????????i（j）=u*t（j）*exp（-af*t（j））/L;
????end
????tm1=sqrt（L*C）;
????im1=u*sqrt（C/L）*exp（-1.0）;
else
????p1=-af+afd;
????p2=-af-afd;
????for?j=1:2*nm
????????i（j）=u*（exp（p1*t（j））-exp（p2*t（j）））/（L*（p1-p2））;
????end
????tm1=log（p2/p1）/（p1-p2）;
????im1=u*（exp（p1*tm1）-exp（p2*tm1））/（L*（p1-p2））;
end
im=i（1）;
for?j=1:nm
????if（i（j）>im）
????????im=i（j）;
????end
end
disp（[‘實際電流峰值Im為: