Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
14
Добавлен:
05.06.2015
Размер:
14.43 Кб
Скачать

Нижняя:

function [sum] = darbu(f,n,a,b)

sh=(b-a)/n;

sum=0;

for i=0:n-1

sum=sum+sh*feval(f,(a+i*sh));

end

end

Верхняя:

function [sum] = darbu(f,n,a,b)

sh=(b-a)/n;

sum=0;

for i=1:n

sum=sum+sh*feval(f,(a+i*sh));

end

end

Нижняя

>> darbu(@(x) x+1,5,0,5)

ans =

15

>> darbu(@(x) x+1,10,0,5)

ans =

16.2500

>> darbu(@(x) x+1,100,0,5)

ans =

17.3750

>> darbu(@(x) x+1,1000,0,5)

ans =

17.4875

>> darbu(@(x) x^2,5,0,1)

ans =

0.2400

>> darbu(@(x) x^2,10,0,1)

ans =

0.2850

>> darbu(@(x) x^2,100,0,1)

ans =

0.3284

>> darbu(@(x) x^2,1000,0,1)

ans =

0.3328

Верхняя

>> darbu(@(x) x+1,5,0,5)

ans =

20

>> darbu(@(x) x+1,10,0,5)

ans =

18.7500

>> darbu(@(x) x+1,100,0,5)

ans =

17.6250

>> darbu(@(x) x+1,1000,0,5)

ans =

17.5125

>> darbu(@(x) x^2,5,0,1)

ans =

0.4400

>> darbu(@(x) x^2,10,0,1)

ans =

0.3850

>> darbu(@(x) x^2,100,0,1)

ans =

0.3384

>> darbu(@(x) x^2,1000,0,1)

ans =

0.3338