Tuesday, February 15, 2011

Sagemath valentine plot


only a day late ;)


def f(x,y):
return x^2 + (y-(x^2)^(1/3))^2

P2=contour_plot(f,(-1,1),(-1,1.5),cmap='Reds_r')

P2.show()

No comments: