beta0 = 0
beta_max = beta0+max(gamma) %radiani

for i=1:21
    beta(i)=beta0+gamma(i);
end
beta

for i=1:21
    s(i)=(cos(beta(i))-cos(beta_max)+(sin(beta_max)-sin(beta(i)))*cot(delta_rad))*l2;
end
s
st = e1*(tan(beta_max)-tan(min(gamma)))