1> help itermeth
error: help: `itermeth' not found


1> pwd
ans = c:\MyOctave


2> help itermeth
`itermeth' is a function from the file c:\MyOctave\itermeth.m

ITERMETH    Un metodo iterativo generale
  [X, ITER] = ITERMETH(A,B,X0,NMAX,TOL,P) cerca di
  risovere iterativamente il sistema di equazioni
  lineari A*X=B su X. La matrice A di N-per-N coef-
  ficienti deve essere non singolare ed il termine
  noto B deve avere lunghezza N. Se P='J' viene usato
  il metodo di Jacobi, se P='G' viene invece selezio-
  nato il metodo di Gauss-Seidel. Altrimenti, P e'
  una matrice N-per-N non singolare che gioca il
  ruolo di precondizionatore in un metodo di
  Richardson a parametro dinamico.
  Il metodo si arresta quando il rapporto
  fra la norma del residuo corrente e quella del
  residuo iniziale e' minore di TOL e ITER e' il
  numero di iterazioni effettuate. NMAX prescrive
  il numero massimo di iterazioni consentite. Se P
  non viene precisata, viene usato il metodo di
  Richardson non precondizionato con parametro di
  rilassamento uguale a 1.

Additional help for built-in functions and operators is
available in the on-line version of the manual.  Use the command
`doc <topic>' to search the manual index.

Help and information about Octave is also available on the WWW
at http://www.octave.org and via the help@octave.org
mailing list.


3> A = hilb(13);


4> cond(A)
ans = 1.9487e+018


5> b = A * ones(13,1);


6> x = ones(13,1);


7> [xiter, iter] = itermeth(A, b, zeros(13,1), 10000, 1e-3, 'J')
xiter =

  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
   NaN
   NaN
   NaN
   NaN

iter =  304


8> [xiter, iter] = itermeth(A, b, zeros(13,1), 10000, 1e-3, 'G')
xiter =

   1.11235
   0.60005
   0.95398
   1.16379
   1.23100
   1.21824
   1.16550
   1.09485
   1.01814
   0.94156
   0.86826
   0.79970
   0.73642

iter =  34


9> [xiter, iter] = itermeth(A, b, zeros(13,1), 10000, 1e-9, 'G')
xiter =

   1.00028
   0.99653
   1.00557
   1.00959
   0.99070
   0.98590
   0.99315
   1.00376
   1.01158
   1.01362
   1.00902
   0.99814
   0.98182

iter =  10000


10> err_rel = norm(x - xiter)/norm(x)
err_rel =  0.0096728


11> [xiter, iter] = itermeth(A, b, randn(13,1), 10000, 1e-3, 'J')
xiter =

  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
  -Inf
   NaN
   NaN
   NaN
   NaN
   NaN
   NaN

iter =  304


12> [xiter, iter] = itermeth(A, b, randn(13,1), 10000, 1e-3, 'J')