function [u,xi,tau]=brezdim_toplota(nxi,T,ntau);
%BREZDIM_TOPLOTA resuje brezdimenzionalni problem 
%prevajanja toplote
%[u,xi,tau]=BREZDIM_TOPLOTA(nxi,T,ntau)
%nxi je stevilo delilnih tock v prostorski smeri
%T je koncni cas
%ntau je stevilo delilnih tock v casovni smeri
  xi=linspace(0,1,nxi+2);
  dxi=xi(2)-xi(1);
  if nargin<3
  	%izberemo optimalni ntau
	lambda=1/6;
	dtau=lambda*dxi^2;
	ntau=ceil(T/dtau)%stevilo delilnih tock zaokrozimo navzgor
	tau=linspace(0,T,ntau);
	dtau=tau(2)-tau(1);
  else
	tau=linspace(0,T,ntau);
	dtau=tau(2)-tau(1);
	lambda=dtau/dxi^2;
  end
  u=zeros(ntau,nxi+2);
  u(:,nxi+2)=ones(ntau,1);
  for i=1:ntau-1
    for j=2:nxi+1
	  u(i+1,j)=lambda*u(i,j-1)+(1-2*lambda)*u(i,j)+lambda*u(i,j+1);
    end;
  end
