istr='(x-2)^2/(2^2) + (y+1)^2/(3^2) - 1';
ezplot(istr,[-2 6 -5 3])
axis square
grid
title(['Figure 25.23: ' istr])