курсачи / курсач 3 вар / finckurs

clear all, close all,
x0 = 10;
y0 = 10;
r = 5;%По условию окружность меньше ячейки
w=10;
h=10;
xx0=2;

peresy=0;
dy=0.01;
yy0=1;%координаты угла ячейки

x = x0-r:0.01:x0+r;
y=yy0;
ds=0;
y1 = sqrt(r^2 - (x - x0).^2)+y0;
y2 = -sqrt(r^2 - (x - x0).^2)+y0;
plot(x,y1,x,y2); axis square
rectangle('Position',[xx0 yy0 w h]);
dx=0.01;
sum=0;
x=xx0;
f = @mykavadrat;
mykavadrat(xx0,yy0);
function summ = mykavadrat(xx0,yy0)
    while(x<xx0+w)
    if(x>xx0)&&(x<xx0+w)
        if(x>x0-r)&&(x<x0+r)
            y1 = sqrt(r^2 - (x - x0)^2)+y0;
            if(y1<yy0)
                while (y1<yy0)&&(x<x0+r)
                    y1 = sqrt(r^2 - (x - x0)^2)+y0;
                    x=x+dx;
                end
            end
            y2 = -sqrt(r^2 - (x - x0).^2)+y0;
            if(y1<yy0+h)
                if(y2>yy0)
                    ds=dx*(y1-y2);
                    sum=sum+ds;
                    rectangle('Position',[x y2 dx (y1-y2)]);
                end
                if(y2<yy0)
                    ds=dx*(y1-yy0);
                    sum=sum+ds;
                    rectangle('Position',[x yy0 dx (y1-yy0)]);
                end
            end
            %важно что трубка меньше ячейки иначе нужно рассмотреть еще случий
            if(y1>yy0+h)
                if(y2<yy0)
                    
                    ds=dx*(yy0+h-yy0);
                    if ds>0
                        sum=sum+ds;
                        rectangle('Position',[x yy0 dx (yy0+h-yy0)]);
                    end
                end
                if(y2>yy0)
                    ds=dx*(yy0+h-y2);
                    if ds>0
                        sum=sum+ds;
                        rectangle('Position',[x y2 dx (yy0+h-y2)]);
                    end
                end
            end
            
        end
    end
    x=x+dx;
end
end




% xx1=xx0;
% yy1=yy0+h;
% sum2=0;
% ds=0;
% x=xx1;
% rectangle('Position',[xx1 yy1 w h]);
% axis([0 40 0 40]);
% %!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% while(x<xx1+w)
%     if(x>xx1)&&(x<xx1+w)
%         if(x>x0-r)&&(x<x0+r)
%             y1 = sqrt(r^2 - (x - x0)^2)+y0;
%             if(y1<yy1)
%                 while (y1<yy1)&&(x<x0+r)
%                     y1 = sqrt(r^2 - (x - x0)^2)+y0;
%                     x=x+dx;
%                 end
%             end
%             y2 = -sqrt(r^2 - (x - x0).^2)+y0;
%             if(y1<yy1+h)
%                 if(y2>yy1)
%                     
%                     ds=dx*(y1-y2);
%                     if ds>0
%                         sum2=sum2+ds;
%                         pos = [x y2 dx (y1-y2)];
%                         rectangle('Position',pos,'EdgeColor','r');
%                     end
%                 end
%                 if(y2<yy1)
%                     ds=dx*(y1-yy1);
%                     if ds>0
%                         sum2=sum2+ds;
%                         pos = [x yy1 dx (y1-yy1)];
%                         rectangle('Position',pos,'EdgeColor','r');
%                     end
%                 end
%             end
%             %важно что трубка меньше ячейки иначе нужно рассмотреть еще случий
%              if(y1>yy1+h)
%                 if(y2<yy0)
%                     
%                     ds=dx*(yy1+h-yy0);
%                     if ds>0
%                         sum=sum+ds;
%                         rectangle('Position',[x yy0 dx (yy1+h-yy1)]);
%                     end
%                 end
%                 if(y2>yy0)
%                     ds=dx*(yy1+h-y2);
%                     if ds>0
%                         sum=sum+ds;
%                         rectangle('Position',[x y2 dx (yy1+h-y2)]);
%                     end
%                 end
%             end
%         end
%     end
%     x=x+dx;
% end
% %!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% xx1=xx0+w;
% yy1=yy0;
% sum3=0;
% ds=0;
% x=xx1;
% rectangle('Position',[xx1 yy1 w h]);
% axis([0 40 0 40]);
% while(x<xx1+w)
%     if(x>xx1)&&(x<xx1+w)
%         if(x>x0-r)&&(x<x0+r)
%             y1 = sqrt(r^2 - (x - x0)^2)+y0;
%             if(y1<yy1)
%                 while (y1<yy1)&&(x<x0+r)
%                     y1 = sqrt(r^2 - (x - x0)^2)+y0;
%                     x=x+dx;
%                 end
%             end
%             y2 = -sqrt(r^2 - (x - x0).^2)+y0;
%             if(y1<yy1+h)
%                 if(y2>yy1)
%                     
%                     ds=dx*(y1-y2);
%                     if ds>0
%                         sum3=sum3+ds;
%                         pos = [x y2 dx (y1-y2)];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%                 if(y2<yy1)
%                     ds=dx*(y1-yy1);
%                     if ds>0
%                         sum3=sum3+ds;
%                         pos = [x yy1 dx (y1-yy1)];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%             end
%             %важно что трубка меньше ячейки иначе нужно рассмотреть еще случий
%             if(y1>yy0+h)
%                 if(y2<yy0)
%                     
%                     ds=dx*(yy0+h-yy0);
%                     if ds>0
%                         sum=sum+ds;
%                         pos = [x yy0 dx yy0+h-yy0];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%                 if(y2>yy0)
%                     ds=dx*(yy0+h-y2);
%                     if ds>0
%                         sum=sum+ds;
%                         pos = [x y2 dx yy0+h-y2];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%             end
%         end
%     end
%     x=x+dx;
% end
% %!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% xx1=xx0+w;
% yy1=yy0+h;
% sum4=0;
% ds=0;
% x=xx1;
% rectangle('Position',[xx1 yy1 w h]);
% axis([0 40 0 40]);
% while(x<xx1+w)
%     if(x>xx1)&&(x<xx1+w)
%         if(x>x0-r)&&(x<x0+r)
%             y1 = sqrt(r^2 - (x - x0)^2)+y0;
%             if(y1<yy1)
%                 while (y1<yy1)&&(x<x0+r)
%                     y1 = sqrt(r^2 - (x - x0)^2)+y0;
%                     x=x+dx;
%                 end
%             end
%             y2 = -sqrt(r^2 - (x - x0).^2)+y0;
%             if(y1<yy1+h)
%                 if(y2>yy1)
%                     
%                     ds=dx*(y1-y2);
%                     if ds>0
%                         sum3=sum3+ds;
%                         pos = [x y2 dx (y1-y2)];
%                         rectangle('Position',pos,'EdgeColor','g');
%                     end
%                 end
%                 if(y2<yy1)
%                     ds=dx*(y1-yy1);
%                     if ds>0
%                         sum3=sum3+ds;
%                         pos = [x yy1 dx (y1-yy1)];
%                         rectangle('Position',pos,'EdgeColor','g');
%                     end
%                 end
%             end
%             %важно что трубка меньше ячейки иначе нужно рассмотреть еще случий
%              if(y1>yy0+h)
%                 if(y2<yy0)
%                     
%                     ds=dx*(yy0+h-yy0);
%                     if ds>0
%                         sum=sum+ds;
%                         pos = [x yy0 dx yy0+h-yy0];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%                 if(y2>yy0)
%                     ds=dx*(yy0+h-y2);
%                     if ds>0
%                         sum=sum+ds;
%                         pos = [x y2 dx yy0+h-y2];
%                         rectangle('Position',pos,'EdgeColor','b');
%                     end
%                 end
%             end
%         end
%     end
%     x=x+dx;
% end

% while(x<xx1+w)
%     if(x>xx1)&&(x<xx1+w)
%         if(x>x0-r)&&(x<x0+r)
% y1 = sqrt(r^2 - (x - x0)^2)+y0;
% y2 = -sqrt(r^2 - (x - x0).^2)+y0;
% if(y1<yy1+h)
% if(y2>yy1)
% ds=dx*(y1-y2);
% sum=sum+ds;
% rectangle('Position',[x y2 dx (y1-y2)]);
% end
% if(y2<yy1)
% ds=dx*(y1-yy1);
% sum=sum+ds;
% rectangle('Position',[x yy0 dx (y1-yy1)]);
% end
% end
% if(y1>yy1+h)
%
%     ds=dx*(yy1+h-y2);
%     sum=sum+ds;
%     rectangle('Position',[x y2 dx (yy1+h-y2)]);
% end
%
%         end
%     end
%     x=x+dx;
% end






%
%  peres1x=0;%сделать массив
% peres1y=0;
% peres2x=0;
% peres2y=0;
% peres34x=zeros(1,2);
% peres34y=zeros(1,2);
% peres56x=zeros(1,2);
% peres56y=zeros(1,2);
% pixel=zeros(2,2);

%  yperes=sqrt(r^2 - (x - x0).^2)+y0;
% i=1;
% while(y<yy0+h)
%         if(yperes-y<dy)
%         peres1x=xx0;
%         peres1y=yperes;
%         end
%         y=y+0.01;
% end
% yperes=-sqrt(r^2 - (x - x0).^2)+y0;
% y=yy0;
% while(y<yy0+h)
%         if(yperes-y<dy)
%         peres2x=xx0;
%         peres2y=yperes;
%         end
%         y=y+0.01;
% end
%
% y=yy0+h;
% i=1;
% x=xx0;
% while(x<xx0+w)
%
%     yperes=sqrt(r^2 - (x - x0).^2)+y0;
%         if(y-yperes<dy)
%        peres34x(i)=x;
%        peres34y(i)=yperes;
%        i=2;
%         end
%        if(peres34x(2)- peres34x(1)<dy)
%        peres34x(2)=0;
%        peres34y(2)=0;
%
%         end
%         x=x+0.01;
% end
% x=xx0;
% y=yy0;
% while(x<xx0+w)
%
%     yperes=-sqrt(r^2 - (x - x0).^2)+y0;
%         if(y-yperes<dy)
%        peres56x(i)=x;
%        peres56y(i)=yperes;
%        i=2;
%         end
%        if(peres34x(2)- peres34x(1)<dy)
%        peres56x(2)=0;
%        peres56y(2)=0;
%
%         end
%         x=x+0.01;
% end