hi all, maybe my previous post is a bit confusing. but i really need someone to answer me as soon as possible.

here is what i was trying to say:
if i were to solve a poisson equation (linear problem) using mglin code, in my main program, i pass the right hand side component of the poisson equation to the function mglin, and then the mglin function will return an array.
i wonder if this array is the solution of the possion equation already ?
Or do i need to call the function mglin again and again until the returned array is converged (which i think is not possible)?
Or do i need to modify the mglin code so that there is a new for-loop for convergent criteria?

without modifying the mglin code, by calling mglin function once, the resulting array i got seems to be not correct to my sample answer to the poisson i'm trying to solve.

so if anyone used the code before, please let me know what i should do if there is any trick of using the code.

many many thanks in advance.