Hi Sayandev,
As you know, iterative methods need an initial guess. But this is the
user's *choice*: you can use a flat start as you mention, or any other
information you may have (for instance, from previously solved cases that
are similar to the one at hand---as it happens when solving N-1 conti
Beg your pardon the missing snippet from the previous mail is as follows:
%% initial state
%V0= ones(size(bus, 1), 1);%% flat start
V0 = bus(:, VM) .* exp(sqrt(-1) * pi/180 * bus(:, VA))
vcb = ones(size(V0)); %% create mask of voltage-controlled
buses