for k=170:4:290
  eval(sprintf('!pnmrotate -90 forward%03d.ppm > forwardr%03d.ppm',k,k))
end%for

