1.priklad intervalovy odhad strednej hodnoty 99% sigma znama

>> rand('state',2)					% nastavenia generatora
>> S2=9;						% rozptyl
>> n=25;						% rozsah nahodneho vyberu
>> x=rand(1,n)*30;					% nahodny vyber
>> S=sqrt(S2);						% odmocnina rozptylu
>> F=inline('0.5*erf(sqrt(0.5)*x)+0.5');		% distribucna funckia normaloveho rozdelenia
>> F0995=inline('0.5*erf(sqrt(0.5)*x)+0.5-0.995')	% 
>> u=fzero(F0995,0)					% 
>> xx=mean(x);						% priemer nahodneho vyberu
>> d=S*u/sqrt(n);					% pomocny vypocet
>> [xx-d,xx+d]						% interval spolahlivosti
	

2.priklad

>> f=inline('x.^3.*exp(-x)+x.^4');			% zadefinovanie funkcie
>> x=-2:0.01:2;						% interval
>> grid on						% mriezka
>> plot(x,f(x))						% graf fx na x
>> hold on						% podrzi graf
>> ezplot('x.^3.*exp(-x)+x.^4',[-2,2])			% graf cez ezplot
>> fzero(f,[-0.0005,0.001])				% hodnota korena cez fzero
>> a=-0.0005; b=0.001; 					% interval na ktorom ma hladat koren
>> for k=1:40 						% kolko krokov taka je presna napr 2^(-40) je dakych na -13
>> s=(a+b)/2; if f(s)<0 a=s; else b=s; end, end, s	% bisekcia, vypise koren