Как вычислить аналитическую функцию функции, которая включает суммирование, используя цикл for.
Я создал 3 цикла для цикла: один для 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