Usage: |
{A,V,ctrl}=itera(A0,FI,b,N,eps,itmax,expl)
|
Input: |
| A0 | p x p matrix, starting value for A;
if 0 is entered I^{-1} is taken
if p>1 and dim(A0)==1, I^{-1}*A0 is taken
|
| FI | p x p matrix, Fisher-Information
|
| b | numeric, clipping height
|
| N | scalar, MC-sample size or integration grid-points
|
| eps | numeric, exactitude
|
| itmax | scalar, maximal number of iterations
|
| expl | numeric, threshold for the changes in absolute
value of A: beyond this value convergence is uncertain ~ 4
|
Output: |
| A | p x p matrix, Lagrange-Multiplier solving (1) |
| V | p x p matrix of corresponding covariance |
| ctrl | scalar, tells if convergence "happened" |