syms pm
f = 8998028467180289/2251799813685248 - ((pm - 12844417258516453/1125899906842624)^2 + 4055644986144483/147573952589676412928)^(1/2)/2 - pm/2;
k = double(solve(f==0));
if isempty(k)
disp(11)
else
disp(22)
end
方程f=0,就是一个无解的方程,如下图:
直接求解得到的结果是,k = Empty sym: 0-by-1
将k转为实数,然后使用isempty函数判断就好了。