normest
stima la norma 2 di una matrice. Fallisce sulla seguente matrice 4x4:
>> A=toeplitz([-2 1 0 0])
A =
-2 1 0 0
1 -2 1 0
0 1 -2 1
0 0 1 -2
>> norm(A)
ans =
3.6180
>> normest(A)
ans =
2.6180
Il problema dipende dalla scelta del
vettore iniziale per l'algoritmo delle potenze.quad
calcola numericamente un integrale con la formula di Simpson adattativa. Fallisce nel
calcolare il seguente integrale:
>> quad('x.^2.*exp(-x.^2)',-10,10)
ans =
0.8862
>> quad('x.^2.*exp(-x.^2)',-20,20)
ans =
2.3412e-07
>> quad('x.^2.*exp(-x.^2)',-30,30)
ans =
3.6975e-18
>> quad('x.^2.*exp(-x.^2)',-31,30)
ans =
0.8862
Il problema è nella stima dell'errore quando l'intervallo di integrazione è simmetrico.