function po=lege01(m,N)    
if m==0
    po=1;
elseif m==1
    po=[2,-1];
else
    v0=[0,1];
    v1=[2,-1];
    for n=1:(m-1) 
        ap0=[0,v0];
        ap1=2*[v1,0]-[0,v1];
        po=((2*n+1)*ap1-n*ap0)/(n+1);
        v0=[0,v1];
        v1=po;
    end
end
if N>m+1
    ap=zeros(1,N-m-1);
    po=[ap,po];
end   


