-
Notifications
You must be signed in to change notification settings - Fork 1
/
batman1.py
23 lines (18 loc) · 832 Bytes
/
batman1.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
from __future__ import division
from numpy import sqrt
from numpy import meshgrid as mg
from numpy import arange
import matplotlib.pyplot as plot
x1 = arange(-7.25, 7.25, 0.01)
y1 = arange(-5, 5, 0.01)
x, y = mg(x1, y1)
eq1 = ((x/7)**2*sqrt(abs(abs(x)-3)/(abs(x)-3))+(y/3)**2*sqrt(abs(y+3/7*sqrt(33))/(y+3/7*sqrt(33)))-1)
eq2 = (abs(x/2)-((3*sqrt(33)-7)/112)*x**2-3+sqrt(1-(abs(abs(x)-2)-1)**2)-y)
eq3 = (9*sqrt(abs((abs(x)-1)*(abs(x)-.75))/((1-abs(x))*(abs(x)-.75)))-8*abs(x)-y)
eq4 = (3*abs(x)+.75*sqrt(abs((abs(x)-.75)*(abs(x)-.5))/((.75-abs(x))*(abs(x)-.5)))-y)
eq5 = (2.25*sqrt(abs((x-.5)*(x+.5))/((.5-x)*(.5+x)))-y)
eq6 = (6*sqrt(10)/7+(1.5-.5*abs(x))*sqrt(abs(abs(x)-1)/(abs(x)-1))-(6*sqrt(10)/14)*sqrt(4-(abs(x)-1)**2)-y)
equation=[eq1,eq2,eq3,eq4,eq5,eq6]
for f in equation:
plot.contour(x, y, f, [0])
plot.show()