Вычислить площадь фигуры, ограниченной кривой

f=@(x)sin(log(x));%анонимная функция 
%построение графика для выбора интервала интегрирования 
h=ezplot(f,[-pi,pi]); 
grid on,hold on 
E = 1e-7;% заданная точность 
 
%выбираем как будем закрашивать фигуру 
figure1 = input('Введите левую границу Фигуры='); 
figure2 = input('Введите правую границу Фигуры= '); 
disp('Найдем координаты х, в которых функция пересекает ось абсцисс: '); 
n = input('Введите число корней уравнения n = '); 
disp('Нашли координаты х, в которых функция пересекает осью абсцисс: '); 
for i=1:n 
 
a = input('Введите левую границу отрезка, содержащего корень, a = '); 
b = input('Введите правую границу отрезка, содержащего корень, b = '); 
syms x 
R = solve(f) 
end 
title('y=sin(log(x))'); %добавление названия графика 
 
%РАСЧЕТ 
%разбиваем на 2 
%считаем левую и правую площадь через рунге(симпсон) 
%ГРАФИК и расчет ЛЕВОЙ ФИГУРЫ 
f=@(x)sin(log(x)); 
E = 1e-7;% заданная точность 
a=0.04; 
b=1.0; 
disp('Площадь левой фигуры plochadLevoy'); 
plochadLevoy = Runge(f,a,b,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadLevoy ='); 
disp(abs(plochadLevoy)); 
 
%закрашиваем фигуру 
x=a:0.0001:b; 
plot(a,0.04, b,0.04) 
fill(x,f(x),'m-') 
title('y=sin(log(x))') 
 
%ГРАФИК и расчет ПРАВОЙ ФИГУРЫ 
f=@(x)sin(log(x)); 
E = 1e-7;% заданная точность 
a=1.0; 
b=2.1; 
disp('Площадь правой фигуры plochadPravoy'); 
plochadPravoy = Runge(f,a,b,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadPravoy'); 
disp(abs(plochadPravoy)); 
%закрашиваем фигуру 
a=1.0; 
b=2.1; 
x=1.0:0.1:b; 
fill(x,0,'r') 
%y=0:0.1:2.1; 
%fill(y,f(y),'r') 
title('y=sin(log(x))') 
 
%найдем общую площадь 
disp('Площадь общей фигуры ЛЕВАЯ ПЛОЩАДЬ+ПРАВАЯ'); 
Obchaya=plochadPravoy-plochadLevoy

f=@(x)sin(log(x));%анонимная функция 
%построение графика для выбора интервала интегрирования 
h=ezplot(f,[-pi,pi]); 
grid on,hold on 
eps = 1e-7;% заданная точность 
%пЕРВАЯ ТОЧКА ПЕРЕСЕЧЕНИЯ 
%A= 
 
 
 
b=solve(f);%2 точка пересечения с осью ОХ 
%выбираем как будем закрашивать фигуру 
figure2 = input('Введите правую границу Фигуры= '); 
b=solve(f);%2 точка пересечения с осью ОХ 
%РАСЧЕТ 
%разбиваем на 2 
%считаем левую и правую площадь через рунге(симпсон) 
%ГРАФИК и расчет ЛЕВОЙ ФИГУРЫ 
disp('Площадь левой фигуры plochadLevoy'); 
E=1e-7; 
plochadLevoy = Runge(f,a,b,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadLevoy ='); 
disp(abs(plochadLevoy)); 
%закрашиваем фигуру 
x=a:0.01:b; 
 
fill(x,f(x),'m-') 
%ГРАФИК и расчет ПРАВОЙ ФИГУРЫ 
disp('Площадь правой фигуры plochadPravoy'); 
plochadPravoy = Runge(f,b,figure2,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadPravoy'); 
disp(abs(plochadPravoy)); 
x=b:0.01:figure2;%вектор значений х 
x1=[x figure2];% 
y1=[f(x) 0]; 
fill(x1,y1,'r') 
title('y=sin(log(x))') 
%найдем общую площадь 
disp('Площадь общей фигуры ЛЕВАЯ ПЛОЩАДЬ+ПРАВАЯ'); 
Obchaya=abs(abs(plochadPravoy)+abs(plochadLevoy ))

f=@(x)sin(log(x));%анонимная функция 
%построение графика для выбора интервала интегрирования 
h=ezplot(f,[-pi,pi]); 
grid on,hold on 
eps = 1e-7;% заданная точность 
%пЕРВАЯ ТОЧКА ПЕРЕСЕЧЕНИЯ 
%A= 
n = input('Введите число корней уравнения n = '); 
for i=1:n 
 
a = input('Введите левую границу отрезка, содержащего корень, a = '); 
b = input('Введите правую границу отрезка, содержащего корень, b = '); 
 
[ R(i) k(i) ] = Newton(f,eps,a,b) 
 
end 
disp('Нашли координаты х, в которых функция пересекает осью абсцисс: '); 
R 
A=R 
TOCHKA1PERESECHENIYA=double(R) 
 
%выбираем как будем закрашивать фигуру 
figure2 = input('Введите правую границу Фигуры= '); 
TOCHKA2PERESECHENIYA=solve(f);%2 точка пересечения с осью ОХ 
TOCHKA2PERESECHENIYA2=double(TOCHKA2PERESECHENIYA) 
%РАСЧЕТ 
%разбиваем на 2 
%считаем левую и правую площадь через рунге(симпсон) 
%ГРАФИК и расчет ЛЕВОЙ ФИГУРЫ 
disp('Площадь левой фигуры plochadLevoy'); 
E=1e-7; 
plochadLevoy = Runge(f,a,b,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadLevoy ='); 
disp(abs(plochadLevoy)); 
%закрашиваем фигуру 
x=a:0.01:b; 
y=f(x); 
y1(:,1)=[]; 
ynew=[0 y1]; 
fill(ynew,f(x),'m-') 
%ГРАФИК и расчет ПРАВОЙ ФИГУРЫ 
disp('Площадь правой фигуры plochadPravoy'); 
plochadPravoy = Runge(f,b,figure2,E); %обращение к m-функции Runge 
disp('Приближенное значение опр. интеграла (м-д Симпсона с применение правила Рунге)'); 
disp('plochadPravoy'); 
disp(abs(plochadPravoy)); 
x=b:0.01:figure2;%вектор значений х 
x1=[x figure2];% 
y1=[f(x) 0]; 
fill(x1,y1,'r') 
title('y=sin(log(x))') 
%найдем общую площадь 
disp('Площадь общей фигуры ЛЕВАЯ ПЛОЩАДЬ+ПРАВАЯ'); 
Obchaya=abs(abs(plochadPravoy)+abs(plochadLevoy ))