Metoda Simpsona – kod w Octave / MatLab

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:

Pozdrawiamy! 🙂

Dodaj komentarz

Twój adres email nie zostanie opublikowany. Pola, których wypełnienie jest wymagane, są oznaczone symbolem *