function [L] =legendre(n)
%Antaa n. Legendren polynomin
if n==0
    Le=0;
elseif n==1
    Le=[2 -1];
else
L=zeros(2,n+1);
L(1,end)=1;                 %tehdään L0
L(2,end-1:end)=[2 -1];   %tehdään L1
for k=2:n
    a=conv(L(2,:),L(k,:));
    a=(2*k+1)*a(end-n:end);
    b=-k*L(k-1,:);
    L(k+1,:)=(a+b)/(k+1);
end
end

