1(* 2 Copyright (C) 2017 M.A.L. Marques 3 4 This Source Code Form is subject to the terms of the Mozilla Public 5 License, v. 2.0. If a copy of the MPL was not distributed with this 6 file, You can obtain one at http://mozilla.org/MPL/2.0/. 7*) 8 9(* type: gga_exc *) 10 11lag_a1 := 0.041106: 12lag_a2 := 2.626712: 13lag_a3 := 0.092070: 14lag_a4 := 0.657946: 15 16lag_f0 := s-> lag_a1 * s^lag_a2/(1 + lag_a3 * s^lag_a2)^lag_a4: 17lag_f := x-> lag_f0(X2S*x): 18 19f := (rs, zeta, xt, xs0, xs1) -> gga_exchange(lag_f, rs, zeta, xs0, xs1): 20