
(require (quote calc-ext))

(defun calc-Need-calc-funcs nil (byte-code "" [nil] 1))

(defun calc-inc-gamma (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "    # #,  ' #, #" [arg calc-is-inverse calc-is-hyperbolic calc-binary-op "gamG" calcFunc-gammaG "gamQ" calcFunc-gammaQ "gamg" calcFunc-gammag "gamP" calcFunc-gammaP] 10)) 117] 3))

(defun calc-erf (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "  # #" [arg calc-is-inverse calc-unary-op "erfc" calcFunc-erfc "erf" calcFunc-erf] 6)) 432] 3))

(defun calc-erfc (arg) (interactive "P") (byte-code " !" [arg nil calc-invert-func calc-erf] 3))

(defun calc-beta (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "#" [arg calc-binary-op "beta" calcFunc-beta] 4)) 747] 3))

(defun calc-inc-beta nil (interactive) (byte-code "\"" [nil calc-do (lambda nil (byte-code "  !B# !B#" [calc-is-hyperbolic calc-enter-result 3 "betB" calcFunc-betaB calc-top-list-n "betI" calcFunc-betaI] 9)) 911] 3))

(defun calc-bessel-J (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "#" [arg calc-binary-op "besJ" calcFunc-besJ] 4)) 1153] 3))

(defun calc-bessel-Y (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "#" [arg calc-binary-op "besY" calcFunc-besY] 4)) 1322] 3))

(defun calc-bernoulli-number (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "  # #" [arg calc-is-hyperbolic calc-binary-op "bern" calcFunc-bern calc-unary-op] 6)) 1491] 3))

(defun calc-euler-number (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "  # #" [arg calc-is-hyperbolic calc-binary-op "eulr" calcFunc-euler calc-unary-op] 6)) 1714] 3))

(defun calc-stirling-number (arg) (interactive "P") (byte-code "\"" [nil calc-do (lambda nil (byte-code "  # #" [arg calc-is-hyperbolic calc-binary-op "str2" calcFunc-stir2 "str1" calcFunc-stir1] 6)) 1934] 3))

(defun calc-utpb nil (interactive) (byte-code "\"" [nil calc-prob-dist "b" 3] 3))

(defun calc-utpc nil (interactive) (byte-code "\"" [nil calc-prob-dist "c" 2] 3))

(defun calc-utpf nil (interactive) (byte-code "\"" [nil calc-prob-dist "f" 3] 3))

(defun calc-utpn nil (interactive) (byte-code "\"" [nil calc-prob-dist "n" 3] 3))

(defun calc-utpp nil (interactive) (byte-code "\"" [nil calc-prob-dist "p" 2] 3))

(defun calc-utpt nil (interactive) (byte-code "\"" [nil calc-prob-dist "t" 2] 3))

(defun calc-prob-dist (letter nargs) (byte-code "\"" [calc-do (lambda nil (byte-code "  	P	P!!DS\"\"#4 	P	P!!DS\"\"#" [nargs letter calc-is-inverse calc-enter-result "ltp" append intern "calcFunc-ltp" calc-top-n 1 calc-top-list-n 2 "utp" "calcFunc-utp"] 16)) 2699] 3))

(defun calcFunc-gamma (x) (byte-code "!
 \"\"!" [x math-numberp math-reject-arg numberp calcFunc-fact math-add -1] 6))

(defun math-gammap1-raw (x &optional fprec nfprec) (byte-code " 	!	[!!\"N !
\"=  !\"
# \"!\"\"!K \"
#\") !X \"`   \"! \"\"\"!%!\"*" [fprec calc-internal-prec nfprec x xplus1 t xinv lnx math-float math-lessp-float calcFunc-re math-neg math-div math-pi math-mul math-gammap1-raw math-add (float -1 0) math-sin-raw (float 1 0) math-realp (float 736276 0) math-overflow 1 math-ln-raw math-sqrt-two-pi math-exp-raw math-gamma-series math-sub (float 5 -1) math-sqr (float 0 0) 2] 38))

(defun math-gamma-series (sum x xinvsqr oterm n) (byte-code "=
 	\"!
8!
8S\"\"!\"\"	\"	\"9 [ \"VL ![ \"\\%+" [calc-display-working-message sum bn n term x next calc-internal-prec xinvsqr lots math-do-working "gamma" math-bernoulli-number math-mul math-div-float math-float 1 * 2 math-add math-nearly-equal calc-record-why "*Gamma computation stopped early, not all digits may be valid" math-gamma-series] 20))

(defun calcFunc-gammaP (a x) (byte-code "\" ł]  	:? 	@> 	\"!* \"	!D 	!!D 	VD 	WQ 	\"\"] 	!	\"
\")" [x a math-current-gamma-value equal (var inf var-inf) (float 1 0) math-inexact-result (bigpos bigneg frac float cplx polar) math-reject-arg numberp math-numberp math-num-integerp integerp math-trunc 0 20 math-sub 1 calcFunc-gammaQ calcFunc-gamma math-div calcFunc-gammag] 15))

(defun calcFunc-gammaQ (a x) (byte-code "\" ʂ  	:? 	@> 	\"!* \"	!D 	!!D 	VD 	W \\
T	Wy \"
\"\"=u \"R !!\")!+ 	!	\"\")" [x a n sum term calc-internal-prec calc-display-working-message math-current-gamma-value equal (var inf var-inf) (float 0 0) math-inexact-result (bigpos bigneg frac float cplx polar) math-reject-arg numberp math-numberp math-num-integerp integerp math-trunc 0 20 (float 1 0) (float 1 0) math-normalize 1 math-div math-mul math-add lots math-do-working "gamma" calcFunc-exp math-neg calcFunc-gamma calcFunc-gammaG] 21))

(defun calcFunc-gammag (a x) (byte-code "\" 	!q  	:? 	@>! 	\":?, @>3 \"
\\	!!	!!X !	!\"\"b 	\"o j 	!	\"\")!" [x a calc-internal-prec math-current-gamma-value equal (var inf var-inf) calcFunc-gamma math-inexact-result (bigpos bigneg frac float cplx polar) math-reject-arg numberp (bigpos bigneg frac float cplx polar) math-normalize 2 math-float math-negp calcFunc-re math-lessp-float math-add-float (float 1 0) math-inc-gamma-series math-sub math-inc-gamma-cfrac] 21))

(setq math-current-gamma-value nil)

(defun calcFunc-gammaG (a x) (byte-code "\" Ƃo  	:? 	@> 	\":?* @>1 \"
\\	!!	!!V !	!\"\"i a 	!	\"\"m 	\")!" [x a calc-internal-prec math-current-gamma-value equal (var inf var-inf) (float 0 0) math-inexact-result (bigpos bigneg frac float cplx polar) math-reject-arg numberp (bigpos bigneg frac float cplx polar) math-normalize 2 math-float math-negp calcFunc-re math-lessp-float math-add-float math-abs-approx (float 1 0) math-sub calcFunc-gamma math-inc-gamma-series math-inc-gamma-cfrac] 19))

(defun math-inc-gamma-series (a x) (byte-code ":\" @>? @= 8= !% =, ɂL 	!\"\"!
\\	\"	$))!\"" [x a calc-internal-prec start (bigpos bigneg) float 1 0 math-zerop (float 0 0) math-mul math-exp-raw math-sub math-ln-raw math-normalize 2 math-div (float 1 0) math-inc-gamma-series-step] 14))

(defun math-inc-gamma-series-step (sum term a x) (byte-code "=
 	\"
\"\"
\"	\"	\"* 0 
$)" [calc-display-working-message sum a term x next lots math-do-working "gamma" math-add (float 1 0) math-div math-mul math-nearly-equal math-inc-gamma-series-step] 11))

(defun math-inc-gamma-cfrac (a x) (byte-code ":\" @>? @= 8= !% =2 	/ 
!K 
!\"\"!
&	\"" [x math-current-gamma-value a (bigpos bigneg) float 1 0 math-zerop calcFunc-gamma math-mul math-exp-raw math-sub math-ln-raw math-inc-gamma-cfrac-step (float 1 0) (float 0 0) (float 1 0) (float 1 0) (float 1 0) (float 0 0)] 18))

(defun math-inc-gamma-cfrac-step (a0 a1 b0 b1 n fac g a x) (byte-code "	
\"	\"	\"\"\"\"\"\"\"	\"\"\"	\"\"\"!c 	

	&	 \"\"=~ \"
\"  	
	&	)*" [ana n a anf fac a0 a1 b0 b1 x g next calc-display-working-message math-sub math-mul math-add (float 1 0) math-zerop math-inc-gamma-cfrac-step math-div (float 1 0) lots math-do-working "gamma" math-nearly-equal] 31))

(defun calcFunc-erf (x) (byte-code "\" ĂT \" ƂT :8 @>?5 @=2 8=5 !; =B T  !!\"	\"\")" [x math-current-gamma-value equal (var inf var-inf) (float 1 0) (neg (var inf var-inf)) (float -1 0) (bigpos bigneg) float 1 0 math-zerop math-sqrt-pi math-to-same-complex-quad math-div calcFunc-gammag (float 5 -1) math-sqr math-to-complex-quad-one] 11))

(defun calcFunc-erfc (x) (byte-code "\" Ă' !!  !\"	\")' !\"" [x math-current-gamma-value equal (var inf var-inf) (float 0 0) math-posp math-sqrt-pi math-div calcFunc-gammaG (float 5 -1) math-sqr math-sub 1 calcFunc-erf] 10))

(defun math-to-complex-quad-one (x) (byte-code "!= !!=% 8!8!E& " [x car-safe polar math-complex cplx math-abs 1 2] 9))

(defun math-to-same-complex-quad (x y) (byte-code "!=M 	!=; 8!! 	8!$ 	88!4 	8!7 	8EJ 8!I 	!J 	o !n 	!=h 	8!	8Ek 	!o 	" [y x car-safe cplx math-negp 1 math-neg 2] 15))

(defun calcFunc-beta (a b) (byte-code "!' \"
! 
\"

	\"	\"\"\")D 
!4 
\"D !
!\"
\"!\"" [a am b math-num-integerp math-add -1 math-numberp math-reject-arg numberp math-div 1 math-mul calcFunc-choose calcFunc-beta calcFunc-gamma] 19))

(defun calcFunc-betaI (x a b) (byte-code "!
 Ƃ \" ɂ 	!$ 	!$ 	!D 
!6 
!6 
!@ 
\"A ΂ 
!S 
!S 
!Z ς 	!?h 	\" 
!?v 
\"   	
\"	
#\")" [x a b t math-current-beta-value math-zerop (float 0 0) math-equal-int 1 (float 1 0) math-num-integerp math-negp math-reject-arg range (float 1 0) (float 0 0) math-numberp numberp math-inexact-result calcFunc-beta math-div calcFunc-betaB] 23))

(defun calcFunc-betaB (x a b) (byte-code "!
 ͂\" 	
\"!?& \"	!?4 	\"
!?B 
\"	!O 	\"
!\\ 
\"
! 
!o 
\" 
!: : 
!AA\"U : 
!W 	! 
!  \\	\"T
W 
Z\"\"\"	\"\"\"= \" +)!	\"\"	!B	!	\"B	!:6:3	!AA\"UB:B	!W[	O	
\"\"
	#\" \\!	!
!	!\"
\"!\"\"!	\"	
\"\"\"\"U	
#\"	\"		
\"
	\"#\"
\"\"))!" [x a b calc-symbolic-mode calc-internal-prec i term sum calc-display-working-message math-current-beta-value t bt math-zerop (float 0 0) math-equal-int 1 calcFunc-beta math-numberp math-reject-arg numberp nonzerop math-num-integerp math-negp range math-trunc 20 math-compare-bignum -1 math-floatp math-inexact-result math-mul math-normalize 2 0 math-div math-add lots math-do-working "beta" math-pow math-sub calcFunc-betaB math-float math-exp-raw math-ln-raw (float 1 0) math-compare (float 1 0) (float 2 0) math-beta-cfrac] 77))

(setq math-current-beta-value nil)

(defun math-beta-cfrac (a b x) (byte-code "	
\"	\"	\"\"\"\"	
&+" [qab a b qap qam x math-add (float 1 0) (float -1 0) math-beta-cfrac-step (float 1 0) math-sub (float 1 0) math-div math-mul (float 1 0) (float 1 0) (float 1 0)] 18))

(defun math-beta-cfrac-step (az bz am bm m qab qap qam a b x) (byte-code "	\"	\"	\"\"\"\"\"\"
	\"\"
\"\"
	\"	\"\"\"\"\"\"\"!\"\"
\"\"\"= \"\"  \"
\"	\"&." [two-m m d b x qam a ap az am bp bz bm d2 qab qap app bpp next calc-display-working-message math-mul (float 2 0) math-div math-sub math-add math-neg lots math-do-working "beta" math-nearly-equal math-beta-cfrac-step (float 1 0) (float 1 0)] 43))

(defun calcFunc-besJ (v x) (byte-code "!
 \"	! 	\"
^
\\	!!!!	!C !? ւ@ ׂ !?o !\"	!\"%	\"\"\")!! ! !	\"! !	\"= 	!= 	!	!\"U 	!	!	\"TW \"\"\" -\"U\"\"!\\\"	\"
SV\"\"\"!8Vr\"\"h\"\"
?
\"U%\"\"\".)!)" [v x calc-internal-prec start j bjm bj two-over-x bjp t jsum nil sum ans math-numberp math-reject-arg numberp 8 math-normalize 3 math-float math-zerop (float 1 0) (float 0 0) math-inexact-result math-num-integerp math-div 1 calcFunc-fact math-mul math-besJ-series 0 (float -25 -2) math-sqr math-pow 2 math-negp math-trunc math-oddp math-neg calcFunc-besJ math-besJ0 math-besJ1 math-compare math-abs-approx -1 math-sub 100 range logior math-isqrt-small * 40 (float 0 0) (float 0 0) (float 1 0) 10 (float 1 -10) (float 1 -10) (float 1 -10) (float 1 -10) math-add] 60))

(defun math-besJ-series (sum term k zz vk) (byte-code "=
 	\"
T\"\"
\"\"	\"	\"3 ; 
%)" [calc-display-working-message sum k vk term zz next lots math-do-working "besJ" math-add 1 math-div math-mul math-nearly-equal math-besJ-series] 13))

(defun math-besJ0 (x &optional yflag) (byte-code "?
 	!! 	!!} 	!\"Un 	\"
!	\"\"\"!N A!@B	\"!A\"@
\"\"\"\".} 	!\"\"\")" [yflag x z y xx a1 a2 sc t math-negp calcFunc-re math-besJ0 math-neg math-compare (float 8 0) math-abs-approx -1 math-div (float 8 0) math-sqr math-add (float (bigneg 164 398 785) -9) math-poly-eval ((float (bigpos 211 887 93 2) -16) (float (bigneg 639 370 73 2) -15) (float (bigpos 407 510 734 2) -14) (float (bigneg 627 628 98 1) -12) (float 1 0)) ((float (bigneg 152 935 934) -16) (float (bigpos 161 95 621 7) -16) (float (bigneg 651 147 911 6) -15) (float (bigpos 765 488 430 1) -13) (float (bigneg 995 499 562 1) -11)) math-sin-cos-raw math-mul math-sqrt (float (bigpos 722 619 636) -9) math-sub ((float (bigneg 456 52 849 1) -7) (float (bigpos 17 233 739 7) -5) (float (bigneg 418 442 121 1) -2) (float (bigpos 407 196 516 6) -1) (float (bigneg 354 590 362 13) 0) (float (bigpos 574 490 568 57) 0)) ((float 1 0) (float (bigpos 712 532 678 2) -7) (float (bigpos 853 264 927 5) -5) (float (bigpos 718 680 494 9) -3) (float (bigpos 985 532 29 1) 0) (float (bigpos 411 490 568 57) 0))] 27))

(defun math-besJ1 (x &optional yflag) (byte-code "!!
 	? !!! !\"U \"
!\"\"\"!	S A!@Bg !g @!A!B\"!A\"@
\"\"\"\". !\"\"\"\")" [x yflag z y xx a1 a2 sc t math-negp calcFunc-re math-neg math-besJ1 math-compare (float 8 0) math-abs-approx -1 math-div (float 8 0) math-sqr math-add (float (bigneg 491 194 356 2) -9) math-poly-eval ((float (bigneg 19 337 240) -15) (float (bigpos 174 520 457 2) -15) (float (bigneg 496 396 516 3) -14) (float 183105 -8) (float 1 0)) ((float (bigpos 412 787 105) -15) (float (bigneg 987 228 88) -14) (float (bigpos 96 199 449 8) -15) (float (bigneg 873 690 2 2) -13) (float (bigpos 995 499 687 4) -11)) math-sin-cos-raw math-mul math-sqrt (float (bigpos 722 619 636) -9) math-sub ((float (bigneg 606 36 16 3) -8) (float (bigpos 826 44 157) -4) (float (bigneg 439 611 972 2) -3) (float (bigpos 531 968 423 2) -1) (float (bigneg 235 59 895 7) 0) (float (bigpos 232 614 362 72) 0)) ((float 1 0) (float (bigpos 397 991 769 3) -7) (float (bigpos 394 743 944 9) -5) (float (bigpos 474 330 858 1) -2) (float (bigpos 178 535 300 2) 0) (float (bigpos 442 228 725 144) 0))] 33))

(defun calcFunc-besY (v x) (byte-code " ! \"	! 	\"
^
\\	!!!!!?Y  \"!	\"A\"!	\"\"@\") !!} !t !	\"!z !	\" = 	! = 	! 	!	!	\"	TW \"\"\"		 -)!)" [v x calc-internal-prec sc t j bym by two-over-x byp math-inexact-result math-numberp math-reject-arg numberp 8 math-normalize 3 math-float math-num-integerp math-sin-cos-raw math-mul math-pi math-div math-sub calcFunc-besJ math-neg math-negp math-trunc math-oddp calcFunc-besY 0 math-besY0 1 math-besY1 2 nil] 39))

(defun math-besY0 (x) (byte-code "!\"U* !	\"	\"\"!!\"\"\")I !!E !\"!!\"\"I \"" [x y t math-compare math-abs-approx (float 8 0) -1 math-sqr math-add math-div math-poly-eval ((float (bigpos 733 622 284 2) -7) (float (bigneg 757 792 632 8) -5) (float (bigpos 129 988 87 1) -2) (float (bigneg 36 598 123 5) -1) (float (bigpos 65 834 62 7) 0) (float (bigneg 389 821 957 2) 0)) ((float 1 0) (float (bigpos 244 30 261 2) -7) (float (bigpos 647 472 474) -4) (float (bigpos 438 466 189 7) -3) (float (bigpos 648 499 452 7) -1) (float (bigpos 269 544 76 40) 0)) math-mul (float (bigpos 772 619 636) -9) math-besJ0 math-ln-raw math-negp calcFunc-re math-neg (cplx 0 2)] 22))

(defun math-besY1 (x) (byte-code "!\"U3 !	\"	\"\"\"!!\"\"\"\"\")T !!P !\"!!\"\"!T \"" [x y t math-compare math-abs-approx (float 8 0) -1 math-sqr math-add math-mul math-div math-poly-eval ((float (bigpos 935 937 511 8) -6) (float (bigneg 726 922 237 4) -3) (float (bigpos 551 264 349 7) -1) (float (bigneg 139 438 153 5) 1) (float (bigpos 439 527 127) 4) (float (bigneg 943 604 900 4) 3)) ((float 1 0) (float (bigpos 885 632 549 3) -7) (float (bigpos 605 42 102) -3) (float (bigpos 2 904 245 2) -2) (float (bigpos 367 650 733 3) 0) (float (bigpos 664 419 244 4) 2) (float (bigpos 57 958 249) 5)) (float (bigpos 772 619 636) -9) math-sub math-besJ1 math-ln-raw 1 math-negp calcFunc-re math-neg (cplx 0 2)] 26))

(defun math-poly-eval (x coefs) (byte-code "	@	A 	@
\"\" )" [accum coefs x math-add math-mul] 6))

(defun calcFunc-bern (n &optional x) (byte-code " !?# 	 !    
!\"C 
!- 
\"
:@  
!!!C 
!" [x calc-symbolic-mode n math-zerop math-floatp math-inexact-result math-build-polynomial-expr math-bernoulli-coefs math-num-natnump math-reject-arg natnump math-float math-bernoulli-number math-trunc] 13))

(defun calcFunc-euler (n &optional x) (byte-code "!
 \"	j \"
!
\"
\"	\"\"	\"	!Y ? 	!G  V \"\"\"\"e \"B\"	\". \":~  \" \"\"" [n x n1 coefs fac k x1 x2 calc-symbolic-mode math-num-natnump math-reject-arg natnump math-add 1 math-bernoulli-coefs math-div math-pow 2 -1 math-numberp math-floatp math-inexact-result math-mul math-sub math-build-polynomial-expr calcFunc-collect math-reduce-vec vec mapcar (lambda (c) (byte-code "T	
\"\"\"\"\"" [k fac c x1 x2 math-mul math-sub math-pow] 9)) calcFunc-euler (float 5 -1) (frac 1 2)] 27))

(defun math-bernoulli-coefs (n) (byte-code "	!C	!

	! SYI 
Z\"!\"	:; !< B\" !." [coefs n nn k term coef calc-prefer-frac calcFunc-bern math-trunc nil integerp 0 math-div math-mul math-bernoulli-number math-float nreverse] 11))

(defun math-bernoulli-number (n) (byte-code "\"U U ς Ђ \"	Y 
=A 	\"\"j \\S\"\"
@\"\"AB S\"\"
\"
B
\"B	T. 	#8" [n math-bernoulli-cache-size sum nk fact ofact p math-bernoulli-b-cache calc-prefer-frac t calc-display-working-message math-bernoulli-B-cache % 2 1 (frac -1 2) 0 / nil lots math-do-working "bernoulli B" * math-mul math-add math-div math-sub (frac 1 2) -] 16))

(setq math-bernoulli-b-cache (quote ((frac -174611 (bigpos 0 200 291 698 662 857 802)) (frac 43867 (bigpos 0 944 170 217 94 109 5)) (frac -3617 (bigpos 0 880 842 622 670 10)) (frac 1 (bigpos 600 249 724 74)) (frac -691 (bigpos 0 368 674 307 1)) (frac 1 (bigpos 160 900 47)) (frac -1 (bigpos 600 209 1)) (frac 1 30240) (frac -1 720) (frac 1 12) 1)))

(setq math-bernoulli-B-cache (quote ((frac -174611 330) (frac 43867 798) (frac -3617 510) (frac 7 6) (frac -691 2730) (frac 5 66) (frac -1 30) (frac 1 42) (frac -1 30) (frac 1 6) 1)))

(setq math-bernoulli-cache-size 11)

(defun calcFunc-utpb (x n p) (byte-code " 	

EEF! 	

\"\"#" [math-expand-formulas p x n math-normalize calcFunc-betaI + - 1 math-add math-sub] 8))

(put (quote calcFunc-utpb) (quote math-expandable) t)

(defun calcFunc-ltpb (x n p) (byte-code "	
#\"" [x n p math-sub 1 calcFunc-utpb] 6))

(put (quote calcFunc-ltpb) (quote math-expandable) t)

(defun calcFunc-utpc (chisq v) (byte-code " 	E
EE! 	\"
\"\"" [math-expand-formulas v chisq math-normalize calcFunc-gammaQ / 2 math-div] 7))

(put (quote calcFunc-utpc) (quote math-expandable) t)

(defun calcFunc-ltpc (chisq v) (byte-code " 	E
EE! 	\"
\"\"" [math-expand-formulas v chisq math-normalize calcFunc-gammaP / 2 math-div] 7))

(put (quote calcFunc-ltpc) (quote math-expandable) t)

(defun calcFunc-utpf (f v1 v2) (byte-code " 		
EEE	E
EF!1 		
\"\"\"	\"
\"#" [math-expand-formulas v2 v1 f math-normalize calcFunc-betaI / + * 2 math-div math-add math-mul] 11))

(put (quote calcFunc-utpf) (quote math-expandable) t)

(defun calcFunc-ltpf (f v1 v2) (byte-code "	
#\"" [f v1 v2 math-sub 1 calcFunc-utpf] 6))

(put (quote calcFunc-ltpf) (quote math-expandable) t)

(defun calcFunc-utpn (x mean sdev) (byte-code " 	
EDEEDEE!0 	
\" \"\"!\"\"" [math-expand-formulas mean x sdev math-normalize / + 1 calcFunc-erf - * calcFunc-sqrt 2 math-mul math-add (float 1 0) math-div math-sub math-sqrt-2 (float 5 -1)] 12))

(put (quote calcFunc-utpn) (quote math-expandable) t)

(defun calcFunc-ltpn (x mean sdev) (byte-code " 	
EDEEDEE!0 	
\" \"\"!\"\"" [math-expand-formulas x mean sdev math-normalize / + 1 calcFunc-erf - * calcFunc-sqrt 2 math-mul math-add (float 1 0) math-div math-sub math-sqrt-2 (float 5 -1)] 12))

(put (quote calcFunc-ltpn) (quote math-expandable) t)

(defun calcFunc-utpp (n x) (byte-code " 	
E! 	
\"" [math-expand-formulas x n math-normalize calcFunc-gammaP] 4))

(put (quote calcFunc-utpp) (quote math-expandable) t)

(defun calcFunc-ltpp (n x) (byte-code " 	
E! 	
\"" [math-expand-formulas x n math-normalize calcFunc-gammaQ] 4))

(put (quote calcFunc-ltpp) (quote math-expandable) t)

(defun calcFunc-utpt (tt v) (byte-code " 		
EEE	EF!* 		
!\"\"	\"#" [math-expand-formulas v tt math-normalize calcFunc-betaI / + ^ 2 (float 5 -1) math-div math-add math-sqr (float 5 -1)] 9))

(put (quote calcFunc-utpt) (quote math-expandable) t)

(defun calcFunc-ltpt (tt v) (byte-code "	\"\"" [tt v math-sub 1 calcFunc-utpt] 5))

(put (quote calcFunc-ltpt) (quote math-expandable) t)
