function [rx, p]=ar2(a1,a2,px);
rx=[px -a1*px./(1+a2)]';
p=[-a1*px./(1+a2) (-a2+a1.^2./(1+a2))*px]';