
root = (3-Sqrt[5])/2;

f[a_,b_,c_]:=(
d=(-1+b c)/a;
test1=d-Floor[d];
test2=If[test1==0,0,1];
test3=(a root + b)/(c root+d) - root;
test3=test3 test3;
test4=If[test3<.00000001,0,1];
test5=test2+test4;
test5)

list1=Table[Table[Table[f[a,b,c],{a,1,16}],{b,-16,16}],{c,-16,16}]
list2=Table[Table[Table[{a,b,c}, {a,1,16}],{b,-16,16}],{c,-16,16}];



m={{2,-1},{1,-1}}