Ollie - I think that you almost have the correct implementation for Simpson's algorithm, but you may want to reconsider how you are determining h. Note that the question is telling you what h should be. So is your calculation giving you the same value? You may want to pass h into your function and calculate h from that. Like the following
func = @(x) 2*x + 2*sin(x) - 4;
mysimp(func, 3.1,3,0.1)
and your function mysimp is defined as
function simp = mysimp(fName,x2,x1,h)
m = floor((x2-x1)/h);
% etc.
Note how we calculate m from the inputs and so have m subintervals. Now notice how you are creating the array of x values
x = x1:h:x2;
A problem with this is that you are not guaranteed that the last value in this array will be x2 (it all depends on x1 and the step size h). Instead, you should use the linspace function as
x = linspace(x1,x2,m);
to create an array of m+1 elements starting from x1 and ending at x2 (the step size between each value will be near enough to h). We then calculate the y values as
y = fName(x);
