errorexamining.m

前推潮流算法

% error examining
E=0.0001
g=11;
for c=1:g
    dU(c)=abs(U(c)-U1(c));
    if dU(c)>E
        computing;
    else 
        show;
    end
end