kill(all)$

van=facebook.com/LearningCalculus$

f(a,b):=a^2/(2*a*b^2 - b^3 + 1)$
f(x,y);

for n: 0 thru 25 step 5 do
(
    count: 0,
    for i: 1 thru n do
    (
        for j: 1 thru n do
        (
            if is(integerp(f(i,j))) then
            (
                count: count + 1
            )
        )
    ),
    printf(true,
"Answer Solutions: Possible: ~3d Actual: ~3d ~%",
    n*n, count)
)$

for n: 10 thru 10 do
(
    count: 0,
    for i: 1 thru n do
    (
        for j: 1 thru n do
        (
            if is(integerp(f(i,j))) then
            (
                count: count + 1,
                printf(true,
"Answer i, j, f(i,j): ~3d ~3d ~3d ~%",
                 i,j,f(i,j))
            )
        )
    )
)$

f(2*m,1);

f(5*m,1);