Get_A.m

function A = Get_A(Ope_h, Theta)
    h1 = Ope_h(1); h2 = Ope_h(2); h3 = Ope_h(3);
    Theta_1 = Theta(1); Theta_2 = Theta(2); Theta_3 = Theta(3); 
    Theta_4 = Theta(4); Theta_5 = Theta(5); Theta_6 = Theta(6);
    
    A = [-(3*218^(1/2)*Theta_1)/(308*h1^(1/2)) - (69*218^(1/2)*sign(h1 - h3)^2*(Theta_4 + 1))/(30800*abs(h1 - h3)^(1/2))-...
        (69*218^(1/2)*abs(h1 - h3)^(1/2)*dirac(h1 - h3)*(Theta_4 + 1))/7700, 0, ...
        (69*218^(1/2)*sign(h1 - h3)^2*(Theta_4 + 1))/(30800*abs(h1 - h3)^(1/2)) + (69*218^(1/2)*abs(h1 - h3)^(1/2)*...
        dirac(h1 - h3)*(Theta_4 + 1))/7700;
        0, - (9*218^(1/2)*(Theta_5 + 1))/(3080*h2^(1/2)) - (3*218^(1/2)*Theta_2)/(308*h2^(1/2)) - (27*218^(1/2)*sign(h2 - h3)^2*...
        (Theta_6 + 1))/(12320*abs(h2 - h3)^(1/2)) - (27*218^(1/2)*abs(h2 - h3)^(1/2)*dirac(h2 - h3)*(Theta_6 + 1))/3080,...
        (27*218^(1/2)*sign(h2 - h3)^2*(Theta_6 + 1))/(12320*abs(h2 - h3)^(1/2)) + (27*218^(1/2)*abs(h2 - h3)^(1/2)*dirac(h2 - h3)*...
        (Theta_6 + 1))/3080;
        (69*218^(1/2)*sign(h1 - h3)^2*(Theta_4 + 1))/(30800*abs(h1 - h3)^(1/2)) + (69*218^(1/2)*abs(h1 - h3)^(1/2)*dirac(h1 - h3)*...
        (Theta_4 + 1))/7700, (27*218^(1/2)*sign(h2 - h3)^2*(Theta_6 + 1))/(12320*abs(h2 - h3)^(1/2)) + (27*218^(1/2)*abs(h2 - h3)^...
        (1/2)*dirac(h2 - h3)*(Theta_6 + 1))/3080, - (3*218^(1/2)*Theta_3)/(308*h3^(1/2)) - (69*218^(1/2)*sign(h1 - h3)^2*(Theta_4 +...
        1))/(30800*abs(h1 - h3)^(1/2)) - (27*218^(1/2)*sign(h2 - h3)^2*(Theta_6 + 1))/(12320*abs(h2 - h3)^(1/2)) - (69*218^(1/2)*...
        abs(h1 - h3)^(1/2)*dirac(h1 - h3)*(Theta_4 + 1))/7700 - (27*218^(1/2)*abs(h2 - h3)^(1/2)*dirac(h2 - h3)*(Theta_6 + 1))/3080];
    A = real(A);
end