import numpy as np
import matplotlib.pyplot as plt
plt.figure(figsize=(5,5),dpi=80)
N=1001
a=-1
b=1
t=(b-a)/(N-1)
x = np.arange(a,b+t,t)
y = np.arange(a,b+t,t)
X1=[]
Y1=[]
X2=[]
Y2=[]
for i in range(0,N,1):
for j in range(0,N,1):
if abs(x[i]*x[i]+y[j]*y[j] -1.0) <= 0.01:
X1.append(x[i])
Y1.append(y[j])
if abs(abs(abs(x[i])-0.3)+abs(abs(y[j])-0.3)-0.6) <= 0.001:
X2.append(x[i])
Y2.append(y[j])
plt.xlim(a*1.25, b*1.25)
plt.ylim(a*1.25, b*1.25)
plt.scatter(X1,Y1,color="red")
plt.scatter(X2,Y2,color="blue")
plt.show()
