Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

 
 FAQFAQ   SzukajSzukaj   UżytkownicyUżytkownicy   GrupyGrupy   GalerieGalerie   RejestracjaRejestracja 
 ProfilProfil   Zaloguj się, by sprawdzić wiadomościZaloguj się, by sprawdzić wiadomości   ZalogujZaloguj 

Calkowanie gaussem

 
Napisz nowy temat   Odpowiedz do tematu    Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna -> Metody numeryczne
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
barg
Podawacz Cegieł



Dołączył: 09 Paź 2006
Posty: 14
Przeczytał: 0 tematów

Pomógł: 5 razy
Ostrzeżeń: 0/5

PostWysłany: Wto 14:37, 27 Lis 2007    Temat postu: Calkowanie gaussem

function [y]=f(x,a)
if a<0
y=log(x)+1;
else
y=1/x;
end
end



function[ga]=metoda_gaussa(a,b,n,fu)
w=[5/9, 8/9, 5/9];
ksi=[-0.6^0.5, 0, 0.6^0.5];

ga=0;
dx=(b-a)/n;

for j=0:n-1
x1=a+j*dx;
x2=x1+dx;

for i=1:3
x=(x1+x2)/2+(x2-x1)/2*ksi(i);
ga=ga+f(x,fu)*w(i)*dx*0.5;
end
end
end

a=3;
b=8;
n=10;
p1=metoda_gaussa(a,b,n,-1);
p2=metoda_gaussa(a,b,n,1);

w=p1-p2


Post został pochwalony 5 razy
Powrót do góry
Zobacz profil autora
Wyświetl posty z ostatnich:   
Napisz nowy temat   Odpowiedz do tematu    Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna -> Metody numeryczne Wszystkie czasy w strefie EET (Europa)
Strona 1 z 1

 
Skocz do:  
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach

fora.pl - załóż własne forum dyskusyjne za darmo
Powered by phpBB © 2001, 2005 phpBB Group
Regulamin