%% G ~ g^2/k^3 rule

K_N =  [1.60 1.70 1.80 1.90 2.00 2.10 2.20 2.30 2.40 2.50];
%DG
g_DG   = [0.08 0.20 0.31 0.42 0.55 0.61 0.70 0.79 0.87 0.96];
k_DG   = [0.21 0.31 0.37 0.42 0.47 0.47 0.47 0.52 0.52 0.47];
G_DG   = g_DG.^2./k_DG.^3./K_N;
%SG
g_SG   = [0.08 0.19 0.30 0.40 0.49 0.59 0.68 0.76 0.85 0.93];
k_SG   = [0.26 0.37 0.42 0.42 0.47 0.47 0.52 0.52 0.52 0.47];
G_SG   = g_SG.^2./k_SG.^3./K_N;
%LR
g_LR   = [0.12 0.18 0.25 0.32 0.39 0.47 0.55 0.64 0.72 0.81];
k_LR   = [0.31 0.31 0.37 0.37 0.37 0.37 0.37 0.37 0.37 0.37];
G_LR   = g_LR.^2./k_LR.^3./K_N;
%LD
g_LD   = [0.05 0.11 0.19 0.27 0.35 0.44 0.54 0.63 0.73 0.83];
k_LD   = [0.26 0.26 0.31 0.31 0.31 0.31 0.31 0.31 0.31 0.31];
G_LD   = g_LD.^2./k_LD.^3./K_N;
% nu=0
g_CL   = [0.22 0.27 0.32 0.37 0.44 0.51 0.59 0.67 0.76 0.85];
k_CL   = [0.73 0.73 0.73 0.73 0.63 0.52 0.52 0.52 0.52 0.42];
G_CL   = g_CL.^2./k_CL.^3./K_N;
%% Plot

figure;
% plot(K_N,G_SG); hold on;
% plot(K_N,G_DG); hold on;
% plot(K_N,G_LR); hold on;
% plot(K_N,G_LD); hold on;
% plot(K_N,G_CL,'--k'); hold on;
k = 0.3;
plot(K_N,g_SG.^2/k.^3./K_N); hold on;
plot(K_N,g_DG.^2/k.^3./K_N); hold on;
plot(K_N,g_LR.^2/k.^3./K_N); hold on;
plot(K_N,g_LD.^2/k.^3./K_N); hold on;
plot(K_N,g_CL.^2/k.^3,'--k'); hold on;
% plot(K_N,(G_DG+G_SG+G_LR+G_LD)/4);