tatchaa

clear D=input("Input Blast hole diameter(mm.) ="); N=input("Input Number of Holes in Blast(รู) ="); SG=input("Input Density of rock(T.) ="); ED=input("Input Explosive density(kg.) ="); BH=D/15; B=30*D/1000; S=1.15*B; SD=3*D/1000; C=20*D/1000; SL=0.7*B; BSR=(BH/B); SMS=D/10; L=BH+SD; C=L-SL; V=B*S*BH*N; T=V*SG; Vb=(%pi*D*D*L)/4000; MP=Vb*ED; PF=(MP*N)/(V*1000); printf("Bench Height=%.3f m. \nBunden=%.3f m. \nSpacing=%.3f m. \nSubdrill=%.3f m. \nCharge length=%.3f m. \nStemming=%.3f m. \nBurden stiffness ratio=%.3f \nStemming material size=%.3f \nHole Length=%.3f m. \nCharge length=%.3f m. \nBlast volume=%.3f m^3 \nBlast tonnes=%3f T. \nVolme of blasthole=%.3f m^3 \nMass of explosive per hole=%.3f kg. \nPowder factor=%.3f kg/m^3",BH,B,S,SD,C,SL,BSR,SMS,L,C,V,T,Vb,MP,PF);

-->q=[.2 .5 .9 1.3 1.6]  q  =     0.2    0.5    0.9    1.3    1.6  -->r=[.1 .4 .9 1.2 2.1]  r  =     0.1    0.4    0.9    1.2    2.1    -->plot2d(q,r, style=2) -->plot2d(q,r, style=-3) -->g={.5 .9 1.8 2.6 3.2]  g  =     0.5    0.9    1.8    2.6    3.2  -->h=[.5 .8 1.6 2.3 2.8]  h  =     0.5    0.8    1.6    2.3    2.8  -->plot2d(g,h, style=-1) -->plot2d(g,h, style=3)

-->a=10.0:5.0:50.0  a  =     10.    15.    20.    25.    30.    35.    40.    45.    50.  -->b=3.0:3.0:27.0  b  =     3.    6.    9.    12.    15.    18.    21.    24.    27.  -->plot(a,b)