function [Y] = findExact( t )

Y = 1 - sqrt(1 + t*t*t + 2*t*t + 2*t + 3);