function [V] = heart(x,y)

ix = (y'.*0) + 1;
iy = (x.*0) + 1;

V = ((y'.^4-20.*y'.^3)*iy./100) + (ix*(x.^2));