% 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(n)
% lee n el numero de funciones propias  
syms t 
f=t   
%
%
hold on
dt=0:0.05:1;
df=subs(f,t,dt);
plot(dt,df,'r') 
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
 