Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 2 ] 
Автор Сообщение
 Заголовок сообщения: Построение графиков
СообщениеДобавлено: 17 май 2021, 00:41 
Не в сети
Продвинутый
Зарегистрирован:
24 дек 2019, 00:16
Сообщений: 58
Cпасибо сказано: 18
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Доброго времени суток) нужна помощь с графиками в maple.
Код:
restart;
with(plots):
tau1:=mu:
tau2:=mu/5:
mu:=1:
r0:=1:
#M(lambda,psi,beta0)
H_M:=beta0:
R2_M:=-(2*lambda)/(beta0^2*psi^2)*ln(1-beta0*psi*r/lambda)-2*r/(beta0*psi);
dRdr2_M:=r^2/(R2_M*(lambda-psi*r*H_M)^2);
I1_M:=dRdr2_M+R2_Q/r^2+H_M^2+psi^2*R2_M+lambda^2;
I2_M:=(lambda-psi*r*H_M)^2*R2_M/r^2+dRdr2_M*(R2_M/r^2+H_M^2+psi^2*R2_M+lambda^2);
M:=unapply(Pi*(int((tau1*diff(I1_M,psi)+tau2*diff(I2_M,psi))*r, r=0..r0)), lambda, psi, beta0):
#plot(M(1,psi,1),psi=-0.5..0.5, discont=true, thickness=2,color=col[i]);
#M(lambda,psi,0)
H_Mbeta0:=0:
R2_M0:=r^2/lambda;
dRdr2_M0:=r^2/(R2_M0*(lambda-psi*r*H_Mbeta0)^2);
I1_M0:=dRdr2_M0+R2_M0/r^2+H_Mbeta0^2+psi^2*R2_M0+lambda^2;
I2_M0:=(lambda-psi*r*H_Mbeta0)^2*R2_M0/r^2+dRdr2_M0*(R2_M0/r^2+H_Mbeta0^2+psi^2*R2_M0+lambda^2);
M0:=unapply(Pi*(int((tau1*diff(I1_M0,psi)+tau2*diff(I2_M0,psi))*r, r=0..r0)), lambda, psi, beta0);
display( plot(M(1,psi,-2),psi=-0.5..0.5,discont=true, thickness=2 ), plot(M(1,psi,-1),psi=-0.5..0.5, discont=true, thickness=2 ), plot(M0(1,psi,0), psi=-0.5..0.5, thickness=2), plot(M(1,psi,1),psi=-0.5..0.5,discont=true, thickness=2 ), plot(M(1,psi,2),psi=-0.5..0.5,discont=true, thickness=2 ));

Выполнялось более 4 часов, графики так и не построились (ошибки не выводились, выполнение прервала). Помогите, пожалуйста, это исправить.

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Построение графиков
СообщениеДобавлено: 17 май 2021, 21:42 
Не в сети
Продвинутый
Зарегистрирован:
14 фев 2020, 10:46
Сообщений: 92
Cпасибо сказано: 12
Спасибо получено:
42 раз в 38 сообщениях
Очков репутации: 26

Добавить очки репутацииУменьшить очки репутации
Основная проблема в том, что в строке
M0:=unapply(Pi*(int((tau1*diff(I1_M0,psi)+tau2*diff(I2_M0,psi))*r, r=0..r0)), lambda, psi, beta0);

вы пытаетесь символьно вычислить весьма сложный интеграл ещё и с кучей параметров. Maple просто зависает на этом вычислении. Я заменил явное вычисление на инертную форму. При построении графика Maple считает этот интеграл численно:

restart;
with(plots):
tau1:=mu:
tau2:=mu/5:
mu:=1:
r0:=1:
H_M:=beta0:
R2_M:=-(2*lambda)/(beta0^2*psi^2)*ln(1-beta0*psi*r/lambda)-2*r/(beta0*psi):
dRdr2_M:=r^2/(R2_M*(lambda-psi*r*H_M)^2):
I1_M:=dRdr2_M+R2_Q/r^2+H_M^2+psi^2*R2_M+lambda^2:
I2_M:=(lambda-psi*r*H_M)^2*R2_M/r^2+dRdr2_M*(R2_M/r^2+H_M^2+psi^2*R2_M+lambda^2):
M:=unapply(eval(Pi*Int((tau1*diff(I1_M,psi)+tau2*diff(I2_M,psi))*r, r=0..r0),[lambda=1, beta0=1]),psi);

plot(M,-0.5..0.5, thickness=2, numpoints=40);

Вернуться к началу
 Профиль  
Cпасибо сказано 
За это сообщение пользователю Kitonum "Спасибо" сказали:
Susanna Gaybaryan
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему      Страница 1 из 1 [ Сообщений: 2 ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Построение графиков

в форуме Начала анализа и Другие разделы школьной математики

alinochka

10

2383

12 июн 2014, 22:23

Построение графиков

в форуме Тригонометрия

351w

4

448

04 апр 2018, 13:36

Построение графиков в maxima

в форуме Информатика и Компьютерные науки

ANazarov1980

1

700

05 апр 2018, 23:30

Построение графиков функции

в форуме Пределы числовых последовательностей и функций, Исследования функций

DodDod

3

543

04 апр 2018, 13:52

Построение нескольких графиков

в форуме MathCad

kivikivi_777

0

530

05 май 2014, 16:37

Построение графиков функций

в форуме Пределы числовых последовательностей и функций, Исследования функций

angelaapina

4

577

09 янв 2015, 19:27

Разложение в ряд Фурье и построение графиков

в форуме Ряды Фурье и Интегральные преобразования

351w

15

1075

09 апр 2019, 08:00

Исследование функций и построение графиков

в форуме Пределы числовых последовательностей и функций, Исследования функций

Xaron

4

370

23 дек 2016, 19:25

Построение совмещенных графиков по периоду

в форуме Maple

AlanSbor

7

552

09 апр 2021, 12:09

Построение графиков по результатам расчета в Excel

в форуме MATLAB

URL2000

0

291

05 мар 2020, 20:54


Часовой пояс: UTC + 3 часа [ Летнее время ]



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

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


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

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved