f=inline('1./(1+x.^2)');
x=[-5:0.01:5]; y=f(x); plot(x,y);
hold on;
n=8; xi=linspace(-5,5,n+1); yi=f(xi); p=polyfit(xi,yi,n);
yp=polyval(p,x);
plot(x,yp,'g',xi,yi,'ro');
% nodi di Chebyshev
xc=-5*cos(pi/8*[0:8]);
yc=f(xc); pc=polyfit(xc,yc,n); ypc=polyval(pc,x);
plot(x,ypc,'c',xc,yc,'k+');
hold off;
% interpolazione lineare composita
plot(x,y,x,interp1(xi,yi,x),'r');
hold on;
% la spline cubica interpolatoria con 'not-a-knot end conditions'
plot(x,spline(xi,yi,x),'c')
% numeri (pseudo-) casuali uniformemente distributi tra 0 e 1
rand(1,1e6);
hist(rand(1,1e6))
hold off
hist(rand(1,1e6))
hist(rand(1,100000))
hist(rand(1,10000))
hist(2*rand(1,10000))
% l'istogramma di numeri casuali uniformemente distributi tra -1 e 1
hist(2*rand(1,10000)-1)
% l'istogramma di numeri casuali normalmente distribuiti
% distribuzione di Gauss con media 0 e deviazione standard 1
hist(randn(1,10000))
% istogramma con 30 colonne
hist(randn(1,1e6),30)
hist(randn(1,1e6),100)
hist(2*randn(1,1e6),100)
hold on;
hist(0.5*randn(1,1e6),100)
hold off;
hist(2*randn(1,1e6),100)
hold on
hist(0.5*randn(1,1e6),100,'c')
y=1:5
y+0.01* (2*rand(1,5)-1)a
% 5 lanci di un dado a 6 facce
fix(6*rand(1,5))+1
fix(6*rand(1,5))+1
fix(6*rand(1,5))+1
fix(6*rand(1,5))+1
fix(6*rand(1,5))+1
fix(pi)
fix(-pi)
fix(3.9)
fix(-3.9)
rand('seed',1)  % inizializza il seme del generatore dei numeri casuali
                % (per ottenere con il 'rand' successivo sempre gli 
                % stessi numeri)
rand(1,10)
quit
