/*	theoretical target - binomial series  /%
	EIods(POIo,t,4th stage)
	=	+	E0pds(POIp)						*{		3/2*β^2*sin(Aθpc(POIo(t),t=0))^2 + 3/8*β^4*sin(Aθpc(POIo(t),t=0))^4 + -1/16*β^6*sin(Aθpc(POIo(t),t=0))^6 + 3/128*β^8*sin(Aθpc(POIo(t),t=0))^8}
		-	E0pds(POIp)*λ(Vons(PART))	*{1 + 3/2*β^2*sin(Aθpc(POIo(t),t=0))^2 + 3/8*β^4*sin(Aθpc(POIo(t),t=0))^4 + -1/16*β^6*sin(Aθpc(POIo(t),t=0))^6}

(mathL)/*	theoretical target - binomial series 
	n for (β*sin)^n					  0  2     4     6      8       10     12
	factors : 
	binomial series (John Wallis)	= 1  3/2   3/8  -1/16   3/128  -3/256  7/1024
	Lucas									= 1  3/2  15/8  35/16   
	12Oct2019 Howell - no (cos - 1) term, iterative, non-feedback, based on EIods(POIo,t,4th stage)
		ETpds(POIp)						= 1  3/2  21/8  35/8   455/64
		ETpds(POIp)*λ(Vons(PART))	= 1  1	  5/4	  5/3
(endMath)