psi = calcKernel(grains.meanOrientation)
% compute the ODF with the kernel psi
odf = calcODF(ebsd('Material').orientations,'kernel',psi)
%
%to get polefig with colorbar
%plot odf
figure
h = [Miller(1,0,0,odf.CS),Miller(1,1,0,odf.CS),Miller(1,1,1,odf.CS)];
%plotPDF(odf,h,'antipodal','silent')
plotPDF(odf,h,'colorrange',[0 5],'antipodal','silent')
CLim(gcm,'equal');