// On va crée la fonction qui permet de calculer la valeur de omega0 à partir de Xi et Tr 
function Omega_0 =calcul(Xi,Tr)
    if (Xi<0.7) then
        A=3/Xi;
    end
        if (Xi>0.7) then
        A=6*Xi;
    end
        if (Xi==0.7) then
        A=3;
    end
    Omega_0=A/Tr;
endfunction
