f:read("Gib die Funktion f(x):= ein");

x0:read("Gib die x-Koordinate eines Punktes ein");

y0:f(x0);

ab:diff(f(x),x);

k:ab,x=x0;

g(x,y):=y=k*x+d;

g:g(x0,y0);

l:solve(g,d);

abschnitt:d,l[1];

steigung:k;

y(x):=steigung*x+abschnitt;

print("Tangente y =",y(x));

wxplot2d([f(x),y(x)], [x,x0-5,x0+5], [gnuplot_preamble, "set size ratio 1; set zeroaxis;"])$