for n=1:10
    s=0 ;
    for k=1:n
        s=s+k*(n-k)
    end
    S(n)=s/n
    afficher([n,S(n)])
end
clf
plot(S,"+r")