function [H_new ny]=update_hyperstate(H_old, ny) H_new=zeros(size(H_old)); H_new(1,:,:)=H_old(1,:,:); next_u=zeros(1,size(ny,2)); %generate trajectories for j=1:size(H_old,2) eta=H_old(1,j,1); beta=H_old(1,j,2); best_control=-ny(1,j); for i=1:size(H_old,1)-1 s=randn; pom=sqrt(1+best_control^2); eta=abs(eta+beta*best_control+pom*s); beta=abs(pom*beta+best_control*s); H_new(i+1,j,:)=[eta beta]; if i+1