> eq:=2*y^5-4*y^3+2*y-2*y^2+1=0;
>
> solve(eq,y);
>
> sols := [solve(eq,y)];
> evalf({%});
>
> plot(2*y^5-4*y^3+2*y-2*y^2+1,y=1..1.5);
> eq:=y^5-2*y^3+y-4*y^2+2=0;
> solve(eq,y);
> evalf({%});
> plot(y^5-2*y^3+y-4*y^2+2,y=1..2);
>
>