function values = partial11(T0,TL,k,xmax,tmax,delta_x,delta_t)


t_steps = (tmax/delta_t)+1
x_steps = (xmax/delta_x)+1
xpoints = []
tpoints = []

for x=0:x_steps-1
  xpoints = [xpoints,x];
endfor
xpoints
for t=0:t_steps-1
  tpoints = [tpoints,t];
endfor
tpoints

values = zeros(t_steps,x_steps);
values(:,1) = T0;
values(:,size(values,2)) = TL;
for t=1:t_steps-1
  for x=2:x_steps-1
  values(t+1,x) = values(t,x) + (k^2)*(delta_t)*(delta_x^-2)*(values(t,x+1) - 2*values(t,x) + values(t,x-1));
  endfor
endfor
values;
mesh(xpoints, tpoints, values);
xlabel("x");
ylabel("t");
zlabel("T");
title ("3-D Plot");
endfunction