anchor.diff:        do n1=1,15
            
            ! perturb x^S:
  -         call random_number( LinVec_1 )
  -         LinVec_1 = LinVec_1/real(n1**3)               
  -         
  -         LinVec_2_theta = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
  +         call getRandomVector(LinVec_1)
  +         
  +         LinVec_2_approx = matmul( map_mod1_mod2%dm%mi, LinVec_1 )  ! approximate delta theta^D         
                         
            j=1
            do i=1,Mesh1_O%nnodes
  
positive.code:             j = j+3
           end do
                    
           print *, 'Translational Acceleration: ', n1, TwoNorm( LinVec_2 - LinVec_2_theta )
                                
        end do      
        
  
negative.code:             j = j+3
           end do
                    
           print *, 'Translational Displacement: ', n1, TwoNorm( LinVec_2 - LinVec_2_approx )
                                
        end do      
        
  