% Начални стойности
t = 0:0.02:1;
x = 0:0.01:pi;
N = 35;
clear('A');

% Пресмятаме интегралите и стойностите на Ak
for k = 1:N
   yx = (2*pi*x-x.^2).*(abs(sin((3*x)/2))).*cos(((2*k+1)*x)/2);
   A(k) = 2*trapz(x,yx)/pi;
end

% Анимираме
for m = 1:length(t)
   u = 0;
   for k = 1:N
      u = u+(A(k)*exp(-(2*k+1).^2*t(m)/4)).*cos(((2*k+1)*x)/2);
   end
   plot(x,u);
   axis([0,pi,-10,10])
   M(m) = getframe;
end
