function values = sin_partial(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
sins0 = [T0,T0+20,T0-20];
sinsL = [TL,TL+20,TL-20];
values = zeros(t_steps,x_steps);

t = [0:delta_t:tmax];
A = 20;
f =1;
y = A*sin(2*pi*f*t)+100;



for i=1:size(values,1)
  values(i,1) = y(i);
endfor


t = [0:delta_t:tmax];
A = 20;
f =1;
y = A*sin(2*pi*f*t)+50;

for i=1:size(values,1)
  values(i,size(values,2)) = y(i);
endfor


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