% make orbit
function make_orbit

%     i=(1:22)';
%     v=0.2*(1-2*rand([1 22])');
%     setpv('HCOR',v ,i);
    setpv('HCOR',0.5 ,7
    );

end