Hey Deck
The provided code is a finite difference method (explicit forward-time, centered-space) for solving the one-dimensional heat equation. As per my understanding you're asking how to implement the separation of variables method instead, while keeping the same values for parameters.
The method of separation of variables works by assuming that the solution can be written as a product of two functions, one dependent only on x and the other only on t:
Below is a MATLAB implementation of the separation of variables method, with the same parameters you've provided. Kindly make changes in the code as per your use case (if needed)
L = 20;
alpha = 0.23;
t_final = 60;
n = 20;
T0 = 20;
T1s = 100;
T2s = 0;
dx = L / n;
dt = 2;
x = linspace(0, L, n);
t = 0:dt:t_final;
nt = length(t);
% Precompute the Fourier modes
T = zeros(n, nt);
T(:, 1) = T0;
% Fourier expansion coefficients
for k = 1:n
lambda_k = (k * pi / L)^2;
Xk = sin(k * pi * x / L);
Theta_k = exp(-lambda_k * alpha * t);
% Superimpose the solutions
for j = 1:nt
T(:, j) = T(:, j) + (T1s - T2s) * Xk' * Theta_k(j);
end
end
% Plot the result
figure(1)
mesh(x, t, T.')
xlabel('Position (x)')
ylabel('Time (seconds)')
zlabel('Temperature, U(x,t)')
