% Plot a function to be tested for integration

i=0.0:.01:3.0;
y=30+2.*i.^3-25.*sin(i);

plot(i,y)
xlabel('x')
ylabel('f(x)')
ylim([0 90])
grid
