barg
Podawacz Cegieł
Dołączył: 09 Paź 2006
Posty: 14
Przeczytał: 0 tematów
Pomógł: 5 razy Ostrzeżeń: 0/5
|
Wysł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
|
|