Numerical integration: Difference between revisions

(Added an ActionScript version.)
Line 1,099:
print integrate "i.trapezium "fn2 4 -1 2 ; 2.351014
print integrate "i.simpsons "fn2 4 -1 2 ; 2.447732</lang>
 
=={{header|MATLAB}}==
 
f.m (The funtion you are wanting to approximate an integral too)
<lang MATLAB>
function f=f(x)
f=x;
</lang>
 
Numerical_Integration.m
<lang MATLAB>
function Numerical_Integration(a,b,n)
format long;
h=(b-a)/n;
x1=linspace(a,b,n+1);
LftRecRl=h.*sum(f(x1(1:n)))
 
x2=linspace(a,b,n+1);
RgtRecRl=h.*sum(f(x2(2:n+1)))
 
x3=linspace(a,b,n+1);
MdPtRl=h.*sum(f((x3(1:n)+x3(2:n+1))/2))
 
x4=linspace(a,b,n+1);
TrpRl=h.*sum((f(x3(1:n))+f(x4(2:n+1)))/2)
 
x5=linspace(a,b,n+1);
SmpRl=h/6.*sum(f(x5(1:n))+4*f((x5(1:n)+x5(2:n+1))/2)+f(x5(2:n+1)))
</lang>
 
=={{header|OCaml}}==
Anonymous user