Лабораторная работа: Интеграл Мат лаб (лаба 3)
Описание
Характеристики лабораторной работы
Список файлов
- Интеграл.txt 1,39 Kb
function Integralz;
format long;
% При P= 1 -используется формула средних прямоугольников;
% 2 -используется формула трапеций;
% 3 -используется формула Симпсона.
a=0;
b=1;
E=10^(-6);
n=1000;
f= @ (x) 2*x+exp(x);
Ih=0;
Ih2=0;
Rung=1;
disp('1 -используется формула средних прямоугольников')
disp('2 -используется формула трапеций')
disp('3 -используется формула Симпсона')
P = input('выбираем № = ')
if P==1
while Rung > E
Ih=0;
Ih2=0;
H=(b-a)/n;
h=H/2;
for i=0:2*n-1
x=a+h/2+h*i;
Ih2=Ih2+h*f(x);
end
h=H;
for i=0:n-1
x=a+h/2+h*i;
Ih=Ih+h*f(x);
end
Rung=abs((Ih2-Ih)/(2^2-1));
n=n*2;
end
rung=(Ih2-Ih)/(2^2-1);
end
if P==2
while Rung > E
Ih=0;
Ih2=0;
H=(b-a)/n;
h=H/2;
for i=0:2*n-1
x=a+h*i;
Ih2=Ih2+h*(f(x)+f(x+h))/2;
end
h=H;
for i=0:n-1
x=a+h*i;
Ih=Ih+h*(f(x)+f(x+h))/2;
end
Rung=abs((Ih2-Ih)/(2^2-1));
n=n*2;
end
rung=(Ih2-Ih)/(2^2-1);
end
if P==3
while Rung > E
Ih=0;
Ih2=0;
H=(b-a)/n;
h=H/2;
for i=0:n-1
x=a+2*h*i;
Ih2=Ih2+h*(f(x)+4*f(x+h)+f(x+2*h))/3;
end
h=H;
for i=0:n/2-1
x=a+2*h*i;
Ih=Ih+h*(f(x)+4*f(x+h)+f(x+2*h))/3;
end
Rung=abs((Ih2-Ih)/(2^4-1));
n=n*2;
end
rung=(Ih2-Ih)/(2^4-1);
end
Rung
Integral=Ih
Utochnen_po_Richardsony=Ih2+rung
n=n/2
h=H
h2=H/2
%Встроенная функция
if P==2
m = a:b/n:b;
y = log(10*m);
VSTROEN_F = trapz(m, y)
end
if P==3
VSTROEN_F = quad('log(10*x)', 0.1, 0.1*exp(1), 1e-6)
end