program enshu04cc05                                                         
implicitnone
integer::seireki,joutai,nensuu
real,dimension(12)::kion
real::kion_nen_heikin,kion_nen_souwa
real::kion_nen_bunsan,kion_nen_hensa
real::kion2_nen_souwa,kion2_nen_heikin
real,dimension(12)::kion_tsuki_heikin,kion_tsuki_souwa
real,dimension(12)::kion2_tsuki_heikin,kion2_tsuki_souwa
real,dimension(12)::kion_tsuki_bunsan,kion_tsuki_hensa
real,dimension(12)::seki_souwa,seki_heikin,kaiki_a,kaiki_b,r2
real::seireki_heikin,seireki2_heikin,seireki_bunsan
integer::seireki_souwa,seireki2_souwa
integer::i
doi=1,12
kion_tsuki_souwa(i)=0.0
kion2_tsuki_souwa(i)=0.0
enddo
seki_souwa(5)=0.0
nensuu=0
seireki_souwa=0
seireki2_souwa=0
do
read(*,'(i4,12f5.1)',iostat=joutai)seireki,kion
if(joutai/=0)exit
nensuu=nensuu+1
seireki_souwa=seireki_souwa+seireki
seireki2_souwa=seireki2_souwa+seireki**2
kion_nen_souwa=0.0
kion2_nen_souwa=0.0
doi=1,12
kion_nen_souwa=kion_nen_souwa+kion(i)
kion2_nen_souwa=kion2_nen_souwa+kion(i)**2
kion_tsuki_souwa(i)=kion_tsuki_souwa(i)+kion(i)
kion2_tsuki_souwa(i)=kion2_tsuki_souwa(i)+kion(i)**2
enddo
seki_souwa(5)=seki_souwa(5)+real(seireki)*kion(5)
kion_nen_heikin=kion_nen_souwa/12.0
kion2_nen_heikin=kion2_nen_souwa/12.0
kion_nen_bunsan=kion2_nen_heikin-kion_nen_heikin**2
kion_nen_hensa=sqrt(kion_nen_bunsan)
write(*,'(i4,13f6.1,2f6.2)')seireki,kion,kion_nen_heikin&
,kion_nen_bunsan,kion_nen_hensa
enddo
doi=1,12
kion_tsuki_heikin(i)=kion_tsuki_souwa(i)/real(nensuu)
kion2_tsuki_heikin(i)=kion2_tsuki_souwa(i)/real(nensuu)
kion_tsuki_bunsan(i)=kion2_tsuki_heikin(i)-kion_tsuki_heikin(i)**2
kion_tsuki_hensa(i)=sqrt(kion_tsuki_bunsan(i))
enddo
write(*,'(i4,12f6.2)')nensuu,kion_tsuki_heikin
write(*,'(i4,12f6.2)')nensuu,kion_tsuki_bunsan
write(*,'(i4,12f6.2)')nensuu,kion_tsuki_hensa
seireki_heikin=real(seireki_souwa)/real(nensuu)
seireki2_heikin=real(seireki2_souwa)/real(nensuu)
seireki_bunsan=seireki2_heikin-seireki_heikin**2
seki_heikin(5)=seki_souwa(5)/real(nensuu)
kaiki_a(5)=(seki_heikin(5)-kion_tsuki_heikin(5)*seireki_heikin)&
/seireki_bunsan
kaiki_b(5)=kion_tsuki_heikin(5)-kaiki_a(5)*seireki_heikin
r2(5)=kaiki_a(5)**2*seireki_bunsan/kion_tsuki_bunsan(5)
write(*,'(f8.4,f6.1,f7.3)')kaiki_a(5),kaiki_b(5),r2(5)
endprogramenshu04cc05