diff --git a/run.m b/run.m
index aa53c83..9c4398a 100644
--- a/run.m
+++ b/run.m
@@ -74,12 +74,13 @@ while(k<=kmax && maxD> EPS)
     If0(Balance)=0;
     If2(Balance)=0;
     %Vf0=fsY0\If0;
-    %[fsY0L,fsY0U,fsY0P,fsY0Q,fsY0R]
     Vf0=fsY0Q*(fsY0U\(fsY0L\(fsY0P*(fsY0R\If0))));
     %Vf2=fsY2\If2;
     Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2))));
+    fprintf('迭代时间%f\n',toc);
 end
-fprintf('Fortiscue法计算时间 %f\n',toc);
+FortiscueToc=toc;
+fprintf('Fortiscue法计算时间 %f\n',FortiscueToc);
 Vf1=Vmf1.*exp(1j*Vaf1);
 %%
 (Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3;
@@ -143,7 +144,10 @@ jaco=Jacobi(Balance3P,busNum*3,QGi3P,Vp3m,YdotSinVolt,YdotCosVolt,diag_Volt_Ydot
 [dV, dVangle]=Solv(busNum*3,jaco,dP,dQ);%解出修正量
 [Vp3m, Vp3a]=Modify(Vp3m,Vp3a,dV,dVangle,1);
 fprintf('第 %d 次迭代， 最大不平衡量为 %f\n',k,full(maxD));
+fprintf('迭代时间%f\n',toc);
 end
-fprintf('牛顿法计算时间 %f\n',toc);
+NewtonToc=toc;
+fprintf('牛顿法计算时间 %f\n',NewtonToc);
+fprintf('加速比为%f\n',NewtonToc/FortiscueToc);