function rotation,object,vector,angle

matrix=dblarr(3,3)

x=vector[0]/mag(vector)
y=vector[1]/mag(vector)
z=vector[2]/mag(vector)

matrix[0,0]=cos(angle)+(1-cos(angle))*x^2
matrix[1,0]=(1-cos(angle))*x*y-sin(angle)*z
matrix[2,0]=(1-cos(angle))*x*z-sin(angle)*y
matrix[0,1]=(1-cos(angle))*y*x+sin(angle)*z
matrix[1,1]=cos(angle)+(1-cos(angle))*y^2
matrix[2,1]=(1-cos(angle))*y*z-sin(angle)*x
matrix[0,2]=(1-cos(angle))*z*x-sin(angle)*y
matrix[1,2]=(1-cos(angle))*z*y+sin(angle)*x
matrix[2,2]=cos(angle)+(1-cos(angle))*z^2

return,matrix##object

end