matlab程序：

solution=solve('(4*x^2-3*x*y+2*y^2)-k*x^3*(4*x-7*y)=0', 'y')
f=solution(1);%% f=solution(2);
x=linspace(1,10,25);%%参数范围可选
k=linspace(1,10,25);
[xx,kk]=meshgrid(x,k);
yy=xx;
for i=1:1:length(x)
  for j=1:1:length(k)
  yy(i,j) = subs(f, '[x,k]', [x(i),k(j)]);
end
end

mesh(xx,kk,yy)
xlabel('x');ylabel('k');zlabel('y')

\left\{\left\{y\to \frac{1}{4} \left(-x \sqrt{49 k^2 x^4-10 k x^2-23}-7 k x^3+3 x\right)\right\},\left\{y\to \frac{1}{4} \left(x \sqrt{49 k^2 x^4-10 k x^2-23}-7 k x^3+3 x\right)\right\}\right\}