[1/2 0 0] 100};
Qlab  = {<span class="string">'P'</span> <span class="string">'M'</span> <span class="string">'X'</span> <span class="string">'P'</span> <span class="string">'\Gamma'</span> <span class="string">'X'</span>};

lacuoSpec = lacuo.spinwave(Qlist,<span class="string">'hermit'</span>,false);
lacuoSpec.omega = lacuoSpec.omega*Zc;

lacuoSpec = sw_neutron(lacuoSpec);
lacuoSpec = sw_egrid(lacuoSpec,<span class="string">'component'</span>,<span class="string">'Sperp'</span>);
figure
subplot(2,1,1)
sw_p