diff --git a/tests/TestContactAngle.cpp b/tests/TestContactAngle.cpp
index a6f7f9fe..e3a0f35a 100644
--- a/tests/TestContactAngle.cpp
+++ b/tests/TestContactAngle.cpp
@@ -200,29 +200,29 @@ int main (int argc, char *argv[])
 	printf("-------------------------------- \n");
 	printf("NWP volume = %f \n", nwp_volume);
 	analytical = 2*PI*(CAPHEIGHT*CAPHEIGHT+RADIUS*RADIUS);
-	RelError = (awn-analytical)/analytical;
+	RelError = fabs(awn-analytical)/analytical;
 	printf("Area wn = %f, Analytical = %f, Rel. Error = %f \n", awn,analytical,RelError);
 	analytical = 2*PI*RADIUS*(N-2)-4*PI*RADIUS*(CAPRAD-CAPHEIGHT);
-	RelError = (ans-analytical)/analytical;
+	RelError = fabs(ans-analytical)/analytical;
 	printf("Area ns = %f, Analytical = %f, Rel. Error = %f \n", ans, analytical,RelError);
 	analytical = 2*PI*RADIUS*(CAPRAD-CAPHEIGHT);
-	RelError = (aws-analytical)/analytical;
+	RelError = fabs(aws-analytical)/analytical;
 	printf("Area ws = %f, Analytical = %f, Rel. Error = %f \n", aws, analytical,RelError);
 	analytical = 2*PI*RADIUS*(N-2);
-	RelError = (As-analytical)/analytical;
+	RelError = fabs(As-analytical)/analytical;
 	printf("Area s = %f, Analytical = %f, Rel. Error = %f \n", As, analytical,RelError);
 	analytical = 4*PI*RADIUS;
-	RelError = (lwns-analytical)/analytical;
+	RelError = fabs(lwns-analytical)/analytical;
 	printf("Length wns = %f, Analytical = %f, Rel. Error = %f  \n", lwns, analytical,RelError);
 	analytical = 1.0*RADIUS/CAPRAD;
-	RelError = (efawns-analytical)/analytical;
+	RelError = fabs(efawns-analytical)/analytical;
 	printf("Cos(theta_wns) = %f, Analytical = %f, Rel. Error = %f  \n",efawns,analytical,RelError);
-	analytical = 2.0/RADIUS;
-	RelError = (Jwn-analytical)/analytical;
+	analytical = 2.0/CAPRAD;
+	RelError = fabs(Jwn-analytical)/analytical;
 	printf("Mean curvature (wn) = %f, Analytical = %f, Rel. Error = %f  \n", Jwn, analytical,RelError);
 	printf("Geodesic curvature (wns) = %f, Analytical, Rel. Error = %f  = %f \n", KGwns, 0.0, 0.0);
 	analytical = 1.0/RADIUS;
-	RelError = (KNwns-analytical)/analytical;
+	RelError = fabs(KNwns-analytical)/analytical;
 	printf("Normal curvature (wns) = %f, Analytical = %f, Rel. Error = %f  \n", KNwns, analytical,RelError);
 	printf("-------------------------------- \n");	
 	//.........................................................................