Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

algebroman

funkcja całkująca:
function sr=metoda_prostokatow(a,b,n);
h=(b-a)/n;
sr=0;
for i=0:n-1
xl=a+i*h;
xp=xl+h;
sr=sr+fun((xp+xl)/2);
end%i
sr=sr*h;
endfunction;
wartością, którą funkcja zwraca jest "sr"-nie sprawdzałem czy działa Wink

miejsce przecięcia regułą falsi nie wiem jak znaleźć, trzeba porównać dwie funkcje ale nie wiem jak znaleźć ich punkty przecięcia jeżeli są 2 bo numerycznie można znaleźć chyba jedno m. zerowe w danym przedziale. więcej na [link widoczny dla zalogowanych]
ew. algorytm z wykładu:
function [x,iter]=metoda_regula_falsi(a,b,eps)
fa=f(a);
fb=f(b);
if fa*fb>0
printf("błąd-(a,b)-nie jest przedziałem izolacji");
for iter=1:150
x=a-fa/(fb-fa)*(b-a);
fx=f(x);
if fa*fx>0
if fabs(x-a)<eps && fabs(fx)<10^(-10)
return();
a=x;
fa=fx;
else
if fabs(x-b)<eps && fabs(fx)<10^(-10)
return();
b=x;
fb=fx;
end%iter
printf("błąd-proces rozbieżny");
return();
endfunction
^to jest przepisane z wykładu, nie wiem czy działa
have fun:)