Додатки

program totogravmagne4m

use msflib

intrinsic dsin,dcos,dabs,dsqrt,dlog,datan

real(8),allocatable:: cappa(:,:,:),cappa1(:,:,:)

CHARACTER(10)dat,tim,zon

real(8):: d(3800),d1(3200),XYZZ(47,47,3800),gtole(47,47), &

lambda(47,47),lambdai(3800),ffonty(47,47),ffonty3(47,47), &

gpolepp(47,47),acappa(3800),bbb(3800),zabj(47,47), &

vecsig(3800), acappa1(3800),ccc(3800),gtoler1(521,511), &

ffonty1(47,47), sv(20,20),sv1(40,20),carra(3800),vecsigb(3800), &

ZecsigTa(47,47),ZecsigTg(47,47),vecsiga(3800),lbdai(3800), &

ZecZT(3,521,511),xx22(20),xx11(20),yy22(20),yy11(20),svo(10,10), &

gtolera(243,238),ezabi1(63,63),su1(3800),su(3800), &

zacj(47,47),ZecsigT(3,261,256),fzabi1(3800),gpole(47,47), &

vsL(180,20),svL2(9,20),svL1(20,9),v(3800),zib(47,47),zis(47,47), &

A7(47,47,3800),A2(47,47),A3(47,47),bz1(3800),bz2(3800), &

ffonty31(47,47),Fary1(3800), ffonty32(47,47),Fary2(3800), &

ffonty33(47,47),Fary3(3800),lam1bda(47,47),lam2bda(47,47), &

alam11(3200),alam22(3200),g2tole(47,47),g1tole(47,47), &

g1polepp(47,47), g2polepp(47,47),f1fonty3(47,47),f2fonty3(47,47), &

f1fonty(47,47),f2fonty(47,47),a2cappa(3200),a1cappa(3200),z1abj(47,47), &

z2abj(47,47),b1bb(3200),b2bb(3200),z12abj(47,47) ,z21abj(47,47), &

a21cappa(3200),a12cappa(3200),f1ffonty3(47,23),f2ffonty3(47,23), &

vsLa(180,20),f3fonty3(47,47),f4fonty3(47,47), &

f3fonty(47,47),f4fonty(47,47), g3polepp(47,47), g4polepp(47,47), &

zf12(47,47),zf21(47,47), E21zab(3200),E12zab(3200), &

E2zab(3200),E1zab(3200),ak21zab(3200),ak12zab(3200), &

ak2zab(3200),ak1zab(3200),s21zab(3200),s12zab(3200), &

s2zab(3200),s1zab(3200),F21zab(47,47),F12zab(47,47),gtolerg(47,47), &

F2zab(47,47),F1zab(47,47),G21zab(47,47),G12zab(47,47), &

G2zab(47,47),G1zab(47,47),Fzabi(47,47),Gzabi(47,47),Pzabi(47,47), &

Ezabi(3200),AKzabi(3200),szabi(3200),vecsib(3800)

real(8):: bbx,bby,aay,aax,zz11,xx1,xx2,yy1,yy2,xx01,yy01,ss1,ss2,ss4

real(8):: xypi,xypp,ff,ff1,ff2,oj,ajo,bjo,a,L,rrg,Zrg,tau3,tau4,tau5, &

amja,atja,ajo1,bjo1,bjo2,zz441,zz46,zz447,tau6,ajob,rr,zib2,zis2, &

sum, bbxd,bbyd,zz44a,tau,Zr3,r4,r3,Zr,Z2r2,r2,Z2r,Zr2,tau1,tau2, &

Z2rr,tau7, sum1,sum2,sum3,sum4,r1L,r2L,ALA,rgL,ajob1,ajo2,Z2r3,Z2r31, &

desa, ai,aj,ai13,r3L,r0Lz,r0LL,r1Lz2,r2Lz,r2L2,r3z,zib1,zis1,ga, &

zz11a,zz447a,zz441a, zz44ab,aay1,aax1,SRgtL,SRgt,SRgtB,SRgtR,ajob2, &

sum5,sum6,sum7,sum8,tbb2 ,tbb3,tbb4,tbb5,tbb6,tbb7,z46,bz22,bz23, &

bz24,bz25,bz26,bz27,bz28,bz29,bz30,xypp1,xypp2,xypp3,xypp4, &

taa2 ,taa3,taa4,taa5,tcc2 ,tcc3,tcc4,tcc5, &

A11, A22, B11,B22,B12,ss1a,ss2a,ss4a ,ssa,ss,ssb,ss1b,ss2b,ss4b, &

ss1ab,ss2ab,ss4ab ,ssab

a=3

b=2

c=a+b

print*, ' to print (SIgt, SRIgtL ) '

print*,a,b,c

write(*,*) 'введите a,n'

read (*,*)a,n

write(*,*) 'a=',a,' n=',n

L=1

do i=1,n

L=L+i*a

d(i)=L

end do !20

write(*,*)'L=',L

aax1 = 1.5*4; aay1=1.5*4

DLux =0.

nx=24; ny=24; nx7=47; ny7=47

zz11=0.07; aax = 0.46*3; aay=0.46*3

zz441=1.070+0.;zz46=-2.0; zz447=0.57+0.; zz44a=1.57+0.

bbxd=0.40;bbyd=0.20*0

zz11a=0.32; zz447a=0.82+0.; zz441a=1.320+0.

zz44ab=1.82+0.; z46=2.07+0.

do j2=1, 20

xx22(j2)=aax*j2-0.15*0; xx11(j2)=xx22(j2)-aax

enddo

xx22(20)=xx22(20)+1.50*4; xx11(1)=xx11(1) -1.5*4

do j1=1, 20

yy22(j1)=aay*j1+0.0-0.15*0; yy11(j1)=yy22(j1)-aay

enddo

yy11(1)=yy11(1)-1.50*4; yy22(20)=yy22(20)+1.5*4

do k=1,8

do i=1, 20

do j=1, 20

j3=i+20*(j-1)+400*(k-1)

v(j3)=aax*aay

end do; end do;end do

do k=1,8

do i=1, 20,19

do j=1, 20

j3=i+20*(j-1)+400*(k-1)

v(j3)=aax1*aay

end do; end do;end do

do k=1,8

do i=1, 20

do j=1, 20,19

j3=i+20*(j-1)+400*(k-1)

v(j3)=aax*aay1

end do; end do;end do

do k=1,8

do i=1, 20,19

do j=1, 20,19

j3=i+20*(j-1)+400*(k-1)

v(j3)=aax1*aay1

end do; end do;end do

do i=1, nx7

bbx=0.2*3; xx01=bbx*(i-1)

do j=1, ny7

bby=0.2*3; yy01=bby*(j-1)

do j2=1, 20

! xx2=aax*j2+0.001-0.15*3; xx1=xx2-aax !50

xx2=xx22(j2)+4.75*0; xx1=xx11(j2)+4.75*0

do j1=1, 20

yy2=yy22(j1)+4.75*0;yy1=yy11(j1)+4.75*0

! yy2=aay*j1+0.001-0.15*3; yy1=yy2-aay

ss=gg1(xx1,xx2,yy1,yy2,zz11,xx01-bbxd,yy01+bbyd) &

-gg1(xx1,xx2,yy1,yy2,zz11a,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz11-zz46,xx01-bbxd,yy01+bbyd) &

+gg1(xx1,xx2,yy1,yy2,zz11a-zz46,xx01-bbxd,yy01) ! *0.7

ssb=gg1(xx1,xx2,yy1,yy2,zz11,xx01+bbxd,yy01+bbyd) &

-gg1(xx1,xx2,yy1,yy2,zz11a,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz11-zz46,xx01+bbxd,yy01+bbyd) &

+gg1(xx1,xx2,yy1,yy2,zz11a-zz46,xx01+bbxd,yy01) ! *0.7

ssa=gg1(xx1,xx2,yy1,yy2,zz11a,xx01-bbxd,yy01+bbyd) &

-gg1(xx1,xx2,yy1,yy2,zz447,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz11a-zz46,xx01-bbxd,yy01+bbyd) &

+gg1(xx1,xx2,yy1,yy2,zz447-zz46,xx01-bbxd,yy01) ! *0.7

ssab=gg1(xx1,xx2,yy1,yy2,zz11a,xx01+bbxd,yy01+bbyd) &

-gg1(xx1,xx2,yy1,yy2,zz447,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz11a-zz46,xx01+bbxd,yy01+bbyd) &

+gg1(xx1,xx2,yy1,yy2,zz447-zz46,xx01+bbxd,yy01) ! *0.7

ss4=gg1(xx1,xx2,yy1,yy2,zz447,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz447a,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz447-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz447a-zz46,xx01-bbxd,yy01) ! *0.7

ss4b=gg1(xx1,xx2,yy1,yy2,zz447,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz447a,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz447-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz447a-zz46,xx01+bbxd,yy01) ! *0.7

ss4a=gg1(xx1,xx2,yy1,yy2,zz447a,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz441,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz447a-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz441-zz46,xx01-bbxd,yy01) ! *0.7

ss4ab=gg1(xx1,xx2,yy1,yy2,zz447a,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz441,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz447a-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz441-zz46,xx01+bbxd,yy01) ! *0.7

ss1=gg1(xx1,xx2,yy1,yy2,zz441,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz441a,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz441-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz441a-zz46,xx01-bbxd,yy01) ! *0.7

ss1b=gg1(xx1,xx2,yy1,yy2,zz441,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz441a,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz441-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz441a-zz46,xx01+bbxd,yy01) ! *0.7

ss1a=gg1(xx1,xx2,yy1,yy2,zz441a,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz44a,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz441a-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz44a-zz46,xx01-bbxd,yy01) ! *0.7

ss1ab=gg1(xx1,xx2,yy1,yy2,zz441a,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz44a,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz441a-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz44a-zz46,xx01+bbxd,yy01) ! *0.7

ss2=gg1(xx1,xx2,yy1,yy2,zz44a,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz44ab,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz44a-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz44ab-zz46,xx01-bbxd,yy01) ! *0.7

ss2b=gg1(xx1,xx2,yy1,yy2,zz44a,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,zz44ab,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz44a-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,zz44ab-zz46,xx01+bbxd,yy01) ! *0.7

ss2a=gg1(xx1,xx2,yy1,yy2,zz44ab,xx01-bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,z46,xx01-bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz44ab-zz46,xx01-bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,z46-zz46,xx01-bbxd,yy01) ! *0.7

ss2ab=gg1(xx1,xx2,yy1,yy2,zz44ab,xx01+bbxd,yy01+bbyd) & !V

-gg1(xx1,xx2,yy1,yy2,z46,xx01+bbxd,yy01) & ! *0.7 &

-gg1(xx1,xx2,yy1,yy2,zz44ab-zz46,xx01+bbxd,yy01+bbyd) & !V

+gg1(xx1,xx2,yy1,yy2,z46-zz46,xx01+bbxd,yy01) ! *0.7

! ss2=gxzg(xx1,xx2,yy1,yy2,zz11,zz11a,xx01,yy01+bbyd) ! & !V

! -gxzg(xx1,xx2,yy1,yy2,zz11,zz11a,xx01-bbxd,yy01) *0.7

! ss2a=gyzg(xx1,xx2,yy1,yy2,zz11,zz11a,xx01,yy01+bbyd) ! & !V

! -gyzg(xx1,xx2,yy1,yy2,zz11,zz11a,xx01-bbxd,yy01) *0.7

j3=j2+20*(j1-1);

xyzz(i,j,j3)=ssb-ss*DLux ! /6.672

xyzz(i,j,j3+400)=ssab-ssa*DLux ! /6.672

xyzz(i,j,j3+800)=ss4b-ss4*DLux ! /6.672

xyzz(i,j,j3+1200)=ss4ab-ss4a*DLux ! /6.672

xyzz(i,j,j3+1600)=ss1b-ss1*DLux ! /6.672

xyzz(i,j,j3+2000)=ss1ab-ss1a*DLux ! /6.672

xyzz(i,j,j3+2400)=ss2b -ss2*DLux ! /6.672

xyzz(i,j,j3+2800)=ss2ab-ss2a*DLux ! /6.672

enddo; enddo; enddo; enddo

do i=1,20

do j=1,20

sv(i,j)=0.001

enddo; enddo

print*, ' to print POLE GRAV'

open(602,file='d:\work\temp\testgo602.dat',status='unknown',form='formatted')

write(602,33) gtolerg

rewind 602

read(602,33) gtole !Гравика вариации (1800Р-600р)м

rewind 602

print '(1x,10f7.2)',gtole

close(602)

!43 format(2209f10.5)

do i5=1, nx7

do j5=1,ny7

i=i5+30;j=j5+120

!ZecsigTg(i5,j5)=gtolerg(i,j)

enddo ;end do

open(603,file='d:\work\temp\testgo603.dat',status='unknown',form='formatted')

!write(603,43) gtole

!rewind 603

!read(603,43) gtole !Гравика (Х0;У0)=(+10;+10)

!rewind 603

!print '(1x,10f7.2)',gtole

close(603)

!43 format(2209f10.5)

m=3200

do j3=1,m ! {nachalo ITERACIY}

lambdai(j3)=0.0

do i8=1, nx7

do j8=1, ny7

lambdai(j3)=lambdai(j3)+ (xyzz(i8,j8,j3))

enddo ;enddo;enddo

do j3=1,m

!lambdai(j3)=dsqrt(lambdai(j3)) !Т³льки для В2 - проба

enddo

do i3=1,nx7

do j3=1, ny7

lambda(i3,j3)=0.0

do i8=1,m

lambda(i3,j3)=lambda(i3,j3)+ (xyzz(i3,j3,i8))

enddo ;enddo ;enddo

Nbjo=nx7*ny7

do j3=1,m ! {nachalo ITERACIY}

lbdai(j3)=0.0

do i8=1, nx7

do j8=1, ny7

lbdai(j3)=lbdai(j3)+ (xyzz(i8,j8,j3))*lambda(i8,j8)

enddo ;enddo;enddo

do i8=1,m

sum5=0.0

sum6=0.0

sum7=0.0

sum8=0.0

su1(i8)=0.0

su(i8)=0.0

enddo

do i8=1,m

do i3=1,nx7

do j3=1, ny7

!su1(i8)=su1(i8)+ gtole(i3,j3)*xyzz(i3,j3,i8) &

!/lambda(i3,j3)/lambdai(i8)

!su(i8)=su(i8)+ xyzz(i3,j3,i8)/lambda(i3,j3)/lambdai(i8)

su1(i8)=su1(i8)+ gtole(i3,j3)*xyzz(i3,j3,i8)

enddo ;enddo

sum5=sum5+ lambdai(i8)/lbdai(i8)

sum6=sum6+ (lambdai(i8)/lbdai(i8))**2

sum7=sum7+ su1(i8)/lbdai(i8)

sum8=sum8+ lambdai(i8)/lbdai(i8)**2*su1(i8)

enddo

ajo=0.0

sum3=0.0

sum4=0.0

bjo=0.0

! {*******************}

do i12=1,m

sum3=sum3+su(i12)

sum4=sum4+su1(i12)

ajo=ajo+su1(i12)*su(i12)

bjo=bjo+su(i12)*su(i12) !tauB

enddo

Nbjo=m !Для Li1 !Для Li0.5 отменить

ajo=ajo/Nbjo

bjo=bjo/Nbjo

sum3=sum3/Nbjo

sum4=sum4/Nbjo

sum1=0.0

sum2=0.0

do i3=1,nx7

do j3=1, ny7

! sum1=sum1+lambda(i3,j3)

!sum2=sum2+lambda(i3,j3)*lambda(i3,j3)

enddo ;enddo

!print*, ' to print (SL1 , SL2 ) '

print*,sum1 ,sum2

SRgtL=0.0

SRgt=0.0

do i5=1, nx7

do j5=1,ny7

!SRgt=SRgt+ gtole(i5,j5) ! SRpoLe

! SRgtL=SRgtL+ gtole(i5,j5)*lambda(i5,j5)

enddo ;end do

print*, ' to print (SRgt, SRgtL ) '

! print*,SRgt ,SRgtL

! SRgtB=(SRgt*sum2-sum1*srgtL)/(sum2*nx7*ny7-sum1*sum1) !A0

! SRgtR=(SRgtL-SRgt*sum1/nx7/ny7)/(sum2-sum1*sum1/nx7/ny7) !A1

print*, ' to print (A0, A1 ) '

! print '(10x,2f15.5)', SRgtB,SRgtR

! sum1=sum1/Nbjo

!SRgtB=(SRgt/Nbjo -sum1*ajo)/(1-sum1*bjo) !A0

! SRgtR=(ajo-SRgt*bjo/Nbjo)/(1-sum1*bjo) !A1

! print '(10x,2f15.5)', SRgtB,SRgtR

! SRgtB=(ajo -sum3*sum4)/(bjo-sum3*sum3) !A0

! SRgtR=(sum4*bjo-sum3*ajo)/(bjo-sum3*sum3) !A1

! print '(10x,2f15.5)', SRgtB,SRgtR

SRgtB=(sum8 -sum5*sum7/m)/(sum6-sum5*sum5/m) !A0

SRgtR=(sum7*sum6/m-sum5*sum8/m)/(sum6-sum5*sum5/m) !A1

print '(10x,2f15.5)', SRgtB,SRgtR

write(*,*) '06 08 2011 22 введите a,n'

read (*,*)a,n

write(*,*) 'a=',a,' n=',n

do j=1,20

do i=1,20

j3=i+20*(j-1)

! vecsig(j3)=0.0050

! vecsig(j3+400)=0.00530 !sv(j,i)

! vecsig(j3+800)=0.00560 !sv(j,i)

! vecsig(j3+1200)=0.0060 !sv(j,i)

enddo ; enddo

print*, ' to print 118 '

open(118,file='d:\work\temp\testgo118.dat',status='unknown',form='formatted')

write(118,34) ffonty3

rewind 118

!read(118,43) Ztole2

!rewind 118

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

close(118)

do i=1, nx7

do j=1,ny7/2+1

f1ffonty3(i,j)=f1fonty3(i,2*j-1)

enddo ; enddo

do i=1, nx7

do j=1,ny7/2

f2ffonty3(i,j)=f2fonty3(i,2*j)

enddo ; enddo

print*, ' to print 119 f1'

open(119,file='d:\work\temp\testgo119.dat',status='unknown',form='formatted')

write(119,45) f1ffonty3

rewind 119

!read(118,43) Ztole2

!rewind 118

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

close(119)

print*, ' to print 120 f2'

open(120,file='d:\work\temp\testgo120.dat',status='unknown',form='formatted')

write(120,43) f2ffonty3

rewind 120

!read(118,43) Ztole2

!rewind 118

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

close(120)

write(*,*) '06 08 2011 22 введите a,n'

read (*,*)a,n

write(*,*) 'a=',a,' n=',n

do i=1, nx7

do j20=1,ny7

xypi=0.0

do j3=801,m

! xypi=xypi+xyzz(i,j20,j3)*Vecsig(j3) !{vybrannoe pole} !!! 11

enddo

!gtole(i,j20)=gtole(i,j20)-xypi ! {vybrannoe pole}

enddo ; enddo

open(131,file='d:\work\temp\testgo131.dat',status='unknown',form='formatted')

!write(131,43) gtole !Osredn.pole

!rewind 131

!read(132,43) gtole !a

!rewind 132

!print '(1x,10f7.2)',gtole

close(131)

goto 236

do j3=1,400

vecsiga(j3)=vecsig(j3)

vecsiga(j3+400)=vecsig(j3+400)+vecsiga(j3)*0 !sv(j,i)

vecsiga(j3+800)=vecsig(j3+800)+vecsiga(j3+400)*0 !sv(j,i)

vecsiga(j3+1200)=vecsig(j3+1200)+vecsiga(j3+800)*0 !sv(j,i)

vecsiga(j3+1600)=vecsig(j3+1600)+vecsiga(j3+1200)*0

vecsiga(j3+2000)=vecsig(j3+2000)+vecsiga(j3+1600)*0 !sv(j,i)

vecsiga(j3+2400)=vecsig(j3+2400)+vecsiga(j3+2000)*0 !sv(j,i)

vecsiga(j3+2800)=vecsig(j3+2800)+vecsiga(j3+2400)*0 !sv(j,i)

enddo

236 jo=1

! goto 236

do in=1, 20

do iL=1,20

VSL(iL,in)=0.0

VSLa(iL,in)=0.0

enddo;end do

do in=1, 20

do ik=1,8

do iL=1,20

VSL(iL+20*(ik-1+1),in)=vecsig(iL+20*(in-1)+400*(ik-1)) !+SRgtr

enddo;end do

enddo

do k=1,20

do j3=1, 160

enddo;end do

print*, ' to print 741 '

k=6

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsL(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

!print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 741 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 741 '

open(741,file='d:\work\temp\testgo741.dat',status='unknown',form='formatted')

write(741,23) svL2

rewind 741

read(741,23) svL2

rewind 741

!print '(1x,8f7.3)',svL2

close(741)

print*, ' to print 741 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

print*, ' to print 742 '

k=8

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsL(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

! print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 742 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 742 '

open(742,file='d:\work\temp\testgo742.dat',status='unknown',form='formatted')

write(742,23) svL2

rewind 742

read(742,23) svL2

rewind 742

!print '(1x,8f7.3)',svL2

close(742)

print*, ' to print 742 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

k=10

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsL(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

!print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 743 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 743 '

open(743,file='d:\work\temp\testgo743.dat',status='unknown',form='formatted')

write(743,23) svL2

rewind 743

read(743,23) svL2

rewind 743

!print '(1x,8f7.3)',svL2

close(743)

print*, ' to print 743 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

! 236 jo=1

! goto 954

do in=1, 20

do ik=1,8

do iL=1,20

VSLa(iL+20*(ik-1+1),in)=vecsig(iL+20*(in-1)+400*(ik-1)) !+SRgtr

enddo;end do

enddo

do k=1,20

do j3=1, 160

enddo;end do

print*, ' to print 745 '

k=12

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsLa(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

!print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 745 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 745 '

open(745,file='d:\work\temp\testgo745.dat',status='unknown',form='formatted')

write(745,23) svL2

rewind 745

!read(745,23) svL2

!rewind 745

!print '(1x,8f7.3)',svL2

close(745)

print*, ' to print 745 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

print*, ' to print 746 '

k=14

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsLa(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

! print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 746 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 746 '

open(746,file='d:\work\temp\testgo746.dat',status='unknown',form='formatted')

write(746,23) svL2

rewind 746

!read(746,23) svL2

!rewind 746

!print '(1x,8f7.3)',svL2

close(746)

print*, ' to print 746 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

k=16

do j3=1,180; j=(j3-1)/20+1;i=j3-20*(j-1);svL2(10-j,i)=vsLa(j3,k) !v matlab plan ne rokirovatj

svL1(i,j)=svL2(10-j,i) ! v f90

enddo

!print '(1x,10f7.3)',((svL1(i,j),i=1,10),j=1,8)

print*, ' to print 747 '

! print '(1x,10f7.3)',((svL1(i,j),i=11,20),j=1,8)

print*, ' to print 747 '

open(747,file='d:\work\temp\testgo747.dat',status='unknown',form='formatted')

write(747,23) svL2

rewind 747

!read(747,23) svL2

!rewind 747

!print '(1x,8f7.3)',svL2

close(747)

print*, ' to print 747 '

!print '(1x,8f7.3)',((svL2(i,j),i=1,8),j=1,20)

954 jo=1

! 236 jo=1

do ii=1,m

d(ii)=vecsig(ii) !+SRgtr

d1(ii)=vecsig(ii)

enddo

print*, ' to print 716 '

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(716,file='d:\work\temp\testgo716.dat',status='unknown',form='formatted')

write(716,22) sv

rewind 716

!read(716,22) sv

!rewind 716

!print '(1x,10f7.3)',sv

close(716)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

print*, ' to print 736 '

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(736,file='d:\work\temp\testgo736.dat',status='unknown',form='formatted')

write(736,22) sv

rewind 736

!read(716,22) sv

!rewind 716

!print '(1x,10f7.3)',sv

close(736)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+800) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 349 '

! print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(349,file='d:\work\temp\testgo349.dat',status='unknown',form='formatted')

write(349,22) sv

rewind 349

!read(349,22) sv

!rewind 349

!print '(1x,10f7.3)',sv

close(349)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+800) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 738 '

! print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(738,file='d:\work\temp\testgo738.dat',status='unknown',form='formatted')

write(738,22) sv

rewind 738

!read(349,22) sv

!rewind 349

!print '(1x,10f7.3)',sv

close(738)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+1600) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 348 '

print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(348,file='d:\work\temp\testgo348.dat',status='unknown',form='formatted')

write(348,22) sv

rewind 348

!read(348,22) sv

!rewind 348

!print '(1x,10f7.3)',sv

close(348)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+1600) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 740 '

print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(740,file='d:\work\temp\testgo740.dat',status='unknown',form='formatted')

write(740,22) sv

rewind 740

!read(348,22) sv

!rewind 348

!print '(1x,10f7.3)',sv

close(740)

!print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+2400) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 519 '

! print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(519,file='d:\work\temp\testgo519.dat',status='unknown',access='sequential',form='formatted')

write(519,22) sv

rewind 519

! read(519,22) sv

! rewind 519

close(519)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+2400) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 752 '

! print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(752,file='d:\work\temp\testgo752.dat',status='unknown',access='sequential',form='formatted')

write(752,22) sv

rewind 752

! read(519,22) sv

! rewind 519

close(752)

! print '(1x,10f7.3)',((sv(i,j),i=1,10),j=1,10)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+400) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 520 '

print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(520,file='d:\work\temp\testgo520.dat',status='unknown',access='sequential',form='formatted')

write(520,22) sv

rewind 520

!read(519,22) sv

!rewind 519

close(520)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+400) !v matlab plan ne rokirovatj

sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 737 '

print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(737,file='d:\work\temp\testgo737.dat',status='unknown',access='sequential',form='formatted')

write(737,22) sv

rewind 737

!read(519,22) sv

!rewind 519

close(737)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+1200) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 521 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(521,file='d:\work\temp\testgo521.dat',status='unknown',access='sequential',form='formatted')

write(521,22) sv

rewind 521

!read(531,22) sv

!rewind 531

close(521)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+1200) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 739 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(739,file='d:\work\temp\testgo739.dat',status='unknown',access='sequential',form='formatted')

write(739,22) sv

rewind 739

!read(531,22) sv

!rewind 531

close(739)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+2000) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 522 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(522,file='d:\work\temp\testgo522.dat',status='unknown',access='sequential',form='formatted')

write(522,22) sv

rewind 522

!read(519,22) sv

!rewind 519

close(522)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+2000) !v matlab plan ne rokirovatj ! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 751 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(751,file='d:\work\temp\testgo751.dat',status='unknown',access='sequential',form='formatted')

write(751,22) sv

rewind 751

!read(519,22) sv

!rewind 519

close(751)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+2800) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 523 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(523,file='d:\work\temp\testgo523.dat',status='unknown',access='sequential',form='formatted')

write(523,22) sv

rewind 523

!read(531,22) sv

!rewind 531

close(523)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d1(j3+2800) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 753 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(753,file='d:\work\temp\testgo753.dat',status='unknown',access='sequential',form='formatted')

write(753,22) sv

rewind 753

!read(531,22) sv

!rewind 531

close(753)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+2400) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 524 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(524,file='d:\work\temp\testgo524.dat',status='unknown',access='sequential',form='formatted')

write(524,22) sv

rewind 524

!read(519,22) sv

!rewind 519

close(524)

do j3=1,400; j=(j3-1)/20+1;i=j3-20*(j-1);sv(j,i)=d(j3+2800) !v matlab plan ne rokirovatj

! sv1(i,21-j)=sv(j,i) ! v f90

enddo

print*, ' to print 525 '

! print '(1x,10f7.3)',((sv1(i,j),i=1,10),j=1,10)

print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(525,file='d:\work\temp\testgo525.dat',status='unknown',access='sequential',form='formatted')

write(525,22) sv

rewind 525

!read(531,22) sv

!rewind 531

close(525)

print*, ' A697 '

do ii=1,m

!d(ii)=vecsig(ii) - vecsiga(ii)

d(ii)=vecsig(ii) !+SRgtr

enddo

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3200) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 526 '

print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

!print '(1x,10f7.2)',((sv1(i,j),i=11,20),j=11,20)

open(526,file='d:\work\temp\testgo526.dat',status='unknown',access='sequential',form='formatted')

write(526,44) svo

rewind 526

!read(519,22) sv

!rewind 519

close(526)

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3300) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 527 '

! print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(527,file='d:\work\temp\testgo527.dat',status='unknown',access='sequential',form='formatted')

write(527,44) svo

rewind 527

!read(531,22) sv

!rewind 531

close(527)

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3400) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 528 '

print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(528,file='d:\work\temp\testgo528.dat',status='unknown',access='sequential',form='formatted')

write(528,44) svo

rewind 528

!read(519,22) sv

!rewind 519

close(528)

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3500) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 529 '

!print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(529,file='d:\work\temp\testgo529.dat',status='unknown',access='sequential',form='formatted')

write(529,44) svo

rewind 529

!read(531,22) sv

!rewind 531

close(529)

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3600) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 530 '

print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(530,file='d:\work\temp\testgo530.dat',status='unknown',access='sequential',form='formatted')

write(530,44) svo

rewind 530

!read(519,22) sv

!rewind 519

close(530)

do j3=1,100; j=(j3-1)/10+1;i=j3-10*(j-1);svo(j,i)=d(j3+3700) !v matlab plan ne rokirovatj

sv1(i,11-j)=svo(j,i) ! v f90

enddo

print*, ' to print 531 '

!print '(1x,10f7.2)',((sv1(i,j),i=1,10),j=1,10)

! print '(1x,10f7.3)',((sv1(i,j),i=11,20),j=11,20)

open(531,file='d:\work\temp\testgo531.dat',status='unknown',access='sequential',form='formatted')

write(531,44) svo

rewind 531

!read(531,22) sv

!rewind 531

close(531)

! do i=1, nx7

!do j20=1,ny7

!xypi=0.0

!do j3=1,m

! xypi=xypi+xyzz(i,j20,j3)*Vecsig(j3) !{vybrannoe pole} !!! 11

! enddo

! ffonty1(i,j20)=ffonty1(i,j20)-xypi ! {vybrannoe pole}

! enddo ;

! enddo

! ! print '(1x,5d12.3)',((ffonty1(i,j),i=24,33),j=1,5)

! print '(1x,5d12.3)',((ffonty1(i,j),i=15,19),j=16,18)

! write(*,*)' OSTATOK Rj NOVOGO PolJA '

! open(118,file='d:\work\temp\testgo118.dat',status='unknown',form='formatted')

!write(118,43) ffonty1

!rewind 118

!read(118,43) Ztole2

!rewind 118

!close(118)

!print*, ' to print ffonty1 '

!print '(1x,10f7.3)',((ffonty1(i,j),i=15,24),j=1,31)

!sum=0

! do i=1, nx7

! do j20=1,ny7

!sum=sum+ dabs(Ffonty1(i,j20))

! enddo ; enddo

!print*, ' to print DS(ffonty1) '

! sum=(sum/nx7/ny7)

! print '(1x,2f10.4)',sum

write(*,*) '02 09 2012 Ferma введите a,n'

read (*,*)a,n

write(*,*) 'a=',a,' n=',n

Call date_and_time(date=dat,time=tim,zone=zon)

print*,dat,tim,' TIME ',zon

GOTO 653

N=50

M=5

k=0

do i=2, N

ai=i

do j=i+1, N+1

aj=j

do i13=j+1,j+i+1

ai13=i13

desa=ai13**M-aj**M-ai**M

if (desa==0) then

k=k+1

print *,k, desa, i,j,i13

END IF

enddo;enddo ;enddo

print*, ' to print DeSa Ferma '

print *,k

653 io6=1

Call date_and_time(date=dat,time=tim,zone=zon)

print*,dat,tim,' TIME ',zon

write(*,*) '02 09 2012 Ferma введите a,n'

read (*,*)a,n

write(*,*) 'a=',a,' n=',n

CONTAINS

function gg (x1,x2,y1,y2,z1,z2,x0,y0);

real(8)::gg

real(8):: xsq111,xsq211,xsq121,xsq221,xsq112,xsq212,xsq122, &

xsq222,ysq111,ysq211,ysq121,x1,x2,y1,y2,z1,z2,x0,y0

real(8):: ysq221,ysq112,ysq212,ysq122,ysq222,t1a,t2a,ty2ln, &

ty1ln,ttx2ln,ttx1ln

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ112=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ212=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ122=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

XSQ222=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ112=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

YSQ122=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

YSQ212=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

YSQ222=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

t1a=z1*(datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1))

t2a=z2*(datan((xsq222+y2-y0)/z2)+datan((xsq112+y1-y0)/z2)- &

datan((xsq122+y2-y0)/z2)-datan((xsq212+y1-y0)/z2))

ty2ln=dlog(dabs(xsq221*xsq122/xsq121/xsq222))*(y2-y0)

ty1ln=dlog(dabs(xsq211*xsq112/xsq111/xsq212))*(y1-y0)

ttx2ln=dlog(dabs(ysq221*ysq212/ysq211/ysq222))*(x2-x0)

ttx1ln=dlog(dabs(ysq121*ysq112/ysq111/ysq122))*(x1-x0)

gg=(2.0*(t1a-t2a)+ty2ln-ty1ln+ttx2ln-ttx1ln)*6.672

end function gg

function gxzg (x1,x2,y1,y2,z1,z2,x0,y0);

real(8)::gxzg

real(8):: ysq111,ysq211,ysq121,x1,x2,y1,y2,z1,z2,x0,y0

real(8):: ysq221,ysq112,ysq212,ysq122,ysq222

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ112=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

YSQ122=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

YSQ212=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

YSQ222=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

gxzg=dlog(dabs((ysq222/ysq122)*(ysq121/ysq221)* &

(ysq211/ysq111)*(ysq112/ysq212)))*6.672

end function gxzg

function gyzg (x1,x2,y1,y2,z1,z2,x0,y0); !250

real(8)::gyzg

real(8):: xsq111,xsq211,xsq121,xsq221,xsq112,xsq212,xsq122, &

xsq222,x1,x2,y1,y2,z1,z2,x0,y0

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ112=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ212=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ122=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

XSQ222=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

gyzg=dlog(dabs((xsq222/xsq122)*(xsq121/xsq221)* &

(xsq211/xsq111)*(xsq112/xsq212)))*6.672

end function gyzg

function gzzg (x1,x2,y1,y2,z1,z2,x0,y0);

real(8)::gzzg

real(8):: xsq111,xsq211,xsq121,xsq221,xsq112,xsq212,xsq122, &

xsq222,x1,x2,y1,y2,z1,z2,x0,y0

real(8)::t1a,t2a

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ112=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ212=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ122=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

XSQ222=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

t1a=(datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1))

t2a=(datan((xsq222+y2-y0)/z2)+datan((xsq112+y1-y0)/z2)- &

datan((xsq122+y2-y0)/z2)-datan((xsq212+y1-y0)/z2))

gzzg=(2.0*(t1a-t2a))*6.672

end function gzzg

function gg1 (x1,x2,y1,y2,z1,x0,y0); !250

real(8)::gg1

real(8):: xsq111,xsq211,xsq121,xsq221,ysq111,ysq211,ysq121, &

x1,x2,y1,y2,z1,x0,y0,ysq221,t1a,ty1ln,ttx1ln,ty2ln,ttx2ln

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t1a=z1*(datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1))

ty2ln=dlog(dabs(xsq221/xsq121))*(y2-y0)

ty1ln=dlog(dabs(xsq111/xsq211))*(y1-y0)

ttx2ln=dlog(dabs(ysq221/ysq211))*(x2-x0)

ttx1ln=dlog(dabs(ysq111/ysq121))*(x1-x0)

gg1=(2.0*t1a+ty2ln+ty1ln+ttx2ln+ttx1ln)*6.672 !280

end function gg1

function gg12 (x1,x2,y1,y2,z1,x0,y0); !250

real(8)::gg12

real(8):: xsq111,xsq211,xsq121,xsq221,ysq111,ysq211,ysq121,ysq221, &

x1,x2,y1,y2,z1,x0,y0,t1a

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t1a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

gg12=-(2.0*t1a)*6.672 !439poslednee

end function gg12

function gg13 (x1,x2,y1,y2,z1,x0,y0); !250

real(8)::gg13

real(8):: xsq111,xsq211,xsq121,xsq221,ysq111,ysq211,ysq121,ysq221, &

x1,x2,y1,y2,z1,x0,y0,t1a,t2a

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t1a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

z1=z1+0.0001

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t2a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

gg13=-(200.0*(t2a-t1a))*6.672*100.0 !439poslednee

end function gg13

function gg14 (x1,x2,y1,y2,z1,x0,y0); !250

real(8)::gg14

real(8):: xsq111,xsq211,xsq121,xsq221,ysq111,ysq211,ysq121,ysq221, &

x1,x2,y1,y2,z1,x0,y0,t1a,t2a,t3a

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t1a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

z1=z1+0.001

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t2a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

z1=z1+0.001

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

t3a=datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1)

gg14=-(2000.0*(t3a-2*t2a+t1a))*6.672*1000.0 !439poslednee !

end function gg14

function ggz1 (x1,x2,y1,y2,z1,z2,x0,y0); !250

real(8)::ggz1

real(8):: xsq111,xsq211,xsq121,xsq221,xsq112,xsq212,xsq122, &

xsq222,ysq111,ysq211,ysq121,x1,x2,y1,y2,z1,z2,x0,y0

real(8):: ysq221,ysq112,ysq212,ysq122,ysq222,t1a,t2a,tya, &

txa,Lozxy

XSQ111=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ211=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

XSQ121=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ221=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

XSQ112=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2) !260

XSQ212=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2)

XSQ122=x1-x0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

XSQ222=x2-x0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

YSQ111=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ121=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ211=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z1*z1)

YSQ221=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z1*z1)

YSQ112=y1-y0+dsqrt((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0)+z2*z2)

YSQ122=y2-y0+dsqrt((x1-x0)*(x1-x0)+(y2-y0)*(y2-y0)+z2*z2)

YSQ212=y1-y0+dsqrt((x2-x0)*(x2-x0)+(y1-y0)*(y1-y0)+z2*z2) !270

YSQ222=y2-y0+dsqrt((x2-x0)*(x2-x0)+(y2-y0)*(y2-y0)+z2*z2)

t1a=z1*z1*(datan((xsq221+y2-y0)/z1)+datan((xsq111+y1-y0)/z1)- &

datan((xsq121+y2-y0)/z1)-datan((xsq211+y1-y0)/z1))

t2a=z2*z2*(datan((xsq222+y2-y0)/z2)+datan((xsq112+y1-y0)/z2)- &

datan((xsq122+y2-y0)/z2)-datan((xsq212+y1-y0)/z2))

tya=(y2-y0)*(y2-y0)*(datan((y2-y0)/(xsq222+z2))-datan((y2-y0)/(xsq122+z2)) &

-datan((y2-y0)/(xsq221+z1))+datan((y2-y0)/(xsq121+z1))) - &

(y1-y0)*(y1-y0)*(datan((y1-y0)/(xsq212+z2))-datan((y1-y0)/(xsq112+z2)) &

-datan((y1-y0)/(xsq211+z1))+datan((y1-y0)/(xsq111+z1)))

txa=(x2-x0)*(x2-x0)*(datan((x2-x0)/(ysq222+z2))-datan((x2-x0)/(ysq212+z2)) &

-datan((x2-x0)/(ysq221+z1))+datan((x2-x0)/(ysq211+z1)))- &

(x1-x0)*(x1-x0)*(datan((x1-x0)/(ysq122+z2))-datan((x1-x0)/(ysq112+z2)) &

-datan((x1-x0)/(ysq121+z1))+datan((x1-x0)/(ysq111+z1)))

Lozxy=(x2-x0)*(y2-y0)*dlog(dabs((xsq222-(x2-x0)+z2)/(xsq221-(x2-x0)+z1)))- &

(x2-x0)*(y1-y0)*dlog(dabs((xsq212-(x2-x0)+z2)/(xsq211-(x2-x0)+z1)))- &

(x1-x0)*(y2-y0)*dlog(dabs((xsq122-(x1-x0)+z2)/(xsq121-(x1-x0)+z1)))+ &

(x1-x0)*(y1-y0)*dlog(dabs((xsq112-(x1-x0)+z2)/(xsq111-(x1-x0)+z1)))

ggz1=(Lozxy-txa-tya-t2a+t1a)*6.672

end function ggz1

end