Test of functions in package WedgeMOI

Load package

Needs["WedgeMOI`"]

aLIGO-ish numerical values for testing

vals= {
rho->2201,
R->0.17,
t->0.2,
f->0.095,
wabh->0,
wabv->-1 Pi/180, (* 1 degree top-heavy vertical wedge on rear face *)
wafh->0,
wafv->0
};

masswedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

masswedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

xCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]

xCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

yCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]

yCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

zCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]

zCOMwedge[R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

xMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

xMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

zMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

zMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

xyMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

xyMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

yzMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

yzMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals

zxMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]

zxMOIwedge[rho,R,t,f,{wabh,wafh},{wabv,wafv}]/.vals