Changeset 19341
- Timestamp:
- 05/05/15 16:45:47 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/psl/p_polynomial_value.m
r19339 r19341 103 103 end 104 104 105 v = zeros ( m, n + 1 ); 105 v = ones ( m, n + 1 ); 106 %v = zeros ( m, n + 1 ); 106 107 107 v(1:m,1) = 1.0;108 % v(1:m,1) = 1.0; 108 109 109 110 if ( n < 1 ) … … 111 112 end 112 113 113 v(1:m,2) = x(1:m,1); 114 v(1:m,2) = x; 115 %v(1:m,2) = x(1:m,1); 114 116 115 117 for i = 2 : n 116 118 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) ) ... 119 123 / ( i ); 120 124 121 125 end 122 126 123 return124 end127 % return 128 %end
Note:
See TracChangeset
for help on using the changeset viewer.