function PlotSolution(model)
  
    x=model.x;
    y=model.y;
    x0=model.x0;
    y0=model.y0;
    
    %%
    for s=1 : model.S
         for k=1: model.I(s)
           if model.r(s,(k))==0
               p(k)='b';
           else
              p(k)='r'; 
           end
              
        plot(x(s,k),y(s,k),'o',...
            'LineWidth',1,...
            'MarkerSize',4,...
            'MarkerEdgeColor',p(k));
        hold on;
          
         end
         
     %%
    axis equal
    plot(x0(s),y0(s),'ks',...
        'LineWidth',1,...
        'MarkerSize',6,...
        'MarkerFaceColor','yellow');
    
   % hold off;
    end
end