无法从 sym 转换为 logical

clear;
clc;
xf = 0.85
xs = 0
ys = 0
p = 30
v = 0.25
e = 46
syms x y g g1 g2 r r1 r2;
g = atan((y-ys).*(pinv(x-xs)));
g1 = atan((y-ys).*(pinv(x+xf-xs)));
g2 = atan((y-ys).*(pinv(x-xs-xf)));
if(eval(g)<0)
    g = pi-g;
end
if(eval(g1)<0);
    g1 = pi-g1;
end
if(eval(g2)<0);
    g2 = pi-g2;
end
r = ((x-xs).^2+(y-ys).^2).^0.5;
r1 = ((x+xf-xs).^2+(y-ys).^2).^0.5;
r2 = ((x-xs-xf).^2+(y-ys).^2).^0.5;
ux = -(p*(1+v)/e).*((1-2*v).*((r1.*r2).^0.5.*cos((g1+g2)/2)-r.*cos(g))-r.^2.*pinv((r1.*r2).^0.5).*sin(g).*sin(g-(g1+g2)/2));
uy = -(p*(1+v)/e).*(2*(1-v).*((r1.*r2).^0.5.*sin((g1+g2)/2)-r.*sin(g))+r.*sin(g).*(1-r.*pinv((r1.*r2).^0.5).*cos(g-(g1+g2)/2)));
ax = diff(ux,x)
ay = diff(uy,y)
axy = 0.5*(ax+ay);
bx = e/((1+v)*(1-2*v))*((1-v)*ax+v*ay);
by = e/((1+v)*(1-2*v))*((1-v)*ay+v*ax);
bxy = e/((1+v)*(1-2*v))*(1-2*v)*axy;
c1 = 0.5*(bx+by)+(bxy.^2+0.25*((bx-by).^2))^0.5;
c2 = 0.5*(bx+by)-(bxy.^2+0.25*((bx-by).^2))^0.5;
cxy = 0.5*(c1-c2);
x = 0.85:0.01:4;
y = 0.85:0.01:4;
[x,y] = meshgrid(x,y);
c1 = eval(c1);
c2 = eval(c2);
cxy = eval(cxy);
figure
% fig1 = mesh(x,y,c1);
%fig1 = contourf(x,y,c1);
fig2 = contourf(x,y,c2);
% fig1 = contourf(x,y,cxy);
% fig2 = mesh(x,y,c2);
% fig3 = mesh(x,y,cxy);
xlabel("x轴")
ylabel("y轴")

知道是if那里出了问题，但是不知道怎么解

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。