Hi Lawrence Ko,
my first advice is to go through the documentation a bit. Here are the basic steps:
1: define crystal symmetry
cs = symmetry('m-3m','mineral','Al')
2: import your orientations given you have your Euler angles in the variables
ph1, Phi, phi2ori = orientation('Euler',phi1,Phi,phi2,cs)
3: compute an ODF from the orientations
I hope this helps. Note all commands have a bunch of options. These are documented in the help.
Ralf.