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

, ,

>> syms x a b;

>> f=sym(1/(x^2+1));

>> g=sym(1/(x^2+a^2));

>> h=sym(1/((x+a)^2+b^2));

>> I=int(f,x)

I =

atan(x)

>> I=int(g,x)

I =

atan(x/a)/a

>> I=int(h,x)

I =

atan((a + x)/b)/b

=y

ex=y+-sqrt(1+y2) , т.к. ex>0, то ex=y+sqrt(1+y2) => x=ln(x+sqrt(1+x2)). Заменяем x->y и y->x :

arsh(x)=ln(x+sqrt(1+x2)).

Ln(x+sqrt(a2+x2)=arsh(x/a)+ln(a).

function int_left=s_l(f,a,b,n)

syms x;

delta_x=(b-a)/n;

s1=subs(f,x,a);

for i=1:n-1

s_l=s1+subs(f,x,a+delta_x*i)*delta_x;

end

disp(s_l)

end