import numpy as np
import matplotlib.pyplot as plt
 
def plot_circle(center=(3, 3),r=2):
  x = np.linspace(center[0] - r, center[0] + r, 5000)
  y1 = np.sqrt(r**2 - (x-center[0])**2) + center[1]
  y2 = -np.sqrt(r**2 - (x-center[0])**2) + center[1]
     
  plt.plot(x, y1, c='k')
  plt.plot(x, y2, c='k')
  plt.show()

import matplotlib.pyplot as plt
plt.xlim(0, 15)
plt.ylim(0, 15)
plot_circle((5, 5),r=3)