net = cho_load('beamgap3.net');
dq = cho_dc(net);
cho_display(net,dq);

dy = cho_dq_view(dq,net,'b','y')
dtheta = cho_dq_view(dq,net,'b','rz')
