Code covered by the BSD License  

Highlights from
Chebfun V4

image thumbnail

Chebfun V4

by

 

30 Apr 2009 (Updated )

Numerical computation with functions instead of numbers.

Editor's Notes:

This file was selected as MATLAB Central Pick of the Week

findexps_test
function pass = findexps_test
% Tests automatic exponent computation for functions 
% that diverge
% Nick Hale, Nov 2009

tol = chebfunpref('eps');
for k = 1:10   
    f = chebfun(@(x) 1./x.^k,[0 1],'blowup',1);
    g = chebfun(@(x) 1./x.^k,[-1 0],'blowup',1);
    pass(k) = ~(abs(f.exps(1)+k) + abs(g.exps(2)+k));
end

f = chebfun(@(x) x.^10,[0,1],'blowup','on');
pass(k+1) = (f.exps(1) == 0);

Contact us