% Example taken from Chapter 8 of N. Trefethen's book 
% "Approximation Theory and Approximation Practice" (SIAM, 2013)

clear all, close all

f = chebfun('1./(1+25*x.^2)');
hold off, chebellipseplot(f,'linewidth',1)
hold on, plot([.2i -.2i],'xr','markersize',12)
axis equal, xlim(1.25*[-1 1]); ylim(.5*[-1 1]), grid on
FS = 'fontsize';
title('Chebfun ellipse for 1/(1+25x^2)',FS,9)
