function [U,x,t]=palica(temp,L,cas,podatki,nx,nt);
%PALICA racuna in prikaze porazdelitev toplote v palici
%[U,x,t]=PALICA(temp,L,cas,podatki,nx,nt);
%U je porazdelitev toplote, x delitev palice, t casovnega intervala,
%temp=[T0,T1]; temperatura na zacetku in koncu palice
%L dolzina
%cas, do katerega gledamo
%podatki=[ro,c,k]
%k konduktivnost, ro gostota, c specificna toplota
%nx st. delitev v smeri x, nt st. delitev po casu
kappa=podatki(3)/prod(podatki(1:2));
nedimcas=kappa*cas/L^2;
if nargin<6
	[U1,x1,t1]=brezdim_toplota(nx,nedimcas);
else
	[U1,x1,t1]=brezdim_toplota(nx,nedimcas,nt);
end;
x=L*x1;
kappa=podatki(3)/prod(podatki(1:2));
t=L^2/kappa*t1;
t(end)=cas;
U=temp(1)+(temp(2)-temp(1))*U1;
clf;
size(U)
size(x)
size(t)
surf(x,t,U);
shading flat;
