================================================================

L 20 d 3 gru alpha 1.70 n2a 0.0  :#GRU alpha = 1.7000 <lmax>/|X| = 0.0072 phi(alpha) = 0.0012 quot = 6.1803
L 20 d 3 gru alpha 1.75 n2a 0.0  :#GRU alpha = 1.7500 <lmax>/|X| = 0.0030 phi(alpha) = 0.0000 quot = inf

L 20 d 3 gru alpha 1.72 n2a 0.01 :#GRU alpha = 1.7200 <lmax>/|X| = 0.0120 phi(alpha) = 0.0112 quot = 1.0692
L 20 d 3 gru alpha 1.73 n2a 0.01 :#GRU alpha = 1.7300 <lmax>/|X| = 0.0036 phi(alpha) = 0.0000 quot = inf
L 20 d 3 gru alpha 1.74 n2a 0.01 :#GRU alpha = 1.7400 <lmax>/|X| = 0.0080 phi(alpha) = 0.0012 quot = 6.4698
L 20 d 3 gru alpha 1.75 n2a 0.01 :#GRU alpha = 1.7500 <lmax>/|X| = 0.0027 phi(alpha) = 0.0000 quot = inf

L 20 d 3 gru alpha 1.76 n2a 0.05 :#GRU alpha = 1.7600 <lmax>/|X| = 0.0178 phi(alpha) = 0.0213 quot = 0.8322
L 20 d 3 gru alpha 1.77 n2a 0.05 :#GRU alpha = 1.7700 <lmax>/|X| = 0.0040 phi(alpha) = 0.0000 quot = inf

L 20 d 3 gru alpha 1.76 n2a 0.10 :#GRU alpha = 1.7600 <lmax>/|X| = 0.0054 phi(alpha) = 0.0003 quot = 17.2083
L 20 d 3 gru alpha 1.80 n2a 0.10 :#GRU alpha = 1.8000 <lmax>/|X| = 0.0040 phi(alpha) = 0.0000 quot = inf

L 20 d 3 gru alpha 1.76 n2a 0.20 :#GRU alpha = 1.7600 <lmax>/|X| = 0.0190 phi(alpha) = 0.0190 quot = 1.0000
L 20 d 3 gru alpha 1.80 n2a 0.20 :#GRU alpha = 1.8000 <lmax>/|X| = 0.0044 phi(alpha) = 0.0000 quot = inf



----------------------------------------------------------------
GRU alpha_c = 1.74 1.75 1.76 1.76 1.76

----------------------------------------------------------------
beta_c = 
0.143678
0.142857
0.142045
0.142045
0.142045

----------------------------------------------------------------
T_c =
6.960008
7.000007
7.040023
7.040023
7.040023

----------------------------------------------------------------
Tc0=6.96
for Tca in 7.000007 7.040023 7.040023 7.040023
do
	echo "($Tca - $Tc0)/$Tc0" | kbc -np
done

0.005748 / 0.01
0.011498 / 0.05
0.011498 / 0.10
0.011498 / 0.20

0.574800
0.229960
0.114980
0.057490


================================================================

L 40 d 3 gru alpha 1.692 n2a 0.00 :#GRU alpha = 1.6920 <lmax>/|X| = 0.0098 phi(alpha) = 0.0058 quot = 1.7060
L 40 d 3 gru alpha 1.701 n2a 0.01 :#GRU alpha = 1.7010 <lmax>/|X| = 0.0047 phi(alpha) = 0.0002 quot = 19.6245
L 40 d 3 gru alpha 1.741 n2a 0.05 :#GRU alpha = 1.7410 <lmax>/|X| = 0.0055 phi(alpha) = 0.0002 quot = 34.6528
L 40 d 3 gru alpha 1.767 n2a 0.10 :#GRU alpha = 1.7670 <lmax>/|X| = 0.0040 phi(alpha) = 0.0001 quot = 27.3321
L 40 d 3 gru alpha 1.778 n2a 0.20 :#GRU alpha = 1.7780 <lmax>/|X| = 0.0037 phi(alpha) = 0.0000 quot = 1883.6206




----------------------------------------------------------------
GRU alpha_c = 1.692 1.701 1.741 1.767 1.778

----------------------------------------------------------------
beta_c = 
0.147754
0.146972
0.143596
0.141483
0.140607

----------------------------------------------------------------
T_c = 
6.768006
6.804017
6.963982
7.067987
7.112021

----------------------------------------------------------------
Tc0=6.768006
for Tca in 6.804017 6.963982 7.067987 7.112021
do
	echo "($Tca - $Tc0)/$Tc0" | kbc -np
done

0.005321
0.028956
0.044323
0.050830

0.005321 / 0.01
0.028956 / 0.05
0.044323 / 0.10
0.050830 / 0.20

0.532100
0.579120
0.443230
0.254150

================================================================
L 50 d 3 gru alpha 1.693 n2a 0.00  :#GRU alpha = 1.6930 <lmax>/|X| = 0.0022 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.693 n2a 0.001 :#GRU alpha = 1.6930 <lmax>/|X| = 0.0031 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.705 n2a 0.01  :#GRU alpha = 1.7050 <lmax>/|X| = 0.0035 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.702 n2a 0.011 :#GRU alpha = 1.7020 <lmax>/|X| = 0.0034 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.74  n2a 0.05  :#GRU alpha = 1.7400 <lmax>/|X| = 0.0024 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.743 n2a 0.051 :#GRU alpha = 1.7430 <lmax>/|X| = 0.0025 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.78  n2a 0.10  :#GRU alpha = 1.7800 <lmax>/|X| = 0.0031 phi(alpha) = 0.0000 quot = inf
L 50 d 3 gru alpha 1.8   n2a 0.20  :#GRU alpha = 1.8000 <lmax>/|X| = 0.0012 phi(alpha) = 0.0000 quot = inf

L=50 d=3

alpha_gru n2a

1.693 0.000
1.693 0.001
1.705 0.010
1.702 0.011
1.74  0.050
1.743 0.051
1.78  0.100
1.8   0.200

----------------------------------------------------------------
beta_c

0.147667
0.147667
0.146628
0.146886
0.143678
0.143431
0.140449
0.138889

----------------------------------------------------------------
T_c

6.771994
6.771994
6.819980
6.808001
6.960008
6.971994
7.120022
7.199994

----------------------------------------------------------------
Tc0=6.771994
for Tca in 6.771994 6.771994 6.819980 6.808001 6.960008 6.971994 7.120022 7.199994
do
	echo "($Tca - $Tc0)/$Tc0" | kbc -np
done

0.000000
0.000000
0.007086
0.005317
0.027763
0.029533
0.051392
0.063201

----------------------------------------------------------------
Now divide by a

(N/A)
0.000000
0.708600
0.483364
0.555260
0.579078
0.513920
0.316005