taylorPicture := proc(zentrum, eps, n::posint) local a, b, k, t, tp, q, p; a := zentrum - eps; b := zentrum + eps; for k to n do t := taylor(f(x), x = zentrum, k); tp := unapply(convert(t, polynom), x); q[k] := plot(tp(x), x = a .. b, y = -4 .. 4, numpoints = 1000, color = blue, args[4 .. nargs]) od; q := plots[display]([seq(q[k], k = 1 .. n)], insequence = true); p := plot(f(x), x = a .. b, y = -4 .. 4, color = red, args[4 .. nargs]); plots[display]([q, p]) end

Die Funktion, von der die Taylorreihe entwickelt werden soll, muss zunächst als Funktion f definiert werden
Der Aufruf hat dann die Signatur

taylorPicture(<Entwicklungs-Zentrum>, <Zeichenbereich in Entfernung vom Zentrum>, <Entwicklungsgrad>)

Aufruf:
> f:=x->sin(x);
> taylorPicture(Pi,4.5*Pi,32);