clc, clearvars
digits(60)
for t=0:64
for x=0:64
z=t-x;
if z>=0
y=round(log(vpa(2)^x+1-vpa(2)^z)/log(vpa(7)));
if vpa(2)^x+1==vpa(7)^y+vpa(2)^z && y>0
disp([x y z vpa(2)^x+1])
end
end
end
end
explores up to x+z=64 and finds
x y z LHS and RHS
[ 3, 1, 1, 9.0]
[ 6, 2, 4, 65.0]
|
Posted by Charlie
on 2021-07-27 11:29:33 |