/**
	Fraga, 30/10/2008
	Calcula el vector x
	tal que Ax = s, donde A es una matrix tridiagonal
	cuyos elementos estan almacenados en 
	a[], b[] y c[]
**/
void tridiagonal( float a[], float b[], float c[], float s[],
	float x[], int n )
{
	int i, j;
	float r;

	for( i=1; i<n; i++ ) {
		j = i - 1;
		r = a[i]/b[j];
		b[i] -= r*c[j];
		s[i] -= r*s[j];
	}
	j = n - 1;
	x[ j ] = s[j]/b[j];
	for( i=j-1; i>=0; i-- )
		x[i] = (s[i] - c[i]*x[i+1])/b[i];
}