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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
//Example 30.4
n=2;//Shell
number_subshell=0;//To store value of number of subshells
for i=0:1:n-1
number_subshell=number_subshell+1;
end
//To calculate number of rows required for the table
rows=0;
for l=0:1:number_subshell-1
rows=rows+2*(2*l+1)
end
//Defining string matrices
A=repmat(' ',[rows 4]);//To list n, l, m_l and m_s
B=repmat(' ',[rows 3]);//To list subshell, number of electrons in each subshell and total number of electrons
C=['n' 'l' 'm_l' 'm_s' 'Subshell' 'Total in subshell' 'Total in shell'];//Header row
//Function to return subshell character
function [s]=subshell(l)
if l==0
s='s';
elseif l==1
s='p';
elseif l==2
s='d';
elseif l==3
s='f';
elseif l==4
s='g';
elseif l==5
s='h';
elseif l==6
s='i';
end
endfunction
//Function defined for upto l=6
r=1;//Row number
total=0;//To store total number of electrons in the shell
for l=0:1:number_subshell-1
m_l=l;
B(r,1)=string(n)+subshell(l);//Subshell
B(r,2)=string(2*(2*l+1));//Number of electrons in subshell
total=total+(2*(2*l+1));
for i=1:1:(2*(2*l+1))
A(r,1)=string(n);//n
A(r,2)=string(l);//l
A(r,3)=string(m_l)//m_l
if modulo(i,2)==0
m_l=m_l-1;
end//m_l (decreases by 1 after every two rows)
if modulo(r,2)==0
A(r,4)='-1/2';
else
A(r,4)='+1/2';
end//m_s (positive and negative sign alternates with every row)
r=r+1;
end
end
B(1,3)=string(total);//Total number of electrons in shell
disp([C;A B])//Combining the matrices
//Openstax - College Physics
//Download for free at http://cnx.org/content/col11406/latest
|
