#ifdef HAVE_STDLIB_H
#include <stdlib.h>
#endif
#include "yagi.h"

void mprove(double **a,double **alud,int n,int *indx,double *b,double *x)
{
	int j,i;
	double sdp;
	double *r;

	r=dvector(1,n);
	for (i=1;i<=n;i++) {
		sdp = -b[i];
		for (j=1;j<=n;j++) sdp += a[i][j]*x[j];
		r[i]=sdp;
	}
	lubksb(alud,n,indx,r);
	for (i=1;i<=n;i++) x[i] -= r[i];
	free_dvector(r,1,n);
}