% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
function y= seriefourier_disc(n)
syms t 
f=t*(heaviside(t)-heaviside(t-1/2))
%
dt=0:0.01:1;
df=subs(f,t,dt);
plot(dt,df,'*r') 
 %
hold on
%
 u=0;
for i=1:n
    c(i)=double(int(sin(i*pi*t)*f,0,1)/int((sin(pi*i*t))^2,0,1));  
   u=u+c(i)*sin(pi*i*t); 
end
%
aproximacion=vpa(u,4) 
%
 ezplot(u,[0,1])
 hold off
 