Mercurial > octave-nkf

comparison scripts/general/quadgk.m @ 11985:bdc383a457fb release-3-2-x

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix typo in previous change
author John W. Eaton <jwe@octave.org>
date Thu, 11 Jun 2009 07:13:19 +0200
parents b4986fecdc53
children b18a50c56144
comparison
equal deleted inserted replaced
11984:b4986fecdc53 11985:bdc383a457fb
193 ## g(t) = t / (1 - t^2) 193 ## g(t) = t / (1 - t^2)
194 ## g'(t) = (1 + t^2) / (1 - t^2) ^ 2 194 ## g'(t) = (1 + t^2) / (1 - t^2) ^ 2
195 ## waypoint transform is then 195 ## waypoint transform is then
196 ## t = (2 * g(t)) ./ (1 + sqrt(1 + 4 * g(t) .^ 2)) 196 ## t = (2 * g(t)) ./ (1 + sqrt(1 + 4 * g(t) .^ 2))
197 if (!isempty (waypoints)) 197 if (!isempty (waypoints))
198 trans = @(x) (2 * x) ./ ((1 + sqrt(1 + 4 * x .^ 2)); 198 trans = @(x) (2 * x) ./ (1 + sqrt(1 + 4 * x .^ 2));
199 subs = [-1; trans(waypoints); 1]; 199 subs = [-1; trans(waypoints); 1];
200 else 200 else
201 subs = linspace (-1, 1, 11)'; 201 subs = linspace (-1, 1, 11)';
202 endif 202 endif
203 h = 2; 203 h = 2;