eps
format long
eps
2^(-52)
1 + eps
if (1 + eps == 1) disp('si') else disp('no') end
if (1 + eps/2 == 1) disp('si') else disp('no') end
eps/2
realmin
realmin^2
pwd()
x = 1.e-15
((x + 1) -1)/x
x = [-2 -1 0 1 2]
y = 2 .^ x
format
y = 2 .^ x
plot(x,y)
x = linspace(-2,2)
x = linspace(-2,2);
size(x)
length(x)
y = 2 .^ x;
length(y)
plot(x,y)
hold on
y3 = 3 .^ x;
plot(x,y3, 'r')
axis equal
axis equal off
yexp = exp(x);
plot(x,yexp, 'g')
axis on
grid on
exp(1)
format long
e
pi
% occorre caricare la funzione costante
% http://www.dm.unibo.it/~achilles/calc/programmi/N_A.m
urlwrite("http://www.dm.unibo.it/~achilles/calc/programmi/N_A.m", "N_A.m")
N_A
hold off
x = linspace(1-1.e-5,1+1.e-5);
y = (x - 1).^3;
yy = x.^3 - 3*x.^2 + 3*x -1;
plot (x,y)
hold on
plot (x,yy, 'r')
hold off
quit
# Octave 3.6.4, Mon Oct 14 13:00:58 2013 W. Europe Daylight Time <unknown@achilles-TOSH>

