%% Display Polynomial Fit w/ Reference Line
 
% Polynomial
imgSize = size(SourceImg);
p = polyfit(y, xNew, 3);
x = linspace(-50, imgSize(1), imgSize(1)+50);
f1 = polyval(p, x);
figure(2);
plot(f1, 'r', 'LineWidth', 10);
yticks([]);
xticks([]);
title('Severity of Spinal Curvature');
xlim([min(y)-50 max(y)+50]);
if max(abs(xNew)) >= min(abs(xNew))
   ylim([-max(abs(xNew))-500 max(abs(xNew))+500]);
else
   ylim([-min(abs(xNew))-500 min(abs(xNew))+500]);
end
hold on
 
% Zero line
f2 = polyval(0, x);
plot(f2, 'r');
hold on
