eval(['load ',CurrentSystem,'.mat']);

% make initial values for param from input data
%ii=no_gen+1;
%last_one='Q';
%for i=1:k_params
% if i<=no_gen-1,
%        param(i)=bus_p(i+1);
% end;
% if (i>no_gen-1) & (i<=no_gen+no_pv-1),
%        param(i)=bus_p(ii);
%        ii=ii+1
% end;
% if (i>no_gen+no_pv-1) & (i<=no_gen+no_pv-1+2*no_pq),
%        if last_one=='Q',
%                param(i)=bus_p(ii);
%                last_one='P';
%        else 
%                param(i)=bus_q(ii);
%                last_one='Q';
%                ii=ii+1;
%        end;
% end;
%end;

%p=zeros(k_params,1);

lfsetsta;
lfsetpar;
v=zeros(length(x),1);