function y=f(x);y=x^2*sqrt(x)
endfunction
n=input("n = ")
u=0;v=0;x=1
for k=0:(n-1)
    u=u+1/n*f(x)
    x=x+1/n
    v=v+1/n*f(x)
end
afficher([u,v])