function polyplot(P,I) % polyplot(P,I) % graph polynomial P on interval I global poly_p poly_p=P; fplot(@f,I); function y=f(x) global poly_p y=polyval(poly_p,x);