f=inline('816*x^3-3835*x^2+6000*x-3125');

ezplot(f,1.43,1.71);
hold on;
plot(5/3,0,'*',25/16,0,'o',25/17,0,'x');
Title('Zeros of 816x^3-3835x^2+6000x-3125');
ylabel('P(x)');
%axis([.75, 2.25, -2, 2]);
grid;