function fouriermain
clear all;
clc;

f1=@(t) cos(t);
f2=@(t) sin(t);
f3=@(t) sin(t).*cos(t);

n=5;

[a0 a b]=fouriercoeff(f3,n);

disp('a_0=')
disp(a0)
disp('a_n')
disp(a)
disp('b_n')
disp(b)
end

function [a0 a b] = fouriercoeff(f,n,T)
if nargin==2
    T=2*pi; % normalerweise 2Pi-periodisch
end

w=2*pi/T;

a0=2/T*quad(f,0,T);
ak=@(k) (2/T*quad(@(t) (f(t).*cos(k*w*t)),0,T));
bk=@(k) (2/T*quad(@(t) (f(t).*sin(k*w*t)),0,T));

a=zeros(1,n);
b=a;
for k=1:n
    a(k)=ak(k);
    b(k)=bk(k);
end

end