kill(all)$

van=facebook.com/LearningCalculus$

for n: -7 thru 7 do
    if n # -1 then
    (
        left:bfloat(n^(n+1)),
        rite:bfloat((n+1)^n),
        printf(true,"        i=~2d: ~9f <? ~9f ~a~%",
               n, left, rite, is(left < rite))
    )

    else if n = -1 then
         printf(true,"        i=~2d: ~9f <? ~a ~a~%",
               -1, -1, "      1/0", "undefined");

van=facebook.com/LearningCalculus$ fpprec:20$

f(x):=x^(x+1)-(x+1)^x;

newton(t):=bfloat(subst(t,x,x-f(x)/diff(f(x),x)))$

newton(2.2931662874118610315b0);