alpha = 0.1;
Cl_deltaA =DE2019.Cl_deltaA_0 + DE2019.Cl_deltaA_alpha * alpha + DE2019.Cl_deltaA_alpha2 * alpha^2;
Cl_deltaR =DE2019.Cl_deltaR_0 + DE2019.Cl_deltaR_alpha * alpha + DE2019.Cl_deltaR_alpha2 * alpha^2;
Cm_deltaA =DE2019.Cm_deltaA_0 + DE2019.Cm_deltaA_alpha * alpha + DE2019.Cm_deltaA_alpha2 * alpha^2;
Cm_deltaE = DE2019.Cm_deltaE_0 + DE2019.Cm_deltaE_alpha * alpha + DE2019.Cm_deltaE_alpha2 * alpha^2;
Cn_deltaA =DE2019.Cn_deltaA_0 + DE2019.Cn_deltaA_alpha * alpha + DE2019.Cn_deltaA_alpha2 * alpha^2;
Cn_deltaR =DE2019.Cn_deltaR_0 + DE2019.Cn_deltaR_alpha * alpha + DE2019.Cn_deltaR_alpha2 * alpha^2;
Ba = [[Cl_deltaA,0,Cl_deltaR];[Cm_deltaA,Cm_deltaE,0];[Cn_deltaA,0,Cn_deltaR]]
Ba =
0.1474 0 0.0090
0.0142 -1.8535 0
0 0 0.0759