Index: /issm/trunk-jpl/src/m/psl/p_polynomial_value.m
===================================================================
--- /issm/trunk-jpl/src/m/psl/p_polynomial_value.m	(revision 19340)
+++ /issm/trunk-jpl/src/m/psl/p_polynomial_value.m	(revision 19341)
@@ -103,7 +103,8 @@
   end
 
-  v = zeros ( m, n + 1 );
+  v = ones ( m, n + 1 );
+  %v = zeros ( m, n + 1 );
 
-  v(1:m,1) = 1.0;
+%  v(1:m,1) = 1.0;
 
   if ( n < 1 )
@@ -111,14 +112,17 @@
   end
 
-  v(1:m,2) = x(1:m,1);
+  v(1:m,2) = x;
+  %v(1:m,2) = x(1:m,1);
 
   for i = 2 : n
  
-    v(1:m,i+1) = ( ( 2 * i - 1 ) * x(1:m,1) .* v(1:m,i)   ...
-                -  (     i - 1 ) *             v(1:m,i-1) ) ...
+    %v(1:m,i+1) = ( ( 2 * i - 1 ) * x(1:m,1) .* v(1:m,i)   ...
+    %            -  (     i - 1 ) *             v(1:m,i-1) ) ...
+    v(:,i+1) = ( ( 2 * i - 1 ) * x .* v(:,i)   ...
+                -  (     i - 1 ) *             v(:,i-1) ) ...
                 /  (     i     );
  
   end
  
-  return
-end
+%  return
+%end