syms x yresult=[];for k=-1:0.025:1 result=[result solve(subs((x^2+y^2+1)^3-x^2 * y^3,k))];end这样求出来必然会出现虚数解，不知是否是您想要的%subs是将x=k带入方程,得到关于y的方程，然后solve求解关于y的方程
