2014 dxdy logo

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

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




Начать новую тему Ответить на тему
 
 Maple долго строит график
Сообщение17.05.2023, 19:56 
Аватара пользователя


21/09/13
136
Уфа
Доброго времени суток!
Хотел посмотреть, как выглядит график функции
$$
f(x)=\int\limits_0^x \frac{dt}{t\ln^2(t)\sqrt{x-t}}.
$$
Функция $f$ для малых положительных $x$ определена корректно, поскольку функция $(t \ln^2t)^{-1}$ интегрируема, например, на $[0,1/2]$.
Пытаюсь выполнить в Maple 2022 такую команду
  1. h := evalf(Int(1/(t*ln(t)^2*(x - t)^(1/2)), t = 0 .. x, epsilon = 0.001)) 
  2. plot(h, x = 0 .. 1/2) 

Но Maple только снизу в левом углу пишет "Вычисление...".
Точность вычисления стоит очень низкая, но даже это не помогает, хотя кажется, что график должен выглядить очень просто.

 Профиль  
                  
 
 Re: Maple долго строит график
Сообщение17.05.2023, 23:51 
Заслуженный участник


12/07/07
4452
Ваш вариант быстро выполнился в Maple 7.
В Maple 15 у меня сработал вариант без evalf:
Код:
> h := Int(1/(t*ln(t)^2*(x - t)^(1/2)), t = 0 .. x, epsilon = 0.001):
> plot(h, x=0..1/2, thickness=2, color=black);

Вложение:
Maple15_pic.PNG
Maple15_pic.PNG [ 8.9 Кб | Просмотров: 340 ]


-- Wed 17.05.2023 23:03:33 --

В Maple 2021 ситуация аналогична Maple 15: без evalf быстро рисует, а с evalf - отображает (в нижнем левом углу) "Вычисление...".

 Профиль  
                  
 
 Re: Maple долго строит график
Сообщение18.05.2023, 10:42 
Аватара пользователя


21/09/13
136
Уфа
Большое спасибо. Без evalf действительно вычисляет. Я вчера еще добавил команду Digits:=5, и он тоже отобразил график (evalf вроде бы был на месте). Хотя я этот evalf специально для Maple'а добавил, по моему представлению это должно было ему облегчить вычисления :-)

 Профиль  
                  
 
 Re: Maple долго строит график
Сообщение18.05.2023, 11:30 
Заслуженный участник


12/07/07
4452
У меня с Digits:= 5 в Maple 15 также сработало. Но менять точность вычислений без необходимости так себе идея.
Даже со значением Digits по умолчанию без evalf c epsilon = 0.00001 Maple 15 строит картинку. (Дальше лень было проверять.)

Напомню себе.
Int(, epsilon = value) — это инертная форма численного вычисления интеграла (инертная — без самого вычисления).
evalf(Int(, epsilon = value)) — это численное вычисление интеграла.
plot() «вычисляет» инертные формы.
Например
Код:
> h := Int(t^2, t = 0..x):  plot(h, x=0..1);
Построит график.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 4 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group