##############################################################
# Rufen Sie dieses Programm unter Maple mit der Befehlszeile #
#							     #
#	read aufgabe20;					     #
#							     #
# auf. Semikolon nicht vergessen!			     #
##############################################################

# Splinebibliothek laden
readlib(bspline);

# Voraussetzen, dass h positiv ist
assume(h,positive);

# B-Spline b_j definieren
b[j] := bspline(1,'t',[j*h,(j+1)*h,(j+2)*h]);

# Funktion g definieren
g := t^2;

# Koeffizient berechnen
f[j] := 1/h*int(b[j]*g,t=-infinity..infinity) + 
      h/6*int(diff(b[j],t)*diff(g,t),t=-infinity..infinity);

# Ergebnis schoener darstellen
f[j] := factor(f[j]);

# Fehler an der Stelle tau_{j+1}
delta[j] := factor(f[j] - ((j+1)*h)^2);