Происхождение вселенной

Страница 17

c write(*,*) fif*180./pi,xt,r1ef,r2ef

c PAUSE' '

P(1)=(Po(i)+Py(i)*alf*dsin(fif+dfi/2.)+Pza(i)*alf*r1ef

* -Pzb(i)*bett*r2ef)*q*L*dfi*R(i)*sh

c write(*,*) P/(sh*L*R(i)*dfi),xt,fif

c pause

P(2)=Po(i)*q*L*dfi*R(i)*sh

P(3)=Py(i)*alf*dsin(fif+dfi/2.)*q*L*dfi*R(i)*sh

P(4)=Pza(i)*alf*r1ef*q*L*dfi*R(i)*sh

P(5)=-Pzb(i)*bett*r2ef*q*L*dfi*R(i)*sh

do 1000 ik=1,5

Pyc(ik)=Pyc(ik)+P(ik)*dsin(fif+dfi/2.)

Pzc(ik)=Pzc(ik)+P(ik)*dcos(fif+dfi/2.)

Myx(ik)=Myx(ik)+P(ik)*dsin(fif+dfi/2.)*xt

Mzx(ik)=Mzx(ik)+P(ik)*dcos(fif+dfi/2.)*xt

Pxc(ik)=Pxc(ik)+P(ik)*dtan(gamm(i))

1000 continue

3 continue

2 continue

xt=xt+sh

R(i)=R(i)+sh*L*dtan(gamm(i))

1 continue

do 100 ik=1,5

if(dabs(Pyc(ik)).gt.1.d-10) then

Xc(ik)=Myx(ik)/Pyc(ik)

else

write(*,*)ik,'Myx=',Myx(ik)

endif

if(dabs(Pzc(ik)).gt.1.d-10) then

Xcc(ik)=Mzx(ik)/Pzc(ik)

else

write(*,*)ik,'Mzx=',Mzx(ik)

endif

write(*,*) 'Значения аэр. сил в Ц.Д.'

100 write(*,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L

12 format(1x,'Pxc=',f10.2,/,1x,'Pyc=',f10.2/,1x,'Pzc=',f10.2/,

* 1x,'Xyc=',f15.7,/,1x,'Xzc=',f15.7)

open(1,file='aer.res')

write(1,*) 'Начальное значение угла Fi'

write(1,*) fin*180./pi

write(1,*) 'Конечное значение угла Fi'

write(1,*) fi*180./pi

write(1,*) 'Начальная Xнач'

write(1,*) Ln*5.6

write(1,*) 'Конечная Xкон'

write(1,*) Lk*5.6

write(1,*) 'Значение угла BETTA'

write(1,*) bett*180./pi

write(1,*) 'Значение угла ALFA'

write(1,*) alf*180./pi

write(1,*) 'Значение q'

write(1,*) q

write(1,*) 'Значения аэр. сил в Ц.Д.'

do 102 ik=1,5

102 write(1,12) Pxc(ik),Pyc(ik),Pzc(ik),Xc(ik)*L,Xcc(ik)*L

close(1)

stop' '

end