function rim=rotate1(im,angle)

r=(angle/360)*2*pi;

rotmat = [cos(r) sin(r); -sin(r) cos(r)];

[rows,cols]=size(im);
xcoords = 1:cols;
ycoords = 1:rows;
xcarray = repmat(xcoords,[rows 1]);
ycarray = repmat(ycoords,[1 cols]);
xcv = xcarray(:);
ycv = ycarray(:);

pairs = [xcv'; ycv'];

newpairs = rotmat*pairs;
newpairs = round(newpairs);

newpairs(newpairs<1)=1;
newpairs(newpairs>rows)=rows;

rim=zeros(rows,cols);

for i=1:size(newpairs,2)
  rim(newpairs(1,i),newpairs(2,i))=im(pairs(1,i),pairs(2,i));
end

