% numero di anni necessari per calcolare un
% determinante di ordine 31 con il supercomputer
% Titan (20 Peta-flops) mediante lo sviluppo di Laplace
% l'eta' dell'universo: 13,72 miliardi di anni
factorial(31)/(20.e15*60*60*24*365)
% Definire funzioni e plotare i loro grafici
% funzione inline
f = inline('sin(x)./x')
f(pi/2)
2/pi
f([pi/2,pi])
x = linspace(-pi,pi);
y = f(x);
size(y)
plot(x,y)
grid on
axis equal
% function handle
g = @(x)[sin(x)./x]
g(pi/2)
x = linspace(-2*pi, 2*pi);
plot(x, g(x))
x = linspace(-10*pi, 10*pi);
plot(x, g(x))
grid on
hold on
plot(x, sin(x), 'r')
hold off
% m-file di una function: funzione.m
function y = funzione(x)
y = sin(x)./x;
end
funzione(pi/2)
% funzioni di due variabili
f2 = inline('x.^2 + y.^2', 'x', 'y')
f2(1,2)
[x, y] = meshgrid([-1 1 3], [2 4])
z = f2(x,y)
mesh(x,y,z)
xx = linspace(-1,1);
[x, y] = meshgrid(xx, xx);
z = f2(x,y);
mesh(x,y,z)
% con le curve di livello
meshc(x,y,z)
% campo vettoriale gradiente
y = x = -1 : 0.2 : 1;
quiver(x, y, 2*x, 2*y)
yy = xx = -1 : 0.2 : 1;
[x, y] = meshgrid(xx, xx);
quiver(x, y, 2*x, 2*y)
quiver(x, y, -2*x, -2*y)
gg = @(x,y)[x.^2 - y.^2]
meshc(x,y,gg(x,y))
quiver(x, y, 2*x, -2*y)
quit
# Octave 3.6.2, Wed Oct 31 16:57:01 2012 CET <achilles@achilles-TOSH>

