% 定义粒子状态，能级对应的未知数
n_states = {
   
    'n1', '1s5';
    'n2', '1s4';
    'n3', '1s3';
    'n4', '1s2';
    'n5', '2p10';
    'n6', '2p9';
    'n7', '2p8';
    'n8', '2p7';
    'n9', '2p6';
    'n10', '2p5';
    'n11', '2p4';
    'n12', '2p3';
    'n13', '2p2';
    'n14', '2p1';
    'n15', '2s3d';
    'n16', '3p';
    'n17', 'hl';
    'n18', 'Ar2*';
    'n19', 'Ar2+';
    'n20', 'Ar+'
};

% 定义 Qegn 数组
Qegn = zeros(1, 17); % 存储 Qegn 值的数组

% 计算 Qegn 值
Qegn(1) = 2.7e-9 * exp(-11.9 / te); % Qegn(1)
Qegn(2) = 3.5e-9 * exp(-12.3 / te); % Qegn(2)
Qegn(3) = 5.5e-10 * exp(-12.2 / te); % Qegn(3)
Qegn(4) = 8.3e-9 * exp(-12.7 / te); % Qegn(4)
Qegn(5) = 2.0e-9 * exp(-13.0 / te); % Qegn(5)
Qegn(6) = 1.9e-9 * exp(-13.5 / te); % Qegn(6)
Qegn(7) = 2.2e-9 * exp(-13.6 / te); % Qegn(7)
Qegn(8) = 1.1e-9 * exp(-13.9 / te); % Qegn(8)
Qegn(9) = 1.5e-9 * exp(-13.5 / te); % Qegn(9)
Qegn(10) = 1.0e-9 * exp(-13.6 / te); % Qegn(10)
Qegn(11) = 1.0e-9 * exp(-14.0 / te); % Qegn(11)
Qegn(12) = 1.3e-9 * exp(-13.7 / te); % Qegn(12)
Qegn(13) = 7e-10 * exp(-13.8 / te);  % Qegn(13)
Qegn(14) = 2.2e-9 * exp(-13.9 / te); % Qegn(14)
Qegn(15) = 8.0e-9 * exp(-14.5 / te); % Qegn(15)
Qegn(16) = 4.0e-9 * exp(-14.9 / te); % Qegn(16)
Qegn(17) = 4.0e-9 * exp(-15.8 / te); % Qegn(17)

% 定义 Qse 数组
Qse = zeros(20, 20); % 初始化 Qse 数组，20行20列

% 计算 Qse 值
Qse(1, 5) = 3e-7 * exp(-2 / te);  % Qse(1, 5)
Qse(2, 5) = 1.5e-7 * exp(-2 / te); % Qse(2, 5)
Qse(3, 5) = 1.5e-7 * exp(-2 / te); % Qse(3, 5)
Qse(1, 6) = 6e-7 * exp(-2 / te);  % Qse(1, 6)
Qse(1, 7) = 1e-7 * exp(-2 / te);  % Qse(1, 7)
Qse(2, 7) = 4e-7 * exp(-2 / te);  % Qse(2, 7)
Qse(4, 7) = 3e-8 * exp(-2 / te);  % Qse(4, 7)
Qse(1, 8) = 2e-8 * exp(-2 / te);  % Qse(1, 8)
Qse(2, 8) = 2e-7 * exp(-2 / te);  % Qse(2, 8)
Qse(3, 8) = 8e-8 * exp(-2 / te);  % Qse(3, 8)
Qse(4, 8) = 1e-8 * exp(-2 / te);  % Qse(4, 8)
Qse(1, 9) = 2e-7 * exp(-2 / te);  % Qse(1, 9)
Qse(2, 9) = 7e-8 * exp(-2 / te);  % Qse(2, 9)
Qse(4, 9) = 1e-7 * exp(-2 / te);  % Qse(4, 9)
Qse(2, 10) = 7e-8 * exp(-2 / te); % Qse(2, 10)
Qse(1, 11) = 1e-8 * exp(-2 / te); % Qse(1, 11)
Qse(3, 11) = 5e-7 * exp(-2 / te); % Qse(3, 11)
Qse(1, 12) = 2e-8 * exp(-2 / te); % Qse(1, 12)
Qse(2, 12) = 7e-8 * exp(-2 / te); % Qse(2, 12)
Qse(4, 12) = 3e-7 * exp(-2 / te); % Qse(4, 12)
Qse(1, 13) = 1e-8 * exp(-2 / te); % Qse(1, 13)
Qse(2, 13) = 6e-8 * exp(-2 / te); % Qse(2, 13)
Qse(3, 13) = 2e-7 * exp(-2 / te); % Qse(3, 13)
Qse(4, 13) = 6e-8 * exp(-2 / te); % Qse(4, 13)
Qse(4, 14) = 6e-8 * exp(-2 / te); % Qse(4, 14)
Qse(1, 16) = 5e-8 * exp(-3 / te); % Qse(1, 16)
Qse(2, 16) = 5e-8 * exp(-3 / te); % Qse(2, 16)
Qse(3, 16) = 5e-8 * exp(-3 / te); % Qse(3, 16)
Qse(4, 16) = 5e-8 * exp(-3 / te); % Qse(4, 16)
Qse(3, 16) = 5e-8 * exp(-3 / te); % Qse(3, 16)

% 定义温度 tg
tg = 300; % 可以根据需要调整

% 定义 Qtrn 数组
Qtrn = zeros(20, 20); % 初始化 Qtrn 数组，20行20列

% 计算 Qtrn 值
Qtrn(1, 2) = 1e-7 * te^(-0.6);         % Qtrn(1, 2)
Qtrn(3, 4) = 1e-7 * te^(-0.6);         % Qtrn(3, 4)
Qtrn(3, 2) = 2e-8 * te^(-0.6);         % Qtrn(3, 2)
Qtrn(1, 4) = 2e-8 * te^(-0.6);         % Qtrn(1, 4)

% 填充 5 到 14 行和列的值
for i = 5:14
    for j = 5:14
        if i ~= j
            Qtrn(i, j) = 4e-7 * te^(-0.6);   % Qtrn(5-14, 5-14)
        end
    end
end

% 填充 5 到 14 行和第 15 列的值
for i = 5:14
    Qtrn(i, 15) = 2e-6 * exp(-2 / te);  % Qtrn(5-14, 15)
end

Qtrn(15, 16) = 2e-5;                    % Qtrn(15, 16)
Qtrn(16, 15) = 2e-5;                    % 对称关系
Qtrn(17, 15) = 2e-5;                    % Qtrn(17, 15)
Qtrn(17, 16) = 2e-5;                    % Qtrn(17, 16)
Qtrn(15, 17) = 2e-5;                    % 对称关系
Qtrn(16, 17) = 2e-5;                    % 对称关系

% 填充 1 到 4 行和 14 行的值
for i = 1:4
    Qtrn(14, i) = 3e-11 * (tg / 300)^0.5;  % Qtrn(14, 1-4)
    Qtrn(13, i) = 1e-11 * (tg / 300)^0.5;  % Qtrn(13, 1-4)
    Qtrn(12, i) = 3e-11 * (tg / 300)^0.5;  % Qtrn(12, 1-4)
    Qtrn(11, i) = 3e-11 * (tg / 300)^0.5;  % Qtrn(11, 1-4)
    Qtrn(8,  i) = 4e-11 * (tg / 300)^0.5;  % Qtrn(8, 1-4)
    Qtrn(7,  i) = 4e-11 * (tg / 300)^0.5;  % Qtrn(7, 1-4)
    Qtrn(6,  i) = 3e-11 * (tg / 300)^0.5;  % Qtrn(6, 1-4)
    Qtrn(5,  i) = 1.5e-11 * (tg / 300)^0.5; % Qtrn(5, 1-4)
end

% 填充特定的 Qtrn 值
Qtrn(12, 11) = 2e-11 * (tg / 300)^0.5;  % Qtrn(12, 11)
Qtrn(12, 9)  = 2e-11 * (tg / 300)^0.5;  % Qtrn(12, 9)
Qtrn(10, 9)  = 2.5e-11 * (tg / 300)^0.5; % Qtrn(10, 9)
Qtrn(10, 7)  = 1.5e-11 * (tg / 300)^0.5; % Qtrn(10, 7)
Qtrn(9, 8)   = 2e-11 * (tg / 300)^0.5;   % Qtrn(9, 8)
Qtrn(9, 7)   = 1e-11 * (tg / 300)^0.5;   % Qtrn(9, 7)
Qtrn(8, 7)   = 1e-11 * (tg / 300)^0.5;   % Qtrn(8, 7)
Qtrn(8, 6)   = 2e-11 * (tg / 300)^0.5;   % Qtrn(8, 6)
Qtrn(7, 6)   = 2.5e-11 * (tg / 300)^0.5; % Qtrn(7, 6)
Qtrn(6, 5)   = 3e-11 * (tg / 300)^0.5;   % Qtrn(6, 5)

% 填充 16 行和 5 到 14 列的值
for i = 5:14
    Qtrn(16, i) = 1e-10 * (tg / 300)^0.5; % Qtrn(16, 5-14)
    Qtrn(i, 16) = 1e-10 * (tg / 300)^0.5; % 对称关系
end

Qtrn(16, 15) = 1e-10 * (tg / 300)^0.5; % Qtrn(16, 15)
Qtrn(15, 16) = 1e-10 * (tg / 300)^0.5; % 对称关系

% 填充 17 行和 15 到 16 列的值
Qtrn(17, 15) = 1e-10 * (tg / 300)^0.5; % Qtrn(17, 15)
Qtrn(17, 16) = 1e-10 * (tg / 300)^0.5; % Qtrn(17, 16)
Qtrn(15, 17) = 1e-10 * (tg / 300)^0.5; % 对称关系
Qtrn(16, 17) = 1e-10 * (tg / 300)^0.5; % 对称关系

% 定义 Qtr 数组（与 Qtrn 相同）
Qtr = Qtrn;

% 定义 Qote 和 Qota 数组
Qote = zeros(20, 20); % 初始化 Qote 数组
Qota = zeros(20, 20); % 初始化 Qota 数组

% 计算 Qote 值
Qote(1:4, 20) = 2e-7 * exp(-6.2 / te);  % Qote(1-4, 20)
Qote(5:14, 20) = 2e-6 * exp(-4.4 / te); % Qote(5-14, 20)
Qote(15, 20) = 6e-6 * exp(-2.4 / te);   % Qote(15, 20)
Qote(16, 20) = 2e-5 * exp(-2.2 / te);   % Qote(16, 20)
Qote(17, 20) = 2e-4 * exp(-0.5 / te);   % Qote(17, 20)

Qote(20, 17) = 2e-27 * te^(-4.5);        % Qote(20, 17)

% 计算 Qota 值
Qota(1:4, 20) = 5e-10 * (tg / 300)^0.5;   % Qota(1-4, 20)
Qota(5:14, 20) = 5e-10 * (tg / 300)^0.5;  % Qota(5-14, 20)
Qota(15, 20) = 7e-10 * (tg / 300)^0.5;    % Qota(15, 20)
Qota(16, 20) = 7e-10 * (tg / 300)^0.5;    % Qota(16, 20)
Qota(17, 20) = 7e-10 * (tg / 300)^0.5;    % Qota(17, 20)

Qota(20, 17) = 1.5e-28 * (tg / 300)^-2.5;  % Qota(20, 17)
Qota(18, 1:4) = 1e-32;                      % Qota(18, 1-4)

% 计算 Qote(18, 1-4) 和 Qota(18, 1-4)
Qote(18, 1:4) = 1e-8 * exp(-1 / te);        % Qote(18, 1-4)
Qote(18, 19) = 2e-8 * exp(-6.2 / te);       % Qote(18, 19)
Qota(18, 19) = 7e-10 * (tg / 300)^0.5;      % Qota(18, 19)
Qota(20, 19) = 2.5e-31 * (tg / 300)^(-1);   % Qota(20, 19)

% 计算 Qote(19, ...) 和 Qota(19, ...)
Qote(19, 20) = 1e-6 * exp(-2 / te);                             % Qote(19, 20)
Qote(19, 19) = 1e-7 * te^(-0.6) * (tg / 300)^(-0.6);          % Qote(19, 19)
Qote(19, 1:4) = 1e-8 * te^(-0.6) * (tg / 300)^(-0.6);         % Qote(19, 1-4)
Qote(19, 5:14) = 1e-8 * te^(-0.6) * (tg / 300)^(-0.6);        % Qote(19, 5-14)
Qote(19, 15) = 1e-8 * te^(-0.6) * (tg / 300)^(-0.6);          % Qote(19, 15)
Qote(19, 16) = 1e-8 * te^(-0.6) * (tg / 300)^(-0.6);          % Qote(19, 16)
Qote(19, 17) = 6e-8 * te^(-0.6) * (tg / 300)^(-0.6);          % Qote(19, 17)

% 定义 A 数组
A = zeros(20, 20); % 初始化 A 数组，20行20列

% 填充 A 值
A(14, 4) = 4.5e7;   % A(14, 4)
A(13, 4) = 1.5e7;   % A(13, 4)
A(13, 3) = 1.2e7;   % A(13, 3)
A(13, 2) = 1.8e6;   % A(13, 2)
A(13, 1) = 6.4e6;   % A(13, 1)
A(12, 4) = 2.2e7;   % A(12, 4)
A(12, 2) = 8.7e6;   % A(12, 2)
A(12, 1) = 3.8e6;   % A(12, 1)
A(11, 4) = 1.7e7;   % A(11, 4)
A(11, 3) = 1.9e7;   % A(11, 3)
A(10, 2) = 4.0e7;   % A(10, 2)
A(9, 4)  = 5.0e6;    % A(9, 4)
A(9, 2)  = 4.9e6;    % A(9, 2)
A(9, 1)  = 2.5e7;    % A(9, 1)
A(8, 4)  = 1.0e6;    % A(8, 4)
A(8, 3)  = 2.4e6;    % A(8, 3)
A(8, 2)  = 2.5e7;    % A(8, 2)
A(8, 1)  = 5.0e6;    % A(8, 1)
A(7, 4)  = 1.5e6;    % A(7, 4)
A(7, 2)  = 2.2e7;    % A(7, 2)
A(7, 1)  = 9.3e6;    % A(7, 1)
A(6, 1)  = 3.3e7;    % A(6, 1)
A(5, 3)  = 1.0e6;    % A(5, 3)
A(5, 2)  = 5.4e6;    % A(5, 2)
A(5, 1)  = 1.9e7;    % A(5, 1)
A(16, 1) = 2e6;      % A(16, 1) - 根据需要添加
A(16, 2) = 2e6;      % A(16, 2)
A(16, 3) = 2e6;      % A(16, 3)
A(16, 4) = 2e6