Can somebody please help me with this Fourier series File?

Can somebody please help me on how to write this code for this file?

for the first part i have

fig = figure;

fig.WindowState = 'maximized';

function answer = Fourier(f,degree)

syms x;

answer = (1/(2.*pi)).*int(f,-pi,pi); % calculates a_0, the constant term

for n = 1:degree

coefficientOfSin = ((1/pi).*(int(f.*sin(n.*x),-pi,pi)));

coefficientOfCos = ((1/pi).*(int(f.*cos(n.*x),-pi,pi)));

answer = answer+(coefficientOfCos.*(cos(n.*x))+coefficientOfSin.*(sin(n.*x)));

end

end

the second part I have this, but i dont believe its right.

hold on;

syms x;

f=x.^2

fplot((x.^2)+1,[-1,1]);

for counter = 1:4

Tcounter = Fourier(f,counter);

fplot(Tcounter,[-1,1]);

end

axis equal;

axis off;

axis tight;