Matematické Fórum

Zdravím,
mám problém s touto funkcí:

p (x,t) = \sum_{k=0}^{inf} [L*(L*cos(k*pi)*sin(L)-k*pi*cos(L)*sin(k*pi)))/(L^2-k^2*pi^2) * cos(k*pi*c/L*t) * cos(k*pi/L*x]

kde:
      L = 1;
      c = 341;

Mělo by se jednat o rovnici změny akustického tlaku mezi dvěmi stěnami o vzdálenosti L.
Potřeboval bych tuto funkci vypočítat a zobrazit ji v Matlabu.
Měl jsem tento nápad (viz níže), ale je to výpočtově náročné a výsledek se mi zdá pochybný.

L = 1; 
c = 340;
t = 0:0.1:2; 
x = 0:0.1:L;
syms k;

A = ones(length(x),length(t));

for i = 1:(length(x))
    disp(i);
    for j = 1:(length(t))
          A(i,j)  = symsum(   (L*(L*cos(k*pi)*sin(L)-k*pi*cos(L)*sin(k*pi)))/(L^2-k^2*pi^2) * cos(k*pi*c/L*t(j)) * cos(k*pi/L*x(i)), 1, 100    );
    end
end

plot(x, A(:,2))  % vykresleni hodnot akustickeho tlaku třeba v čase t(2) podel osy x

Kdyby měl někdo lepší nápad (jako, že by určitě měl nějaky být), tak budu za každý rád. Díky