h = 0.01;
tau = 0.125;

%ñåãà çàäàâàìå îòêúäå äîêúäå ñå ìåíÿò ïðîìåíëèâèòå è ñúîòâåòíèòå ñòúïêè

x = 0:h:4;      % 0 < x < 4
t = 0:tau:8;    % 0 < t < 8
L = length(x);  % L = 401
M = length(t);  % M = 65 
N = 40;  % ïî óñëîâèå èñêàìå ñàìî ïúðâèòå 40 ÷ëåíà íà ðåäà

% ïðèñâîÿâàíå íà ñòîéíîñòèòå îò Àk(k)
%  for k = 1:L    % öúêúë ïî L = 401

for k = 1:N      % öèêúë ïî N = 40   
    A(k) = Ak(k);
end
    
for n = 1 : M   % M = length(t) 
    % u(x,t) = 1 + x/2 + Sum(Ak*sin(k*pi*x/4)*cos(k*pi*t/4))
    u = 1 + x/2 ; 

%      öèêúë ïî L = 401    
%      for l = 1 : L
%          u = u + A(l)*sin(n*pi*x/4)*cos(n*pi*t(n)/4);
%      end
    
%      öèêúë ïî N =40
     for k = 1 : N    
         u = u + A(k)*sin(n*pi*x/4)*cos(n*pi*t(n)/4);
     end 
     
    plot (x,u)
    axis([0 4 -1 5])
    Z(n)=getframe;
end
movie(Z,10,5);
