Задание
.
Получение
решения в пакете MatLab.
Ход
работы
function
f11_1
a=input('enter
<a>: ');
b=input('enter
<b>: ');
n=input('enter
<n>: ');
I_rectav=RectAv(a,b,n)
I_trap=Trap(a,b,n)
X=a:(b-a)/n:b;
Y=f(X);
I_trapz=trapz(X,Y
)
I_quad=quad(@f,a,b)
function
I_rectav=RectAv(a,b,n)
h=(b-a)/n;
s=0;
for
i=1:n
x=a+i*h-h/2;
s=s+f(x);
end
I_rectav=s*h;
function
I_trap=Trap(a,b,n)
h=(b-a)/n;
s=f(a)/2+f(b)/2;
for
i=1:(n-1)
s=s+f(a+i*h);
end
I_trap=s*h;
function
t=f(x)
t=3*sin(x-4)-x.^5+5;
Результаты
работы программы.
f11_1
enter
<a>: 1
enter
<b>: 1.5
enter
<n>: 10
I_rectav
=
2.037388196825291e-001
I_trap
=
1.975694391833556e-001
I_trapz
=
1.975694391833556e-001
I_quad
=
2.016825235074213e-001
f11_1
enter
<a>: 1
enter
<b>: 1.5
enter
<n>: 100
I_rectav
=
2.017030921501479e-001
I_trap
=
2.016413861689008e-001
I_trapz
=
2.016413861689007e-001
I_quad
=
2.016825235074213e-001
f11_1
enter
<a>: 1
enter
<b>: 1.5
enter
<n>: 500
I_rectav
=
2.016833462540952e-001
I_trap
=
2.016808780101234e-001
I_trapz
=
2.016808780101232e-001
I_quad
=
2.016825235074213e-001
Вывод.
Программа
MatLab
работает правильно.
Список
литературы.
Учебное
пособие. «Информатика пакет MatLab»
СПб Изд. Политехнического университета
2010
22