function chilun(m,z,ha,c)
r=m*z/2;
rb=r*cos(pi/9);
ra=r+ha*m;
rf=r-(ha+c)*m;
sita=2*pi/z;
almax=acos(rb/ra);
alg=acos(rb/r);
ct1=tan(alg)-alg;
ct2=tan(almax)-almax;
bita=sita/4+ct1-ct2;
for i=0:1:z
fai=-bita:0.001:bita;
xa=ra*cos(fai+i*sita);
ya=ra*sin(fai+i*sita);
hold on;
plot(xa,ya,'b');
fail=bita+ct2:0.01:sita-bita-ct2;
xb=rb*cos(fail+i*sita);
yb=rb*sin(fail+i*sita);
hold on;
plot(xb,yb,'r');
for al=almax:-0.001:0;
ct=tan(al)-al;
ri=rb/cos(al);
alfa=-ct+(i+0.25)*sita+ct1;
x2=ri*cos(alfa);
y2=ri*sin(alfa);
plot(x2,y2,'b');
end
for al2=0:0.001:almax
ct=tan(al2)-al2;
ri=rb/cos(al2);
beta=ct+(i-0.25)*sita-ct1;
x3=ri*cos(beta);
y3=ri*sin(beta);
plot(x3,y3,'g');
end
end