%%%% Wykonac wykres kwantylowy (q-q plot) dla
%%%%%% Pareto(1.2)+1 - linia 1: oooo
%%%%%% lognormal(0,1) - linia 3: ++++
%%%%%% Pareto(0.5)+1 - linia 3: ****
n =zeros;
u = zeros;
s = zeros;
m = zeros;
x = zeros;
y1 = zeros;
y2 = zeros;
y3 = zeros;
y4 = zeros;

% wspolrzedna x-owa
for i=1:70
x (i) = -log(1-(i/71));
endfor

%Pareto(1.2)
for i=1:70
u (i) = log( ((rand)^(-1/1.2)));
endfor

y1 = sort (u);


%%% logN(0,1)
for i=1:70
n (i) = randn;
endfor

y3 = sort (n);


%% Pareto(0.5)
for i=1:70
m (i) = log( ((rand)^(-1/0.5)));
endfor

y4 = sort (m);


plot (x, y1, '*', x, y3, 'o', x, y4 ,'@')












