x0=3;
y0=3;
r=2;
n=100;
h=(2*r)/(n-1);
X=x0-r:h:x0+r;
Y1=sqrt(r^2-(X-x0).^2)+y0;
Y2=-sqrt(r^2-(X-x0).^2)+y0
plot(X,Y1,'r-');
hold on
plot (X,Y2,'r-');
axis equal;


alfa=0:2*pi/100:2*pi;
x=r*cos(alfa)+x0;
y=r*sin(alfa)+y0;
plot(x,y,'b-');




ile=0
s=0;
while s<=1000 % ~(s>1000)
s=s+round(rand(1)*10)
ile=ile+1;
end
s
ile


n=100;
s=0;
for i=1:2:n
s=s+1;
end

function [z1,z2]=fun1(x,y,z)
z1=x+y+z;
z2=x*y*z;
s
sum(1:2:n)


A=rand(10,6)*-5;
ile=0;
for i=1:10
for j=1:6
if A(i,j)>0
ile=ile+1
end
end
end
ile

ile=0
s=0;
while s<=1000 % ~(s>1000)
s=s+round(rand(1)*10)
ile=ile+1;
end
s
ile