function best_control=simple_control(eta,beta)
best_control=eta*beta/(1+beta^2);
end