How to calculate cell edge "user" throughput based on distance from the eNodeB?

Hello

I'm trying to calculate throughput for all users in a cell and generate a graph based on power and throughput how the eNodeB should be positioned.

Is there any Matlab formula to calculate user throughput given the pre-defined parameters for LTE bandwidth (20MHz) and 4x4 MIMO configuration?

My current code:

for i=1:N

         Dt(i)=norm(Tr1-[x(i);y(i)]); %distance between base station and users

         Pr(i)=Pt - a * log10(Dt(i)); %RSRP of the all the users in the network

         PL_m(i)=128.1+ 37.6 * log10(Dt(i)); %Calculate Path Loss for Macro cell 

         Prs(i)=Pt - (a * log10(Dt(i))-s(i)); %RSRP of the all the users in the network with shadowing

end