function [u,x,t]=lampodiffE(u0,tloppu,r,c) 
% Sauva olkoon 1:n pituinen
% Esim: u_t=u_xx, alkuehto: u(x,0)=f(x)=4*x-4*x^2, 0-reunaehdot.
%  c=1;r=0.5;  % n=5, h=0.2
% f=inline('4*x.*(1-x)')
% u0=f(0:0.2:1);
% [U,x,t]=lampodiff(u0,0.2,r,c);
%
n=length(u0)-1;
h=1/n;  % x-askel
k=r*h^2/c^2; % t-askel (pisin, jolla stabiilisuus taattu) 
x=0:h:1; t=0:k:tloppu;
m=length(t);
u=zeros(n-1,m); % alueen sisällä vain n-1 x-jakopistettä
u(:,1)= u0(2:n)'; % jätetään päätyjen nollat pois 
M=lampomatr(n,r);
for j=1:m-1
   u(:,j+1)= M*u(:,j);
end
 % Liitetään lopuksi matriisin ylle ja 
 % alle nollarivi, jotta saadaan reuna-arvotkin mukaan.
nollarivi=zeros(1,m);
u=[nollarivi;u;nollarivi];




