三维曲线作图

8*y^2-6*x*y+4*x^2-(4*y-x)^3=0; (隐函数方程)                                                          x<y;                         （不等式）
        z=4*y^2*(y-x)/(4*y-x)^2;         
   
请问大家，这个怎么作三维曲线？ 我的想法如下，但是老报错，不清楚是不是格式问题，请教大家：

x=0:0.001:0.5; %穷举作图
y=solve('(8*y^2-6*x*y+4*x^2-(4*y-x)^3)/(x<y)','x'); %求出满足不等式的隐函数解
z=4*y^2*(y-x)/(4*y-x)^2;
plot3(x,y,z);
grid on;

图已画好，一共用了两个m file，其中一个是主程序，另外一个是用来解隐函数的函数文件，都在rar文件里。隐函数求解用的是fsolve的命令。注意，solve命令只能用来解符号运算方程，数值运算应该使用fsolve命令。当求解每一个x所对应的y的时候，我发现当x超过0.225左右之后x<y的约束就无法满足了。第一幅图是你要画的三维图，我使用view命令把三维图的视角调整了一下，默认的角度效果不是很好。第二幅图中的红线是x=y，而蓝线是隐函数所对应的y，这幅图说明当x足够晓得时候约束才能够满足。请查阅附件后付论坛币。

liuturtle 发表于 2014-11-25 19:37
图已画好，一共用了两个m file，其中一个是主程序，另外一个是用来解隐函数的函数文件，都在ra ...

感谢，对我很有启发