1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
(*
Copyright (C) 2017 M.A.L. Marques
This Source Code Form is subject to the terms of the Mozilla Public
License, v. 2.0. If a copy of the MPL was not distributed with this
file, You can obtain one at http://mozilla.org/MPL/2.0/.
*)
(* Equation (7) from the paper *)
tpss_p := x -> X2S^2*x^2:
tpss_z := (x, t) -> x^2/(8*t):
tpss_alpha := (x, t) -> (t - x^2/8)/K_FACTOR_C:
tpss_qb := (x, t) -> \
9/20 * (tpss_alpha(x, t) - 1)/sqrt(1 + params_a_b*tpss_alpha(x, t)*(tpss_alpha(x, t) - 1)) \
+ 2*tpss_p(x)/3:
(* Equation (10) in all its glory *)
tpss_fxnum := (x, t) -> \
+ (MU_GE + params_a_c*tpss_z(x, t)^tpss_ff(tpss_z(x, t))/(1 + tpss_z(x, t)^2)^2)*tpss_p(x) \
+ 146/2025 * tpss_qb(x, t)^2 \
- 73/405 * tpss_qb(x, t) * sqrt(1/2*(9/25*tpss_z(x, t)^2 + tpss_p(x)^2)) \
+ MU_GE^2/tpss_kappa(x, t) * tpss_p(x)^2 \
+ 2*sqrt(params_a_e)*MU_GE*9/25*tpss_z(x, t)^2 \
+ params_a_e*params_a_mu*tpss_p(x)^3:
tpss_fxden := x -> \
(1 + sqrt(params_a_e)*tpss_p(x))^2:
tpss_fx := (x, t) -> tpss_fxnum(x, t)/tpss_fxden(x):
|
