function animate_impeller(times,sol_vals,omega)
figure
L = 1;
for i = 1:length(times)
   theta = omega*times(i);
   clf
   xlim([-1,1]);
   ylim([-1,1]);
   hold on
   plot([0,L*cos(theta)],[0,L*sin(theta)],'LineWidth',2)
   r = sol_vals(i,1);
   plot([r*cos(theta)],[r*sin(theta)],'ro','MarkerSize',11,'MarkerFaceColor','r');
   pause(0.15)
end

end