PR=0;GR=0;
nbfois=input("nombre de trapèzes = ");
ecart=(3-1)/nbfois;
for k=1:nbfois 
    a=1+(k-1)*ecart;
    b=1+k*ecart;
    PR=PR+(b-a)*(a+b)^2/4;
    GR=GR+(b-a)*(a^2+b^2)/2;
end;
mes="l''aire est comprise entre "+string(PR)+" et " +string(GR);
disp(mes);