using simpson 3/8 rule

 F(x)= x*((1/1+exp(x-5))-(1/1+exp(3-x)))*(1/((E-x)^2+2))
where  x limit from 0 to 5 nm
and E vary from 0 to 3 
can anyone tell me matlab code for this  to solve this integration using simpson 3/8 rule

Answer 1

[Birdman]

By downloading the file from the following link( simpson.m):

https://www.mathworks.com/matlabcentral/fileexchange/33493-simpson-s-1-3-and-3-8-rules

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')