clear all; 
close all;

% where we solve
a = 
b = 
% final time
T = 

% spatial step
n = 
h = (b-a)/n;

% time step
m = 
k = T/m;

% grid
u = zeros(m+1,n+1);
x = a:h:b ;
t = (0:k:T)' ;

% IC
u(1,:) = 

% solve time layers one by one
for j = 1:m
   
   % BC 1
   u(j+1,1) =  
   for i = 2:n
      u(j+1,i) = 
   end
   % BC 2
   u(j+1,n+1) =  
end

figure;
surf(x,t,u);
shading interp;
colorbar;
xlabel('x');
ylabel('t');
zlabel('u(x,t)');

