Numerical integration: Difference between revisions
Content added Content deleted
(Added an ActionScript version.) |
|||
Line 1,099: | Line 1,099: | ||
print integrate "i.trapezium "fn2 4 -1 2 ; 2.351014 |
print integrate "i.trapezium "fn2 4 -1 2 ; 2.351014 |
||
print integrate "i.simpsons "fn2 4 -1 2 ; 2.447732</lang> |
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}}== |
=={{header|OCaml}}== |