SC15

clear
D=input("Hole Diameter(mm)=")
SG=input("Density of rock(t/m^3)=")
N=input("Number of holes=")
ED=input("Explosive density(kg/m^3)=")
BH=D/10
B=30*D/1000
S=1.15*30*D/1000
J=10*D/1000
LC=20*D/1000
T=30*D/1000
BSR=((D/15)/(30*D))*1000
SMS=D/10
L=(BH)+(J)
V=(B)*(S)*(BH)*(N)
BT=(V)*(SG)
Vb=(%pi)*(D^2)*(L)/4000
MP=(Vb)*ED
PF=(MP)*(N)/((V)*1000)
printf("Bench Height(m)=%.3f \nBurden(m)=%.3f \nSpacing(m)=%.3f \nSubdrill(m)=%.3f \nCharge Length(m)=%.3f \nStemming(m)=%.3f \nBSR=%.3f \nSMS=%.3f \nHole length(m)=%.3f \nBlast volume(m^3)=%.3f \nBlast tonnes(T)=%.3f \nVolume of blast hole(m^3)=%.3f \nMass of explosive per hole(kg)=%.3f \nPF(kg/m^3)=%.3f",BH,B,S,J,LC,T,BSR,SMS,L,V,BT,Vb,MP,PF)