%----------------------------------------
% Numeerinen ja symbolinen laskenta 2005
% Harjoitus 6, tehtävä 4
%----------------------------------------

c=1;
L=10;

nx=50;
nt=40;
Tmax=50;
dt=Tmax/(nt-1);

u0=inline('zeros(size(x))','x');
ga=inline('t./(1+t) +t.*exp(-t/5)','t');
gb=inline('zeros(size(t))','t');

[x,t,Z]=CrNich(u0,ga,gb,[0,L],nx,nt,dt);

figure(1), clf
surf(x,t,Z)
xlabel('x'), ylabel('t')

figure(2), clf
Tplot=ceil( linspace(1,nt,15) );
plot(x(:),Z(Tplot,:)')


