Yes, MTEX can do so. The way is as follows:
Assume we have the Santa Fe sample ODF with cubic crystal symmetry
odf = SantaFe
cs = symmetry('cubic')
First define the Euler angles of the fibre
phi1 = linspace(0*degree,90*degree);
Phi = 45*degree;
phi2 = 0*degree;
Next the corresponding orientations
ori =orientation('euler',ph1,Phi,phi2,cs)
In order to evaluate the ODF we use the function
eval and plot the values against phi1
plot(phi1./degree,eval(odf,ori))
adjust axis limimts
xlim([0,90])
ylim([0,1.4])
The plot should look like this

Hope this helps,
Ralf.