f:read("Gib den Funktionsterm f(x) ein");

ab:diff(f,x);

ab2:diff(f,x,2);

K:ab2/(1+ab**2)**(3/2);

K(x):=''K;

x0:read("An welcher Stelle soll der Krümmungskreis bestimmt werden");

Kruemmung:K(x0);

Kruemmungsradius:1/abs(Kruemmung);