%PALICA racuna in prikaze porazdelitev toplote v palici
%[x,t,U]=PALICA(temp,L,cas,podatki,nx,nt);
%x je delitev palice, t casovnega intervala, U porazdelitev toplote
%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
function [x,t,U]=palica(temp,L,cas,podatki,nx,nt);
kappa=podatki(3)/prod(podatki(1:2));
nedimcas=kappa*cas/L^2;
if nargin<6
	[x1,t1,U1]=nedim_toplota(nx,nedimcas);
else
	[x1,t1,U1]=nedim_toplota(nx,nedimcas,nt);
end;
x=L*x1;
kappa=podatki(3)/prod(podatki(1:2));
t=L^2/kappa*t1;
U=temp(1)+(temp(2)-temp(1))*U1;
clf;
surf(t,x,U);
shading flat;
