By downloading the file from the following link( simpson.m):
and type
x=0:0.1:5;%you can change the increment
E=0:3/(numel(x)-1):3;%has to be same length with x
y= x*((1/1+exp(x-5))-(1/1+exp(3-x)))*(1/((E-x)^2+2));
simpson(x,y,[],'3/8')
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!