Guest User

Untitled

a guest
Jul 18th, 2019
59
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. R.<x,y> = PolynomialRing(QQ,2)
  2. var('z')
  3. f1 = 432*x^5+144*x^4*y-504*x^3*y^2-168*x^2*y^3+99*x*y^4+33*y^5+36*x^4*z-13*x^2*y^2*z+y^4*z-540*x^3*z^2-180*x^2*y*z^2+216*x*y^2*z^2+72*y^3*z^2-45*x^2*z^3+108*x*z^4+36*y*z^4+9*z^5
  4. f = f1.subs(x=1,y=x,z=y)
  5. ff(x,y) = f
  6. mylinewidth = 4.5
  7. bb1 = (-4,6)
  8. bb2 = (-3,3)
  9.  
  10. G = Graphics()
  11. G += region_plot([ff>=0],(x,-1.2,1.2),(y,-1.2,1.2), ticks=[[],[]], axes=false, gridlines=False,incol='lightgreen')
  12. G += implicit_plot(ff, bb1, bb2, linewidth=mylinewidth, ticks=[[],[]], axes=false, gridlines=False, color="blue")
  13. G.axes(False)
  14. G.axes_label_color('white')
  15.  
  16. G.show(ticks=[[],[]], axes=False, gridlines=False) # dpi = 100
RAW Paste Data