anchor.diff:      do i=3,(N-1)
         y(i,localt) = a1*(initY(i+2)+initY(i-2)) + a2*(initY(i+1) + initY(i-1)) &
    + a3*initY(i) + a4*initY(i) + a5*(initY(i+1)+initY(i-1)) &
  - + deltaT*deltaT*N*stringForce(i)/Ms
  + - deltaT*deltaT*N*stringForce(i)/Ms
       end do
   
   ! Transfer the bridge element string amplitudes to bridgeY:
  
positive.code:     do i=3,(N-1)
        y(i,localt) = a1*(y(i+2,localt-1)+y(i-2,localt-1)) + a2*(y(i+1,localt-1) + y(i-1,localt-1)) &
   + a3*y(i,localt-1) + a4*y(i,localt-2) + a5*(y(i+1,localt-2)+y(i-1,localt-2)) &
   + deltaT*deltaT*N*stringForce(i)/Ms
      end do
  
      if (hammerDone.eqv..FALSE.) then
  
negative.code:     do i=3,(N-1)
        y(i,localt) = a1*(y(i+2,localt-1)+y(i-2,localt-1)) + a2*(y(i+1,localt-1) + y(i-1,localt-1)) &
   + a3*y(i,localt-1) + a4*initY(i) + a5*(initY(i+1)+initY(i-1)) &
   - deltaT*deltaT*N*stringForce(i)/Ms
      end do
  
  ! Transfer the bridge element string amplitudes to bridgeY:
  