% Get error bound for HW31

clear
close all

x = linspace(0,3,101)';
y = hw3der(x);
plot(x,y)

my = max(y);
h = 2.727272e-1;
bound = h^4/24*my