program main
implicit none
integer i
real::x(i),y(i),sumx=0,sumy=0,xa,ya,b,a,sx=0,sy=0,sxy=0,f,rxy
data(x(i),i=1,20)/0.9,1.2,2.2,2.4,-0.5,2.5,-1.1,0,6.2,2.7,3.2,-1.1,2.5,1.2,1.8,0.6,2.4,2.5,1.2,-0.8
data(y(i),i=1,20)/32,25,20,26,27,24,28,24,15,16,24,30,22,30,24,33,26,20,32,35
do i=1,20
    read(8,*)y(i),x(i)
end do
do i=1,20
    sumx=sumx+x(i)
    xumy=xumy+y(i)
end do
xa=sumx/20
ya=sumy/20
do i=1,20
    sx=sx+(x(i)-xa)**2
    sy=sy+(y(i)-ya)**2
    sxy=sxy+(x(i)-xa)*(y(i)-ya)
end do
sx=sqrt(sx/20)
sy=sqrt(sy/20)
sxy=sxy/20
b=sxy/sx/sx
a=ya-b*xa
print*,'y=',b,'x','+',a 
rxy=xa*b/sy
f=rxy**2*18/(1-rxy**2)
if(f>4.41)write(*,*)'F=',f,'>Fα=4.41，回归方程显著'
if(f<=4.41)write(*,*)'F=',f,'<Fα=4.41，回归方程不显著'
end