如何在Matlab中计算总和

Question

如何使用for循环计算包含求和的函数的解析函数。

我创建了3个for循环，一个为m，一个为t，一个为x，但肯定无法正常工作。它给了我inf而不是一个数字数组，以便我可以对其进行绘图。原始公式是：

L=1;
dx=L/N;
N=10;
t_final=0.1;
 dt=0.01;
 D=1;
 sigma=D*dt/(dx^2);
 x=linspace(0,1,N);
 I=eye(N-2);
 B=zeros(N-2,N-2);
 t=0:dt:t_final

for k=1:100
    for j=1:length(t)
        for p=1:length(x)
            C_exact=C_exact+4*((1/((2*k-1)*pi))*exp(-(2*k-1)^2*pi^2*t(j))*sin((2*k-1)*pi*x(p)));
        end
    end
end