clear;close all;
r = 5;
g = 9.8;
thetaC = pi/6;
thetaB = pi;
xb = r*(thetaB-sin(thetaB));
yb = r*(1-cos(thetaB));
xc = r*(thetaC-sin(thetaC));
yc = r*(1-cos(thetaC));
syms x
a=(1-cos(x))/(cos(thetaC)-cos(x));
a=matlabFunction(a);
a=quadgk(a,thetaC,thetaB)
t = sqrt(r/g)*a