diff --git a/Igor Experiment/Test new seteps vs original.h5xp b/Igor Experiment/Test new seteps vs original.h5xp
index bbadf8c05da0492583e8892cf5e0f3fb5a0e2410..4495edb895791d9d584f8fecd9f88aa91303ba66 100644
Binary files a/Igor Experiment/Test new seteps vs original.h5xp and b/Igor Experiment/Test new seteps vs original.h5xp differ
diff --git a/source/f90/seteps.f90 b/source/f90/seteps.f90
index 20337d9a6662dc9f0b3612b811537c6e90dc24f6..1d797da0b610d57874c61030b4ed0919d8b9ee7f 100644
--- a/source/f90/seteps.f90
+++ b/source/f90/seteps.f90
@@ -60,8 +60,8 @@ subroutine seteps(neps, nos, oscType, osc, epsinf, wn, eps, layers)
           addDrude = addDrude - wp**2 / dcmplx(wn2, y * wn)
         case (extendedDrude) ! Gervais 2002, equ. 2.14
           wp = osc(1,k)
-          yp = osc(2,k)
-          y0 = osc(3,k)
+          y0 = osc(2,k)
+          yp = osc(3,k)
           addDrude = addDrude - dcmplx(wp**2, (y0 - yp) * wn) / dcmplx(wn2, wn * y0)
         case default
           write (*,'(A)') '*** Error in oscillator type! ***'