%dispois_lt.m
%poissoniana (dipendenza dal tempo per vari n)
figure(1);
clear;
k=6;
nt=301;
i1=1;

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