[VA] SRC #009 - Pi Day 2021 Special

[Albert Chan]

For (a), (b), it is simpler to setup as iterative formula, X=F(X), instead of searching in range [3, pi]

>X=3 @ P=1E-10 @ W1=.56714329041
>INTEG(0,X,P,(SIN(IX)/IX*EXP(IX/TAN(IX)))^X) * GAMMA(X)/X^(X-1)
3.1415926536
>INTEG(0,X,P,LOGP1(SIN(IX)/IX*EXP(IX/TAN(IX)))) / W1
3.14159265358

Both (a),(b) already converged to X=PI (confirmed by X=PI, and run integral again)