kill(all)$

fpprintprec:4$

r2d(a):=180*a/%pi$

load(atrig1)$
load(ntrig)$

van=facebook.com/LearningCalculus$

bc:(((sin((π*8/9)))^((-2))*((sin((π*8/9)))^(2)+
     (cos((π*35/18))*sin((π*8/9))*(-1))+
     (cos((π*35/18)))^(2))))^(1/2);

float(bc);

CD(PA):=((((((((sin((π*29/18))*PA*(3)^(1/2)*1/2)+(cos((π*29/18))*(-1))+
                (cos((π*29/18))*PA*1/2)))^((-1))*
             ((((cos((π*35/18))*(sin((π*8/9)))^((-1))*sin((π*11/18))*(3)^(1/2)*1/2)+
                (cos((π*35/18))*cos((π*11/18))*(sin((π*8/9)))^((-1))*1/2))*PA*(3)^(1/2)*(-1/2))+
                (cos((π*29/18))*PA*(3)^(1/2)*(-1/2)))*(-1))+
                (cos((π*35/18))*(sin((π*8/9)))^((-1))*(3)^(1/2)*1/2)))^(2)+
           ((((((sin((π*29/18))*PA*(3)^(1/2)*1/2)+(cos((π*29/18))*(-1))+
                (cos((π*29/18))*PA*1/2)))^((-1))*((((cos((π*35/18))*
                (sin((π*8/9)))^((-1))*sin((π*11/18))*(3)^(1/2)*1/2)+
                (cos((π*35/18))*cos((π*11/18))*(sin((π*8/9)))^((-1))*1/2))*(1+(PA*(-1/2))))+
                (sin((π*29/18))*PA*(3)^(1/2)*1/2))*(-1))+
                (cos((π*35/18))*(sin((π*8/9)))^((-1))*1/2)))^(2)))^(1/2)$

CD(PA) = PA + 1 + bc$
solve(%,PA)$
f(PA):=''(rhs(first(%)))$
f(x)$
subst(PI,%pi,%)$
bfloat(f(x));

newton(t):=bfloat(subst(t,x,x-f(x)/diff(f(x),x)))$
pa:newton(newton(newton(4.098877012128954b-1)));
pa:0.4098877012128954;

phi(PA):=atan((((2+(PA*(-1))))^((-1))*PA*(3)^(1/2)));

phi(pa);

float(r2d(phi(pa))) * degrees;