%dispois_t.m
%poissoniana (dipendenza dal tempo per n fissato e vari lambda)
%mostra le proprietà di scaling nella dipendenza da t
figure(1);
clear;
n=6;
nt=81;
i1=1;

tmax=400;
lambda=0.05;
for k=1:4
	subplot(2,2,k);
	t=0:tmax/(nt-1):tmax;
	i2=size(t,2);
	imax=i2;
	for i=i1:i2
        y(i)=poisson(n,lambda*t(i));
        if y(i)>0.0001 imax=i; end
	end
	ymax=max(y);
	plot(t,y);
	axis([t(1) t(imax) 0 ymax*1.2]);
	[string,errmsg] = sprintf('%-7.2f',lambda);
	string1=['\lambda=',string];
	text(0.5*t(imax),ymax,string1);
	xlabel('t');
	lambda=2*lambda;
	tmax=tmax/2;
end
