--> //initializations

--> Vm=1000 //ml (volume of mould)
 Vm  = 

   1000.


--> Mc=2.1 //Kg (mass of compacted soil)
 Mc  = 

   2.1


--> Ml=1.635 //Kg (mass in loose state)
 Ml  = 

   1.635


--> Dd=1.5 //Mg/m^3 (in-situ density)
 Dd  = 

   1.5


--> G=2.7 //(specific gravity)
 G  = 

   2.7


--> Dw=1 //g/m^3 (density of water)
 Dw  = 

   1.


--> //calculations

--> Dmax=Mc*1000/Vm // (max density)
 Dmax  = 

   2.1


--> Dmin=Ml*1000/Vm // (min density)
 Dmin  = 

   1.635


--> e=((G*Dw)/Dd)-1 //(void ratio)
 e  = 

   0.8


--> emin=(G*Dw-Dmax)/(Dmax-Dw) //(min density)
 emin  = 

   0.5454545


--> emax=(G*Dw-Dmin)/(Dmin-Dw) //(max density)
 emax  = 

   1.6771654


--> Dr=((emax-e)/(emax-emin))*100 //(relative density)
 Dr  = 

   77.507906


--> //results

--> printf('relative density=%f',Dr) //%
relative density=77.507906