
// Question 2
X=[-5:0.01:20]
R=[]

for n=0:5
    R=[R;laguerre(n,X)]
end

scf(0)
clf()

plot2d(X',R',rect=[-5,-10,20,20])

// Question 4
X=[-1:0.01:1]
R=[]

for n=1:5
    L0=laguerre(n-1,X)
    L1=laguerre(n,X)
    LP=derive(L1,-1,0.01)
    S=X.*LP-n*L1+n*L0
    R=[R;S]
end

disp(R)
// Il n'y a que des nombres très proche de 0 dans la matrice R donc cela "démontre" numériquement l'égalité.
