function fp = henon_fp( a, b )

fp = zeros(2,1);
fp(1,1) = (-(1-b)+sqrt((1-b)^2+4*a))/(2*a);
fp(2,1) = b*fp(1,1);