p=2;
q=1.2;
ContourPlot[(Abs[x]/3)^p+(Abs[y]/1)^q==1,{x,-3,3},{y,-3,3}]

p=3;
q=3;
ContourPlot[(Abs[x]/3)^p+(Abs[y]/1)^q==1,{x,-3,3},{y,-3,3}]

c=2;
Water drop= ContourPlot [{((Sqrt[(1-x^2)])+y)((Sqrt[(1-x^2)])+1/(50x^2+c)-y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

a=0.6;
b=1;
c=2;
VAM1=ContourPlot [{((Sqrt[(1-x^2)])+b y)((Sqrt[(1-x^2)])+a/(50x^2+c)-b y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

a=0.3;
b=1.1;
c=1.6;
AGL1=ContourPlot [{((Sqrt[(1-x^2)])+b y)((Sqrt[(1-x^2)])+a/(50x^2+c)-b y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

a=3;
b=1;
c=5;
AAR1=ContourPlot [{((Sqrt[(1-x^2)])+b y)((Sqrt[(1-x^2)])+a/(50x^2+c)-b y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

c=1;
a=b=4;
Heart curve=ContourPlot [{(Sqrt[1-x^2]-4/(a(54x^2+9 Abs [x]+3))+c y)((Sqrt[1-x^2])+4/(b(54x^2+9 Abs [x]+3))-c y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

c=1;
a=8;
b=12;
ACO1=ContourPlot [{(Sqrt[1-x^2]-4/(a(54x^2+9 Abs [x]+3))+c y)((Sqrt[1-x^2])+4/(b(54x^2+9 Abs [x]+3))-c y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

c=1;
a=1/3;
b=1/2;
PHI1=ContourPlot [{(Sqrt[1-x^2]+a/(10x^2+1)-a/21-c y)((Sqrt[1-x^2])-a/(b(54x^2+9 Abs [x]+3))+c y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

c=0.9;
a=0.35;
b=0.3;
PPE1=ContourPlot [{(Sqrt[1-x^2]+a/(5x^4+25x^2+1)-c y)((Sqrt[1-x^2])-b/(5x^4+25Abs [x^3]+1)+c y)==0},{x,-2,2},{y,-2,2}, ContourStyle->Black]

a=0.9;(* aumentar a aumenta el pico superior*)
b=0.3; (* aumentar b profundiza el entrante inferior*)
ARO1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+5Abs[x^4]+25Abs[x]^2))-1.1y)(Sqrt[(1-x^2)]-b/( (1+5Abs[x^4]+25Abs[x]^3))+1.1y)==0,{x,-2,3},{y,-2,3},ColorFunction->Function[{x,y},{}],  ContourStyle->Black]

AJA1=ContourPlot [{(((25Sqrt[(36-x^2)])/22)+(25/(550(x^2/10)+15))-0.8y)(((25Sqrt[(36-x^2)])/25)-(5/(25(x^2/40)+2))+0.8y)==0},{x,-10,10},{y,-10,12},  ContourStyle->Black]

AER1=ContourPlot [{(((25Sqrt[(36-x^2)])/22)+(25/(550(x^2/10)+15))-1.14y)(((25Sqrt[(36-x^2)])/25)-(5/(25(x^2/40)+2))+1.14y)==0},{x,-10,10},{y,-10,12}, ContourStyle->Black]

a=0.3;(* aumentar a aumenta el pico superior*)
b=0.3; (* aumentar b profundiza el entrante inferior*)
AGR1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+5Abs[x^4]+25Abs[x]^2))-1.1y)(Sqrt[(1-x^2)]-a/( b(3+9Abs[x]+54x^2))+1.1y)==0,{x,-2,3},{y,-2,3},ColorFunction->Function[{x,y},{}], ContourStyle->Black]

a=0.40;(* aumentar a aumenta el pico superior*)
b=0.25; (* aumentar b profundiza el entrante inferior*)
ADE1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+10Abs[x^2]))-a/21-1.2y)(Sqrt[(1-x^2)]-b/( (1+10Abs[x^4]+10Abs[x]^3))+b/21+1.2y)==0,{x,-2,2},{y,-2,2},ColorFunction->Function[{x,y},{}],ContourStyle->Black]

AMA1=ContourPlot [{(((25Sqrt[(36-x^2)])/22)+(25/(550(x^2/10)+15))-0.8y)(((25Sqrt[(36-x^2)])/25)-(5/(25(x^2/40)+3))+0.8y)==0},{x,-10,10},{y,-10,12},  ContourStyle->Black]

a=1;(* aumentar a aumenta el pico superior*)
b=1/10; (* aumentar b profundiza el entrante inferior*)
YAU1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+3Abs [x]^1.6))-4a/21-0.75y)(Sqrt[(1-x^2)]-b/( (1+10Abs[x^4]+10Abs[x]^3))+b/21+0.75y)==0,{x,-2,2},{y,-2,3},ColorFunction->Function[{x,y},{}], ContourStyle->Black]

a=6;(* aumentar a rebaja el pico superior*)
b=8; 
c=1.5(* aumentar b suaviza la entrada inferior*)
COL1=ContourPlot[{(Sqrt[(1-(c x)^2)]-4/(a (3+54(c x)^2+9Abs[c x]))+1.2y)(Sqrt[(1-(c x)^2)]+4/(b (3+54(c x)^2+9Abs[c x]))-1.2y)==0},{x,-1,1.5},{y,-1,1.5},AspectRatio->1, ContourStyle->Black]

a=1/2;(* aumentar a aumenta el pico superior*)
b=1/10; (* aumentar b profundiza el entrante inferior*)
VAE1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+10Abs[x^2]))-a/21-0.85y)(Sqrt[(1-x^2)]-b/( (1+10Abs[x^4]+10Abs[x]^3))+b/21+0.85y)==0,{x,-2,2},{y,-2,2},ColorFunction->Function[{x,y},{}],  ContourStyle->Black]

a=0.6;(* aumentar a aumenta el pico superior*)
b=1/10; (* aumentar b profundiza el entrante inferior*)
AME1=ContourPlot[(Sqrt[(1-x^2)]+a/( (1+10Abs[x^1.4]))-a/21-0.80y)(Sqrt[(1-x^2)]-b/( (1+10Abs[x^4]+10Abs[x]^3))+b/21+0.80y)==0,{x,-2,2},{y,-2,2},ColorFunction->Function[{x,y},{}],  ContourStyle->Black]

a=6;(* aumentar a rebaja la entrada inferios*)
b=6; (* aumentar b rebaja el pico superior*)ESHC1=ContourPlot[(Sqrt[(1-x^2)]-4/(a (3+54x^2+9Abs[x]))+0.6y)(Sqrt[(1-x^2)]+4/(b(3+54x^2+9Abs[x]))-0.6y)==0,{x,-2,3},{y,-2,3}]

a=6;(* aumentar a rebaja el pico superior*)
EHSC2=ContourPlot[{(Sqrt[(1-x^2)]+0.6y)(Sqrt[(1-x^2)]+4/(a (3+54x^2+9Abs[x]))-0.6y)==0},{x,-2,3},{y,-2,3},AspectRatio->1,PlotPoints->10]