Classic Fourier Series

[Snorre]

Hi,

another approach:

Code:

#cas
fourier(args):=
BEGIN
  LOCAL argv,argc,f,k;
  argv:=[args];
  argc:=SIZE(argv);
  f:=argv(1);
  k:=argv(argc);
  IF argc=3 THEN f:=zip('unapply',f,argv(2)); END;
  RETURN int(f(t)*e^(-i*k*t),t,-π,π)/π;  // <-- replace that by your a0/ak/bk formulas
END;
#end

This doesn't check types, but number of args, so usage is either fourier(expr,var,k) or fourier(func,k).