clear all
close all
global KingCharles
global KingCharlesIndex
global Big
KingCharles = zeros(10000,1);
KingCharlesIndex = 1;
A = chebfun( @Aphra, [-1,1], 'minsamples', 129 )
%Warning: Function failed to evaluate on array inputs; vectorizing the function may speed up its
%evaluation and avoid the need to loop over array elements. Use 'vectorize' flag in the chebfun
%constructor call to avoid this warning message. 
%> In @chebfun\private\vectorcheck at 36
%%  In @chebfun\private\ctor_adapt at 117
% In chebfun.chebfun at 204
% In treachery at 5
%A = 
%   chebfun column (1 smooth piece)
%       interval          length    endpoint values   
%(       -1,        1)        1          0         0   
%
KingCharlesIndex
%
%KingCharlesIndex =
%
%    133
%
Big = 1.0e77;
B = chebfun( @Benedict, [-1,1], 'minsamples', 129 )
%B = 
%   chebfun column (1 smooth piece)
%       interval          length    endpoint values   
%(       -1,        1)        1          0         0   
%vertical scale =   0 
% Heh.  Same as Aphra.  She has done her job well! 
t = linspace(-1,1,2012);
Bt = Benedict(t');
Bct = B(t');
plot( t', Bct, 'r--', t', Bt, 'k.' )
