Changeset 19341


Ignore:
Timestamp:
05/05/15 16:45:47 (10 years ago)
Author:
adhikari
Message:

minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/psl/p_polynomial_value.m

    r19339 r19341  
    103103  end
    104104
    105   v = zeros ( m, n + 1 );
     105  v = ones ( m, n + 1 );
     106  %v = zeros ( m, n + 1 );
    106107
    107   v(1:m,1) = 1.0;
     108%  v(1:m,1) = 1.0;
    108109
    109110  if ( n < 1 )
     
    111112  end
    112113
    113   v(1:m,2) = x(1:m,1);
     114  v(1:m,2) = x;
     115  %v(1:m,2) = x(1:m,1);
    114116
    115117  for i = 2 : n
    116118 
    117     v(1:m,i+1) = ( ( 2 * i - 1 ) * x(1:m,1) .* v(1:m,i)   ...
    118                 -  (     i - 1 ) *             v(1:m,i-1) ) ...
     119    %v(1:m,i+1) = ( ( 2 * i - 1 ) * x(1:m,1) .* v(1:m,i)   ...
     120    %            -  (     i - 1 ) *             v(1:m,i-1) ) ...
     121    v(:,i+1) = ( ( 2 * i - 1 ) * x .* v(:,i)   ...
     122                -  (     i - 1 ) *             v(:,i-1) ) ...
    119123                /  (     i     );
    120124 
    121125  end
    122126 
    123   return
    124 end
     127%  return
     128%end
Note: See TracChangeset for help on using the changeset viewer.