var('x,y')
cm = colormaps.Blues
def f(x,y) :
return (2*x^2 + 3*x*y+4*y^2)/(3*x^2 + 5*y^2)
def c(x,y):
return f(x,y)
S=plot3d(f,(x,-2,2),(y,-2,2),color = (c,cm), opacity=1, mesh=1)
show(S, aspect_ratio = [1,1,3])
C=contour_plot(f, (x,-2, 2), (y,-2, 2), cmap='Blues',linestyles='solid', colorbar=True)
show(C,figsize=4)