Grandfather Clock

Any ideas how to make the pendulum move more realistically? Currently it moves by the formula:

angular_speed=arctan(sqrt(maximum_height-current_height)))

and the current height is calculated by the formula:

current_height=lenght-length*cos(current_angle)

I thought those were the physically correct formulas, yet somehow it doesn't feel quite realistic.