odf = calcODF(pf);
figure(1)
plotpdf(odf, [Miller(1,1,1), Miller(2,0,0),Miller(2,2,0)],'complete','contour',0:1:15, 'antipodal');colorbar
figure(2)
plot(odf, 'phi2',[0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90]*degree,'contour',0:1:50,'antipodal');colorbar
figure(3)
plot(odf,'surf3','contour',0:1:50);colorbar
caxis(gca,[0,50])