--- a/inst/lombnormcoeff.m
+++ b/inst/lombnormcoeff.m
@@ -23,11 +23,13 @@
 ## @end deftypefn
 
 
-function coeff = lombnormcoeff(T,X,omega)
-tau = atan2( sum( sin( 2.*omega.*T)), sum(cos(2.*omega.*T))) / 2;
-coeff = ( ( sum ( X .* cos( omega .* T - tau ) ) .^ 2 ./ sum ( cos ( omega .* T - tau ) .^ 2 )
-	   + sum ( X .* sin ( omega .* T - tau ) ) .^ 2 / sum ( sin ( omega .* T - tau ) .^ 2 ) )
-	 / ( 2 * var(X) ) );
+function coeff = lombnormcoeff (T, X, omega)
+
+tau = atan2 (sum (sin (2 .* omega .*T)), sum (cos (2 .* omega .* T))) / 2;
+coeff = ((sum (X .* cos (omega .* T - tau)) .^ 2 / sum (cos (omega .* T - tau) .^ 2)  + 
+          sum (X .* sin (omega .* T - tau)) .^ 2 / sum (sin (omega .* T - tau) .^ 2 )) / 
+         (2 * var(X)));
+
 endfunction
 
 %!shared t, x, o, maxfreq