function [ql,qm,qu]=kvartiilit(v)
% Tätä en uskaltanut testata kuin tap, n=10, jolla toimii oikein
% Sattaa toimia muillakin, mutta en jaksa näperrellä tuota
% inhottavaa näpräystä, jos nyt sattuisi yhden indeksin heitto
% jossain tapauksessa.
n=length(v);
v=sort(v);
iM=(n+1)/2;i1=floor(iM);i2=ceil(iM);
qm=(v(i1)+v(i2))/2;
iL=(i1+1)/2;
iU=i1+iL;
i1=floor(iL);i2=ceil(iL);
ql=(v(i1)+v(i2))/2;
i1=floor(iU);i2=ceil(iU);
qu=(v(i1)+v(i2))/2;
