//page 70
clc;funcprot(0);//EXAMPLE 3.11
// Initialisation of Variables
E=12;......//No. of Edges in the octahedral sites of the unit cell
S=1/4;.......//so only 1/4 of each site belongs uniquelyto each unit cell
N=E*S+1;.....//No.of site belongs uniquely to each unit cell
disp(N,"No.of octahedral site belongs uniquely to each unit cell:")