l=101
x=rand(1:l);
g=1:1:length(x);
y=pi.*x-.5*pi.*g;
z=rand(1,l);
zz=-log(z);
alfa=1/2;
%plot(g,zz)
for i=1:l
v(i)=1/(cos(y(i))^(1/alfa));   
ya(i)=(sin(alfa*y(i)))*v(i);
m(i)=(cos((1-alfa)*y(i))/zz(i))^((1-alfa)/alfa);
p(i)=ya(i)*m(i);
end
plot(g,p)