program statvett1

integer::i,n
real,allocatable,dimension(:)::x,y
real:: xmedia,ymedia

read(5,*)n
allocate (x(n),y(n))
do i=1,n
   read(5,*) x(i),y(i)
enddo
xmedia=vmed(x,n)
ymedia=vmed(y,n)

write(6,*)'xmedia=',xmedia
write(6,*)'ymedia=',ymedia
write(6,*)'xvar=',varianza(x,xmedia,n)
write(6,*)'yvar=',varianza(y,ymedia,n)
write(6,*)'xmax=',vmax(x,n)
write(6,*)'ymax=',vmax(y,n)
write(6,*)'xmin=',vmin(x,n)
write(6,*)'ymin=',vmin(y,n)
end program statvett1

real function vmed(v,n)
integer:: n
real ,dimension(1)::v

real::somma

somma=0
do i=1,n
   somma=somma+v(i)
enddo
vmed=somma/n
end function vmed