MATLABで台形近似

N=100;

min=0;
max=1;
t=linspace(min,max,N);


dt=t(2)-t(1);

S=zeros(size(t));

for i=2:1:length(t)
    S(i)=S(i-1)+((func(t(i))+func(t(i-1)))*dt)/2;
end

figure 
plot(t,S)
grid on

function y=func(x)
    y=2*x^4;
end