function [shepp] = SheppLogan(width)

n = -(width-1)/2;

for tap=1:width
  shepp(tap)=-2/(pi*(4*n^2-1));
  n=n+1;
end