有意思的爱心方程

MC冰河

前几天在一个网站上看到爱心方程17*x^2-16*abs(x)*y+17*y^2-225=0，matlab代码ezplot('17*x^2-16*abs(x)*y+17*y^2-225')，还有立体的，
1.[x,y,z]=meshgrid(linspace(-1.5,1.5));
2.f=@(x,y,z)(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3;
3.p=patch(isosurface(x,y,z,f(x,y,z),0)); view(3);axis equal
4.set(p,'FaceVertexCData',jet(size(get(p,'faces'),1)) ,'FaceColor', 'flat', 'EdgeColor', 'none');
我曾经也写过
theta1=(0.5:0.01:1.5)*pi;
rho1=2*theta1;
theta2=(1.5:0.01:2.5)*pi;
rho2=2*(3*pi-theta2);
polar(theta1,rho1,'r')
hold on
polar(theta2,rho2,'r')
也是个爱心

darker50

  建议给出一个结果图，谢谢分享！

绿色子

挺好的，效果也非常好、