function out1 = fcn136 (beta)

cldmg  =   1 ;  % Zahlenwerte für die
l0dl   =   1 ;  % beiden Problemparameter

adl    =   cos (beta) ;
bdl    =   sin (beta) ;
lsdl   =   sqrt ((1-adl)^2 + (l0dl+bdl)^2) ;
sdelta =   (1-adl)/lsdl ;
cdelta =   (l0dl+bdl)/lsdl ;
fcdmg  =   cldmg * (lsdl-l0dl) ;
f      =   adl/2 - fcdmg * (adl*cdelta + bdl*sdelta) ;

out1 = f ;
