GNU Octave/Całkowanie numeryczne

Z Wikibooks, biblioteki wolnych podręczników.

Całkowanie numeryczne[edytuj]

Narysować wykres dystrybuanty 1-wymiarowego rozkładu normalnego, tj. wykres funkcji:

Całki w Octave oblicza się za pomocą quad, na przykład, żeby obliczyć możemy wykonać:

octave:2> quad("sin",0,pi) 
ans = 2
Dystrybuanta rozkładu normalnego .

Żeby obliczyć wartość funkcji F zdefiniujmy funkcję podcałkową:

function [y]=myexp(t)
   y=exp(-t*t/2)/sqrt(2*pi);
endfunction;

Następnie zdefiniujmy dystrybuantę (czyli funkcję F):

function [y]=dystrybuanta(x)
   y=quad("myexp", -inf, x);
endfunction;

Rysujemy wykres dystrybuanty F:

x=(-10:0.1:10);
N=length(x);
y=zeros(1,N);
for (k=1:N)
   y(k)=dystrybuanta(x(k));
endfor;
plot(x,y,"-r;Dystrybuanta N(0,1);");