%begin iteration


for jj=1:100
   
if e==1
   get_start;
   newdata;
   itr=itr+1;
end

if e==0
%    disp('iteration times is: itr=');
%    disp(itr);
   itr=1;
%    disp('Pwi=');
%    pwi
%    disp('Swi=');
%    swi
   
   
   
   pwit=pwi;
   swit=swi;
   phit=phi;
   
   
   break;
end

end