fm_awiphan_bulge.m

%IMF du modèle de Besancon 
%Awiphan, 2016
%Fonction du Bulbe

function pm = fm_awiphan_bulge(m)
pm = zeros(size(m));

global minf msup

% Thin disk
m1 = 0.01;
m2 = 0.079;
m3 = 0.7;

c = 1;

%-----------
% Bulge
%----------
alpha_bd = -0.4;
alpha_MS = -1.5;
alpha = -2.35;


i1 = find(m>=m1 & m<m2 & m>=minf);
i2 = find(m>=m2 & m<m3);
i3 = find(m>=m3 & m<msup);

if (length(i1)>=1)
    pm(i1) = c * m(i1).^(alpha_bd);
end

if (length(i2)>=1)
    pm(i2) = c * m2^(alpha_bd)/m2^(alpha_MS) * m(i2).^(alpha_MS);
end

if (length(i3)>=1)
    pm(i3) = m2^(alpha_bd)/m2^(alpha_MS) * m3^(alpha_MS)/m3^(alpha)*m(i3).^(alpha);
end