Program 6 : afix, defvec

 
      program prog6
      implicit none

      real rms
      external rms
      real vect(13),vect2(13)
      integer i
      
      call afix(vect,2.34,13)
      do i=1,13
         vect2(i)=2.34
      enddo
      print *,'The error of afix is:',rms(vect,vect2,13)

      call defvec(vect,13,2.37,3.02)
      vect2(1)=3.02
      do i=2,13
         vect2(i)=vect2(i-1)+2.37
      enddo
      print *,'The error of defvec is:',rms(vect,vect2,13)

      stop 
      end

      real function rms(a,b,n)
      real a(n),b(n),sum
      integer i,n
      sum=0
      do i=1,n
         sum=sum+(a(i)-b(i))**2
      enddo
      rms=sqrt(sum/n)
      return
      end