blob: a3dc6d537adbe77b71dbeac846af6c2f2f244071 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
//Example 5_13_u1
clc();
clear;
//To determine the crystal structure and indices of plane and lattice parameter of the material
theta1=20.7
peak1=sin((theta1/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 1 is %d\n",peak1)
theta2=28.72
peak2=sin((theta2/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 2 is %.1f\n",peak2)
theta3=35.36
peak3=sin((theta3/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 3 is %.1f\n",peak3)
theta4=41.07
peak4=sin((theta4/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 4 is %.1f\n",peak4)
theta5=46.19
peak5=sin((theta5/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 5 is %.1f\n",peak5)
theta6=50.9
peak6=sin((theta6/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 6 is %.1f\n",peak6)
theta7=55.28
peak7=sin((theta7/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 7 is %.1f\n",peak7)
theta8=59.42
peak8=sin((theta8/2)*%pi/180)^2/0.0308
printf("Indices of the plane corresponding to peak 8 is %.1f\n",peak8)
printf("The material corresonds to bcc structure\n")
lamda=0.07107 //units in nm
theta=29.71 //units in degrees
d400=lamda/(2*sin(theta*(%pi/180))) //units in nm
hkl=16
a=d400*sqrt(hkl) //units in nm
printf("Lattice parameter of the material a=%.4fnm",a)
|
