%genereaza variabile stabile simetrice( Cauchy (z),si S(a,0)), unde a este exponentul caracteristic
%w este variabila aleatoare exponentiala

n=20001;a=.5;
e=rand(1,n);

r=ones(size(e));
for i=1:n
r(i)=(pi/2)*(2*e(i)-1);
end
f=randn(1,n);
for j=1:n
%w(j)=.5*(f(j)/(sin(r(j))))^2;
w(j)=-log(e(j));
end
t=-100:2:100;
%hist(w,t)
for k=1:n
z(k)=tan(r(k));
y(k)=(sin(a*r(k))/(cos(r(k))^(1/a)))*(cos((1-a)*r(k))/w(k))^((1-a)/a);
end
hist(y,t)
%pentru vizualizare dilat vectorul y sau z
%v=4.*y;hist(v,t)

