function [loss]=lee(R,f,h_1,h_2,ty);
% Loss in decibels from the Lee model

% r_km
% f_mhz
% h_1 (m)
% h_2 (m)
% Type: see below


if (ty==1) % "open terrain"
    l0=35.75;
    gam=4.35;
    n=2;
elseif (ty==2) % "suburban"
    l0=49.5;
    gam=3.84;
    n=2;
elseif (ty==3) % philadelphia
    l0=28.59;
    gam=3.68;
    n=3;
elseif (ty==4) % newark
    l0=21.29;
    gam=4.31;
    n=3;
elseif (ty==5) % Tokyo
    l0=43.89;
    gam=3.05;
    n=3;
end
loss=l0+10*n*log10(f)+10*gam*log10(R)-20*log10(h_1)-10*log10(h_2);