function [theta] = theta(U);
global a b r
theta = min(max(a-b*U^r,0),1);