cs = crystalSymmetry('222')

odf = uniformODF(cs)

%%

ori = orientation.rand(cs);

volume(odf,ori,30*degree)*6

f = fibre.rand(cs);

volume(odf,f,30*degree)*6

%%

mtexdata dubna

plot(pf)
%%
odf = calcODF(pf)


%%

pdf = calcPDF(odf,odf.CS.cAxis)
pdf.eval(zvector)



%%

ipdf = calcPDF(odf,[],zvector)

ipdf.eval(odf.CS.cAxis)

%%

plot(ipdf)

%%

[value,pos] = max(pdf,'numLocal',2)
%%

annotate(pos)

%%
plotPDF(odf,odf.CS.cAxis)
