--- a/scripts/polynomial/spline.m
+++ b/scripts/polynomial/spline.m
@@ -125,24 +125,24 @@
           - 3 * (a(n-1,:) - a(n-2,:)) / h(n-2);
 
       c(2:n-1,:) = spdiags ([[e(:); 0], dg, [0; e(:)]],
-			      [-1, 0, 1], n-2, n-2) \ g;
+                              [-1, 0, 1], n-2, n-2) \ g;
     endif
 
     c(1,:) = (3 / h(1) * (a(2,:) - a(1,:)) - 3 * dfs
-	      - c(2,:) * h(1)) / (2 * h(1));
+              - c(2,:) * h(1)) / (2 * h(1));
     c(n,:) = - (3 / h(n-1) * (a(n,:) - a(n-1,:)) - 3 * dfe
-		+ c(n-1,:) * h(n-1)) / (2 * h(n-1));
+                + c(n-1,:) * h(n-1)) / (2 * h(n-1));
     b(1:n-1,:) = diff (a) ./ h(1:n-1, idx) ...
-	- h(1:n-1,idx) / 3 .* (c(2:n,:) + 2 * c(1:n-1,:));
+        - h(1:n-1,idx) / 3 .* (c(2:n,:) + 2 * c(1:n-1,:));
     d = diff (c) ./ (3 * h(1:n-1, idx));
 
   else
 
     g = zeros (n-2, columns (a));
     g(1,:) = 3 / (h(1) + h(2)) ...
-	* (a(3,:) - a(2,:) - h(2) / h(1) * (a(2,:) - a(1,:)));
+        * (a(3,:) - a(2,:) - h(2) / h(1) * (a(2,:) - a(1,:)));
     g(n-2,:) = 3 / (h(n-1) + h(n-2)) ...
-	* (h(n-2) / h(n-1) * (a(n,:) - a(n-1,:)) - (a(n-1,:) - a(n-2,:)));
+        * (h(n-2) / h(n-1) * (a(n,:) - a(n-1,:)) - (a(n-1,:) - a(n-2,:)));
 
     if (n > 4)
 
@@ -157,7 +157,7 @@
       udg(1) = udg(1) - h(1);
       ldg(n - 3) = ldg(n-3) - h(n-1);
       c(2:n-1,:) = spdiags ([[ldg(:); 0], dg, [0; udg(:)]],
-			      [-1, 0, 1], n-2, n-2) \ g;
+                              [-1, 0, 1], n-2, n-2) \ g;
 
     elseif (n == 4)
 
@@ -165,10 +165,10 @@
       ldg = h(2) - h(3);
       udg = h(2) - h(1);
       c(2:n-1,:) = spdiags ([[ldg(:);0], dg, [0; udg(:)]],
-			      [-1, 0, 1], n-2, n-2) \ g;
+                              [-1, 0, 1], n-2, n-2) \ g;
       
     else # n == 3
-	    
+            
       dg = h(1) + 2 * h(2);
       c(2:n-1,:) = g/dg(1);
 
@@ -177,7 +177,7 @@
     c(1,:) = c(2,:) + h(1) / h(2) * (c(2,:) - c(3,:));
     c(n,:) = c(n-1,:) + h(n-1) / h(n-2) * (c(n-1,:) - c(n-2,:));
     b = diff (a) ./ h(1:n-1, idx) ...
-	- h(1:n-1, idx) / 3 .* (c(2:n,:) + 2 * c(1:n-1,:));
+        - h(1:n-1, idx) / 3 .* (c(2:n,:) + 2 * c(1:n-1,:));
     d = diff (c) ./ (3 * h(1:n-1, idx));
 
   endif