#Generierung fuer TA23 (6) mit gleichm"a"siger Fahrbahnausleuchtung
#Wiederaufnahme 11.08.2008
#Wiederaufnahme: 08.09.2008
BEGIN{asta2=-8;aend2=8;asta=-4;aend=4;bsta=-3;bend=5;dphih=0.1862106;EHV=20;angtol=dphih/5;finest=0.1;emin=0.001
for (b=bsta-dphih;b<=bend+dphih;b+=finest*dphih)
 {for (a=asta2-dphih;a<=aend2+dphih;a+=finest*dphih)
  {{esoll=emin}
    {if (a>=asta-angtol && a<=aend+angtol) {esoll=Eopt(EHV, 10, 0.7, b)} else {esoll=emin}} 		#L1--R1 und HV	
    {if (b<=5+dphih+angtol && b>=5-dphih-angtol) {if (a>=-dphih-angtol && a<=dphih+angtol) {esoll=2.5} else {esoll=2.0}}} #L5--R5	
    {printf("%8.3f   %8.3f   %8.3f\n",b,a,esoll)}
  }
  print ""
 }	
}

function Eopt(Eziel,eziel,h,phii) {
 pi=3.141592654;phir=phii*pi/180
 Ephi=Eziel/(eziel/h*cos(0-phir-pi/2-atan2(h/eziel,1)))^2
 return Ephi 
 }


# 10/(25/0.4*cos((0-x+(90+atan(0.4/25)*180/pi))*pi/180))**2