My Geeky Valentine’s Day Gift

February 14, 2013

My Geeky Valentine’s Day Gift

Today, on 14 February 2013, I officially spread the love to the whole internet:

function mysteryFunc()

[x,y,z] = meshgrid(-1.5:0.025:1.5,-1.0:0.025:1.5,-1.0:0.025:1.5);
 rr = 320*((x.^2 + (9*y.^2)/4 + z.^2 - 1).^3 - x.^2.*z.^3 - (9*y.^2.*z.^3)/80);
 p=patch(isosurface(x,y,z,rr,0));
 isonormals(x,y,z,rr, p)
 set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
 daspect([1 1 1])
 view(3)
 camlight; lighting phong

end

Can you guess what this script is about ? To find out about it, please execute it in the Matlab. (note : unfortunately, this script does not work for Octave, I will check it later).

Disclaimer : I found the equation from somewhere in the internet, therefore I’m not the one who invented it.

Permalink