>> radius= 10
>> z=3
>> alphaa=0
>> lamda=7
>> r= 0.1*radius
>> tr=(1/z)*(16*pi)*(r)*(sind((1/3)*(atand(radius/(lamda*r)))))^2
>> trr=tr/10
>> x=[0*trr, 0.2*trr, 0.4*trr, 0.6*trr, 0.8*trr, 1.0*trr]
>> yup=[0*trr, 0.1*trr, 0.11*trr, 0.09*trr, 0.05*trr, 0*trr]
>> ydown=[0*trr, 0.03*trr, 0.02*trr, 0.013*trr, 0.008*trr, 0*trr]
>> h1=line(x,yup)
>> hold
>> h2=line(x,ydown)
>> dir=[0 0 1]
>> alphaconnected=((2/3)*(atand(radius/lamda*r)))-alphaa
>> rotate(h1,dir,alphaconnected)
>> rotate(h2,dir,alphaconnected)
Posts
