horzLine = '____________________________________________';
fprintf('\n         Homework #4, Problem 2b \n%s \n', horzLine);
fprintf('  x    y         (x+y)^15=1          error \n%s\n', horzLine);

n=15;
for x=1:20
    y = 1 - x;
    result = binomial(x, y, n);
    fprintf(' %2d  %3d   %22.15e   %.1e \n', x, y, result, abs(result-1) );
end

fprintf('%s\n  ****  Only works for x below 7! \n\n', horzLine);
