SC15

clear
D=input("Hole diameter in millimeters : ")
SG=input("Density of rock in t/m^3 : ")
E=input("Explosive density in kg/m^3 : ")
N=input("Number of holes in a blast : ")
//BH=Bench Height
BH=D/15
//B=Berden
B=(30*D)/1000
//S=Spacing
S=1.15*B
//J=Subdrill
J=(5*D)/1000
//C=Charge length
C=(20*D)/1000
//T=Stemming
T=B
//SMS=Stemming Material Size
SMS=D/10
//L=Hole length
L=BH+J
//V=Blast volume
V=B*S*BH*N
//Bt=Blasted tonnes
Bt=V*SG
//Vb=Volume of blasthole
Vb=%pi*((D^2)/4000)*L
//MH=Mass of explosive per hole
MH=Vb*E
//PF=Powder factor
PF=(MH*N)/(Bt*1000)
printf("Bench height in meters : %.3f \n",BH)
printf("Burden in meter : %.3f \n",B)
printf("Spacing in meters : %.3f \n",S)
printf("Subdrill in meters : %.3f \n",J)
printf("Charge length in meters: %.3f \n",C)
printf("Stemming in meters : %.3f \n",T)
printf("Stemming material size : %.3f \n",SMS)
printf("Hole length in meters : %.3f \n",L)
printf("Blast volume in m^3 : %.3f \n",V)
printf("Blasted tonnes in tonnes : %.3f \n",Bt)
printf("Volume of blasthole : %.3f \n",Vb)
printf("Mass of explosive per hole (kg) : %.3f \n",MH)
printf("Powder factor : %.3f \n",PF)
//BSR=Burden Stiffness Ratio
BSR=BH/B
if(BSR>=2)
    printf('Good Fragmentation\n')
elseif(BSR>3.5)
    printf('Very Good Fragmentation\n')
else
    printf('Not Good\n')
end