How to solve heat equation on matlab ?

dt=0.1;
dx=1;
D=0.002;%thermal coefficient
s=(D*0.1)/(1)^2; %s=D*dt/dx^2
T(1,1)=0.2;%initial condition 
t=0:0.1:1;
x=0:1:10;

for n=1:length(t-1)
    for j=2:length(x-1)
    T(n+1,j)=s*T(n,j+1)+(1-2*s)*T(n,j)+s*T(n,j-1);
    end

end

but the error I,m getting is in the matrix dimension :

Index exceeds matrix dimensions.

Error in untitled (line 13) T(n+1,j)=s*T(n,j+1)+(1-2*s)*T(n,j)+s*T(n,j-1);

[KSSV]

Your T is a scalar..and you re trying to access it as vector......so the error.

[KSSV]

May be helpful: https://in.mathworks.com/matlabcentral/fileexchange/42604-finite-difference-method-to-solve-heat-diffusion-equation-in-two-dimensions

Answer 1

[Torsten]

Use MATLAB's "pdepe".

Best wishes

Torsten.