求解方程组问题，画图和实际不一致

ywh19860616

各位老师，好
我想求解一个方程组
clear all;
beta=1;
p=2;
a=3;
eta=3;
gamma=4;
alpha=0.5;
thet=2;
syms vr vd
eq1=beta.*(p-(eta.*vr.^2)./2)-gamma.*(p-(eta.*vd.^2)./2)+eta.*vr.*(alpha.*p-beta.*vr+gamma.*vd+a.*(thet-1))
eq2=beta.*(p -(eta.*vd.^2)./2)-gamma.*(p-(eta.*vr.^2)./2)+eta.*vd.*(alpha.*p-a.*thet-beta.*vd+gamma.*vr)
ezplot(eq1)
hold on
ezplot(eq2)

为了简单，我先把除vr，vd变量外的其他变量赋值，从上面图形看，应该有交点，可是下面用solve求根，却提示
没有
syms vr vd
solve('beta*(p-(eta*vr^2)/2)-gamma*(p-(eta*vd^2)/2)+eta*vr*(alpha*p-beta*vr+gamma*vd+a*(thet-1))','beta*(p-(eta*vd^2)/2)-gamma*(p-(eta*vr^2)/2)+eta*vd*(alpha*p-a*thet-beta*vd+gamma*vr)','vr','vd')

Warning: Explicit solution could not be found.
> In solve at 140

ans =

[ empty sym ]

哪里出错了？