# Einfach alle Zeilen nacheinander mit RETURN best\303\244tigen. F\303\274r andere Funktionen:
# unten f:=abs(x) auskommentieren, bei einer anderen das # wegmachen. Oder neues f angeben
with(plots):N:=5: # H\303\266here Werte -> bessere Approximation # f:=piecewise(x<0,-1,x>0,1); # Rechteckfunktion
# f:=x; # S\303\244gezahnfunktion
f:=abs(x); # Dreiecksfunktiona[0]:=int(f,x=-Pi..Pi)/Pi:
for j from 1 to N do:
a[j]:=1/Pi*int(f*cos(j*x),x=-Pi..Pi):
b[j]:=1/Pi*int(f*sin(j*x),x=-Pi..Pi):
od:p1:=plot(f,x=-Pi..Pi,scaling=constrained):p2:=plot(a[0]/2+add(a[j]*cos(j*x)+b[j]*sin(j*x),j=1..N),x=-5...5,colour=blue,scaling=constrained):display([p1,p2]);LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2I1EhRic=LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2I1EhRic=