function cS = monazite

%
cs_Mnz = crystalSymmetry('2/m',[0.967 1 0.921], [90, 103.4, 90]*degree);
N = Miller({1,0,0},{0,1,0},{1,1,0},{0,1,1},{1,1,-1},{0,2,1},{0,1,2},{1,0,-1},{1,0,1},{3,-1,-1},{2,-1,0},{0,-1,-1},{1,-2,-1},cs_Mnz);
dist = [3.6, 4.95, 6.5, 10, 8.2, 10.65, 11.1, 6.75, 7.45, 14.65, 9.9, 6.15, 12.2];
cS = crystalShape(N./dist);