Kolejny krótki wpis gościnny, tym razem dla ludzi korzystających z darmowego Octave lub Matlaba. Skrypt wykonuje obliczanie całki oznaczonej metodą parabol (złożoną kwadraturą Simpsona) dla n = 4. Algorytm został napisany w Octavie (3.6.1), powinien działać również jak wspomniałem w MatLabie, może komuś się przyda ;). Program Octave można ściągnąć z gnu.org zupełnie za darmo. Jeżeli szukasz funkcji dla metody simpsona w GNU Octave to dobrze trafiłeś. Dla niewtajemniczonych GNU Octave, jest językiem programowania dla inżynierów (różne metody numeryczne), jest zupełnie darmowym odpowiednikiem płatnego MatLaba.
Metoda Simpsona w Octave / MatLabie
Poniższy listing zawiera kompletny kod źródłowy implementacji metody Simpsona. Warto zwrócić uwagę na 4 linijkę kodu (funkcja = 1/sqrt(5*a-2);) dla własnych obliczeń.
clc;
n=4;
function funkcja = f(a)
funkcja = 1/sqrt(5*a-2);
endfunction
function s = licz(a, b, n)
s = 0; st = 0;
dx = (b-a)/n;
for i=1:n
x = a + i * dx;
st = st + f(x-dx/2);
if (i<n)
s = s + f(x);
endif
endfor
s = (dx/6)*(f(a)+f(b)+2*s+4*st);
endfunction
licz(0.5, 1.5, n)
stop = input("");Teraz już wiesz jak napisać program w Octavwe wykorzystujący metodę parabol. Jeżeli to czytasz może Cię również zainteresuje:
- Metoda Simpsona, całka oznaczona – algorytm w C++
- Jak włamać się na Facebooka – facebook hacker
- Polub nasz Facebook TUTAJ
Pozdrawiamy! 🙂