x=0:0.01:5; 

y=func(x); 

plot(x,y) 

1:3:18 

plot(x,y,'*') 

newton(3,5) 

newton(3,4) 

newton(3,4) 

newton(3,40)