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

Necran

a=3
b=8
m=20

function [y]=funkcjaF1(x)
y=log(x)+1
end

function [y]=funkcjaF2(x)
y=1/x
end

function[si1,si2]=metoda_simpsona(a,b,m)

n=2*m
h=(b-a)/n
si1=0
si2=0

for i=1:2 Sad

2*m-1)
x=a+i*h;
si1=si1+4*funkcjaF1(x); %dla funkcji 1.
si2=si2+4*funkcjaF2(x); %dla funkcji 2.
end

for i=2:2 Sad

2*m-2)
x=a+i*h;
si1=si1+2*funkcjaF1(x); %dla funkcji 1.
si2=si2+2*funkcjaF2(x); %dla funkcji 2.
end

si1=1/3*h*(funkcjaF1(a)+funkcjaF1(b)+si1); %dla funkcji 1.
si2=1/3*h*(funkcjaF1(a)+funkcjaF1(b)+si2); %dla funkcji 2.

end %metoda_simpsona

[pole1,pole2]=metoda_simpsona(a,b,m)
pole_wlasciwe=pole1-pole2

algebroman

Fajnie wyszło to w pętlach Smile

Wystarczyłoby napisać for i=1:2:2*n-1

czasowstrzymywacz

Strasznie to powalone...