2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




На страницу Пред.  1, 2
 
 Re: Несобственные кратные интегралы в Maple
Сообщение12.04.2014, 17:34 
“Digits = 10”
Код:
> f:= 1/(x^4+y^4)^(1/3):
> I1:= int(f, [x=-1..1, y=-1..1]);
I1 := 3/4*(-3^(1/2)*Pi*2^(1/2)*GAMMA(2/3)*cos(5/12*Pi)+2*3^(1/2)*Pi*2^(1/2)*GAMMA(2/3)*sin(1/12*Pi)-16*hypergeom([1/12, 1/4, 1/3],[13/12, 5/4],-1)*GAMMA(11/12)*GAMMA(3/4)*sin(1/12*Pi)*cos(5/12*Pi))/sin(1/12*Pi)/cos(5/12*Pi)/GAMMA(11/12)/GAMMA(3/4)
> evalf(I1);
                     11.39721951
> eps:= 10^(-12);
> evalf(Int(f, [x=-1..-eps, y=-1..1])) + evalf(Int(f, [x=eps..1, y=-1..1]));
                     11.39721927
[Редактирование: один пропущенный ">" добавил]

 
 
 
 Re: Несобственные кратные интегралы в Maple
Сообщение12.04.2014, 17:37 
Красота. :) Спасибо.

 
 
 
 Re: Несобственные кратные интегралы в Maple
Сообщение12.04.2014, 17:57 
Упс. Я выше копировал и забыл удалить numeric. Зачеркнул.

И еще, на всякий случай: evalf(Int(xxx)) эквивалентен int(xxx, numeric).

 
 
 [ Сообщений: 18 ]  На страницу Пред.  1, 2


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group