/*  solvtd.c    CCMATH mathematics library source code.
 *
 *  Copyright (C)  2000   Daniel A. Atkinson    All rights reserved.
 *  This code may be redistributed under the terms of the GNU library
 *  public license (LGPL). ( See the lgpl.license file for details.)
 * ------------------------------------------------------------------------
 */
void solvtd(double *a, double *b, double *c, double *x, int m)
{
    double s;

    int j;

    for (j = 0; j < m; ++j) {
	s = b[j] / a[j];
	a[j + 1] -= s * c[j];
	x[j + 1] -= s * x[j];
    }
    for (j = m, s = 0.; j >= 0; --j) {
	x[j] -= s * c[j];
	s = (x[j] /= a[j]);
    }
}