clear clc syms x close disp('Graficas de la funciones:') f=input(' Funcion a Graficar: '); P=input(' Intervalos: '); a=P(1); b=P(2); k=0; for i=a:0.01:b k=k+1; x=i; X(k)=x; unos(k)=0; Y(k)=eval(f); end hold on plot(X,Y) plot(X,unos) hold off