Avec MuPAD, maple ou xcas

On définit par exemple la fonction

f₂(x) = x - $${\frac{x^3}{3!}}$$

en tapant f2:=x->(x-x^3/3!); Pour définir une fonction par morceaux, on peut insérer un test if dans la définition.

On calcule la valeur de l'intégrale I = $\int_{0}^{1}$f₂(x)dx en tapant :
float(int(f2(x),x=0..1)); (MuPAD) ou evalf(int(f2(x),x=0..1)); (maple/xcas)
Si f₂ n'admet pas de primitive continue sur [0;1]:
float(hold(int(f2(x),x=0..1)));/evalf(Int(f2(x),x=0..1));
on effectue ainsi directement l'intégrale par une méthode numérique (hold empêche l'évaluation en mupad, et Int est la forme inerte (non évaluée) de l'intégrale en maple, on évite ainsi la recherche d'une primitive de f₂).

Pour tracer le graphe de fonction(s), on utilise la commande plotfunc2d en mupad et plot en maple/xcas Par exemple pour avoir le tracé simultané de y = sin(x) et de y = x, on tape:
plotfunc2d(sin(x),x,x=0..4);/plot([sin(x),x],x=0..4);