diff --git a/shiftingPburg.m b/shiftingPburg.m
new file mode 100644
index 0000000..fe32166
--- /dev/null
+++ b/shiftingPburg.m
@@ -0,0 +1,13 @@
+function [ EEG ] = shiftingPburg( signal,frequency,window,windowShift,NFFT,order)
+%shiftingPburg uses pburg to calculate the power spectrum
+%   Detailed explanation goes here
+    signalShift=frequency*windowShift;
+    signalWindow=frequency*window;
+    shiftProp=window/windowShift;
+    EEG=zeros([fix(floor(size(signal,1)/signalWindow)-1)*shiftProp+1,NFFT/2+1]);
+    for j=1:fix(floor(size(signal,1)/signalWindow-1)*shiftProp+1)
+        EEG(j,:)=pburg(signal((j-1)*signalShift+1:(j-1)*signalShift+signalWindow),order,NFFT);
+    end
+
+end
+