%
% Figures 20.10, 20.11, 20.12
%
function f = func2( x )
a = 9;
b = 10;
dd = ((a*11.36*x(2)-a*b-a*a*x(1))/(a*a+1))^2 + ...
	       ((a*x(1)+b-11.36*x(2))/(a*a+1))^2;
f = dd + 0.25*( x(1) - 2.5 )^2;
%
