//chapter9,Example9_6,pg 239

a=3.615*10^-10

d111=a/sqrt(1+1+1)//for 111 plane 

theta=21.7*(%pi/180)//converting into radian

lam=2*d111*sin(theta)

printf("wavelength of X-rays\n")

disp(lam)