function U=lampodiffB(RAA,L,nx,nt,dt)
% Lämpöyhtälön ratkaisu implisiittisellä menetelmällä.
% Syötteet:
% RAA -- Reuna-ja alkuarvot, funktio, joka palauttaa kolme tulosta...
% Esim:
%   function [vreuna,aa,oreuna]=RAA(x,t)
%   vreuna=1./t; aa=(sin(pi*x/2)).^2; oreuna=t;
% L,nx,nt,dt   kuten aina.
%
% Testaa sitä mukaan kuin kirjoitat rivejä:
% L=1; nx=10; dt=0.1; nt=3;
%  A=lampodiffB(@RAA,L,nx,nt,dt); full(A)    % Kohta 1

dx=1/(nx+1);
r=dt/dx^2;
e=ones(nx,1);
A=spdiags([-e 2*e -e],-1:1,nx,nx);
U=A  % Kohta 1. 
