Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Otchet_po_laboratornym_rabotam_v_srede_programm...doc
Скачиваний:
1
Добавлен:
27.09.2019
Размер:
838.66 Кб
Скачать
  1. Задание

.

Получение решения в пакете MatLab.

  1. Ход работы

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;

  1. Результаты работы программы.

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 работает правильно.

Список литературы.

  1. Учебное пособие. «Информатика пакет MatLab» СПб Изд. Политехнического университета 2010

22

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]