функции и подпрограммы фортран

[Igor1997]

помогите переделать код чтобы присутствовали функции или подпрограммы

program lab4 !апроксимация данных
implicit none
real,dimension(:),allocatable:: func,arg
real::a,b,avgy,avgx,sum1,sum2
integer::i,kol

open(unit=1,file='dataXX.dat')
Print*,'vsego tochek?'
Read(*,*) kol

allocate(func(1:kol))
allocate(arg(1:kol))

do i=1,kol
Read(1,*) arg(i),func(i)
end do

sum1=0
sum2=0

do i=1,kol
sum1=sum1+func(i)
sum2=sum2+arg(i)
end do

avgy=sum1/kol
avgx=sum2/kol

sum1=0
sum2=0

do i=1,kol
sum1=sum1+(arg(i)-avgx)*(func(i)-avgy)
sum2=sum2+(arg(i)-avgx)**2
end do

b=sum1/sum2
a=avgy-b*avgx

Print*, 'a=',a
Print*, 'b=',b

end program lab4