Funzione epsilon

FUNZIONE EPSILON

Scopo: la funzione calcola l'epsilon macchina o l'epsilon relativo a x.

Specifica: l'intestazione della funzione è e=epsilon(x).

Parametri: - x è il valore a cui si vuole calcolare epsilon. È un parametro di ingresso.

- e è il valore restituito dalla funzione, cioè l'epsilon macchina o l'epsilon relativo a x.

È un parametro di uscita.

Descrizione: se viene inserito il parametro di ingresso(x) la funzione calcolerà l'epsilon macchina,  altrimenti calcolerà l'epsilon relativo a x.

Indicatori di errore: Too many input arguments= avviene se si inseriscono più parametri di    ingresso.(Questo indicatore di errore è definito dal Matlab)

Formato non corretto di x. Inserirlo come numero reale = indica che è stato inserito un carattere

Esempio d'uso: >> x=4

>> e=epsilon(x)

e = 8.881784197001252e-016

Situazioni di errore: >> epsilon(3,5)

??? Error using ==> epsilon

Too many input arguments.

>> epsilon('a')

??? Error using ==> epsilon

Formato non corretto di x. Inserirlo come numero reale

Test dei casi funzionanti: come già detto, la funzione può essere utilizzata per due diversi motivi.

Per calcolare l'epsilon macchina: >> e=epsilon 

e = 2.220446049250313e-016

O per calcolare l'epsilon relativo a x: >> x=4

>> e=epsilon(x)

e = 8.881784197001252e-016

Confronto con Matlab: con la funzione EPSILON si ha che l'epsilon macchina è 2.220446049250313e-016  con la funzione EPS del Matlab si ha che l'epsilon macchina è 2.220446049250313e-016.

M-FILE:

function e=epsilon(x)

% EPSILON funzione che calcola l'epsilon macchina o l'epsilon relativo a x.                                                                                                              

% parametri di ingresso:    

x=valore a cui si vuole calcolare epsilon

% parametri di uscita:                                

e=epsilon macchina o epsilon relativo a x

% note:                                                                          

se non c'è input la funzione calcolerà l'epsilon macchina

% esempio d'uso:                                                                                            

1) >> e=epsilon

e = 2.220446049250313e-016

2) >> x=4

>> e=epsilon(x)

e = 8.881784197001252e-016

if (nargin==0)       %controlla se c'è input

e=1;

while(1+e~=1)

e=e/2;

end

e=2*e;

else            %se c'è input calcola l'epsilon relativo a x

if (ischar(x)==1) %Controlla se x è un numero reale

error('Formato non corretto di x. Inserirlo come numero reale');

end

e=eps*abs(x);

end