import numpy as np
import scipy.optimize
from scipy.optimize import fsolve
import math




def cycloid(a,p) :
	x,y = p
	func = a*np.arccos(1-y/a)-sqrt(y*(2*a-y)) - x
	return func



a0 = fsolve(cycloid , 1.0 , (6.0 , 0.0))

print a0