Gauss-Seidel Linear Ecuation Solver (using epsilon value)

/////////////////////////////////////////////////////
// Hecho por XeBuZer0 y licenciado bajo GNU GPL v3 //
// Este programa calcula las soluciones numéricas  //
// de un sistema de ecuaciones arbitrario mediante //
// el método de Gauss-Seidel dado un valor épsilon //
// a alcanzar                                      //
/////////////////////////////////////////////////////
/* ** // **  F v q _ U k r a N a z i s !  ** // ** */
/////////////////////////////////////////////////////

#include <stdio.h>
#define epsilon 0.00000000001

long double VAbsol (long double a);

void main (void){

    int         iter = 0;
    long double     x1=.1, x2=.1, x3=.1;
    register long double    x1t=0, x2t=0, x3t=0;

    while ( VAbsol(x3t - x3) > epsilon || VAbsol( x2t - x2 ) > epsilon || VAbsol(x3t - x3) > epsilon){
        x3t = x3; x2t = x2; x1t = x1;
        x1 =  ( ( 5 - (2*x2t) )        / 7 ) ;
        x2 = -( ( (3*x1t) - x3t  - 4 ) / 5 ) ;
        x3 =  ( ( (5*x2t) - 3)         / 6 ) ;
        printf("Los valores son: x1=%.10Lf, x2=%.10Lf, x3=%.10Lf, Iteración %d\n", x1, x2, x3, ++iter);
    }
}

long double VAbsol (long double a){
    if (a<0) return (-1)*a;
    else return a;
}