SRiemann(f,a,b,n) := block(
    /* les variables locales                    */
    [s,v,r], 
    /* la suite de points de la subdivision     */
    s:makelist(a+(b-a)/n*i,i,1,n),  
    /* les valeurs de la fonction en ces points */
    v:map(f,s),                     
    /* la somme de Riemann                      */
    r:lsum(i*(b-a)/n,i,v),          
    /* une évaluation comme résultat final      */
    float(r)                        
);