Модераторы: Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [Maple] Вложенные циклы 
:(
    Опции темы
PrettyGirl
Дата 23.10.2007, 17:42 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 62
Регистрация: 9.5.2006

Репутация: нет
Всего: нет
Привет всем! Помогите пожалуйста с Maple!
Вот такой ФРАГМЕНТ КОДА для 2-х уравнений работает
Код

for j from 1 to 3 do
> y[1,j] :=y[1,j-1] + h * subs(y1 = y[1,j-1],y2=y[2,j-1],f[1]):
> y[2,j] :=y[2,j-1] + h * subs(y1 = y[1,j-1],y2=y[2,j-1],f[2]): 
>  end do;

начала переделывать для n уравнений,не могу понять в чем ошибка.
Код

y[1,0] := 3;
y[2,0] := 1; 
f[1] := y1^2 / (3 * y2); 
f[2] := 6 * y2 / y1;
h := 0.01;                                                           
for i from 1 to 2 do                                                                                         
for j from 1 to 3 do
y[i,j] :=y[i,j-1] + h * subs(y1 = y[i,j-1],y2=y[i+1,j-1],f[i]): 
end do:
end do;
print(y);

помогите пожалуйста! 
PM MAIL   Вверх
kali
Дата 24.10.2007, 10:08 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 139
Регистрация: 9.11.2006
Где: Минск

Репутация: 18
Всего: 20
Ошибка похоже в end do

попробуй так

Код

y[1,0] := 3;
y[2,0] := 1; 
f[1] := y1^2 / (3 * y2); 
f[2] := 6 * y2 / y1;
h := 0.01;                                                           
for i from 1 to 2 do                                                                                         
for j from 1 to 3 do
y[i,j] :=y[i,j-1] + h * subs(y1 = y[i,j-1],y2=y[i+1,j-1],f[i]): 
od:
od;
print(y);


maple под рукой к сожалению нет, так что проверить немогу.
--------------------
Работая над решением задачи, всегда полезно знать ответ.
PM MAIL Skype   Вверх
aurealat
Дата 9.6.2009, 23:17 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 3
Регистрация: 25.4.2009

Репутация: нет
Всего: нет
Плиз спасите утопающих
Мучаем вложенные циклы в Maple но почему-то совсем ничего не работает 

Код

restart;   
with(plots):    
a:=4: b:=5: 
 for l from 1 to 10 do  for k from 0 to 15 do
m=(Pi*l/a)^2-(Pi*k/b)^2;
plot3d(cos((Pi*l/a)*x)*cos((Pi*k/b)*y),x = -a .. a, y = -b .. b, style = patch, scaling = constrained, axes = boxed);
od:
od;
PM MAIL   Вверх
cardinal
Дата 9.6.2009, 23:39 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Инженер
****


Профиль
Группа: Экс. модератор
Сообщений: 6003
Регистрация: 26.3.2002
Где: Германия

Репутация: 15
Всего: 99
Попробуй end do вместо od
http://books.google.ru/books?id=rt9xKPLCFc...lt&resnum=8


--------------------
Немецкая оппозиция потребовала упростить натурализацию иммигрантов
В моем блоге: Разные истории из жизни в Германии

"Познание бесконечности требует бесконечного времени, а потому работай не работай - все едино".  А. и Б. Стругацкие
PM   Вверх
aurealat
Дата 9.6.2009, 23:41 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 3
Регистрация: 25.4.2009

Репутация: нет
Всего: нет
пробовала при end do тоже не работает
PM MAIL   Вверх
cardinal
Дата 10.6.2009, 00:41 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Инженер
****


Профиль
Группа: Экс. модератор
Сообщений: 6003
Регистрация: 26.3.2002
Где: Германия

Репутация: 15
Всего: 99
Код

restart;   
with(plots):    
a:=4: b:=5: 
for l from 1 to 10 do
   for k from 0 to 15 do
      m=(Pi*l/a)^2-(Pi*k/b)^2;
      plot3d(cos((Pi*l/a)*x)*cos((Pi*k/b)*y),x = -a .. a, y = -b .. b, style = patch, scaling = constrained, axes = boxed)
   end do
end do;

Вот это один в один скопируй, может пойдет...


--------------------
Немецкая оппозиция потребовала упростить натурализацию иммигрантов
В моем блоге: Разные истории из жизни в Германии

"Познание бесконечности требует бесконечного времени, а потому работай не работай - все едино".  А. и Б. Стругацкие
PM   Вверх
aurealat
Дата 10.6.2009, 01:19 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 3
Регистрация: 25.4.2009

Репутация: нет
Всего: нет
Цитата(cardinal @ 10.6.2009,  00:41)
Вот это один в один скопируй, может пойдет...

неа, ничего не рисует, может это глюк самой программы?
PM MAIL   Вверх
cardinal
Дата 10.6.2009, 02:29 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Инженер
****


Профиль
Группа: Экс. модератор
Сообщений: 6003
Регистрация: 26.3.2002
Где: Германия

Репутация: 15
Всего: 99
Цитата(aurealat @  9.6.2009,  23:19 Найти цитируемый пост)
неа, ничего не рисует, может это глюк самой программы? 

Да не он все рисует, но т.к. паузы нет между командами, то увы ничего не видно smile

Можно так сделать (но это вариант для одного значения l)
Код

a:=4: b:=5: 
l:=1:
animate( plot3d, [cos((Pi*l/a)*x)*cos((Pi*k/b)*y),x = -a .. a, y = -b .. b], k=0..15, frames=100);

Покопай в этом направлении, может можно второй параметр в animate добавить...


--------------------
Немецкая оппозиция потребовала упростить натурализацию иммигрантов
В моем блоге: Разные истории из жизни в Германии

"Познание бесконечности требует бесконечного времени, а потому работай не работай - все едино".  А. и Б. Стругацкие
PM   Вверх
Khanty2012
Дата 23.3.2014, 18:55 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 1
Регистрация: 23.3.2014

Репутация: нет
Всего: нет
with(linalg); with(plots); 
a := [0, 0, 0]; b := [1, 1, 1]; 
N := 10; MMM := 5;
t := [[1/3, 1/2, 2/3, 3/4, 4/5], [1/3, 1/2, 2/3, 3/4, 4/5], [1/3, 1/2, 2/3, 3/4, 4/5]]; 

for i to nops(a) do

if i = 1 then reshen[x, i] := piecewise(x <= t[i, 1], x^2, x <= t[i, 2], (1/3)*x, t[i, 2] <= x and x <= t[i, 3], x^3+1/24, t[i, 3] <= x and x <= t[i, 4], .3379629630, t[i, 4] <= x and x <= t[i, 5], x^2-97/432, t[i, 5] <= x and x <= b[i], 4487/10800, 0); 
p[x, i] := piecewise(a[i] <= x and x <= b[i], 1, 0) fi; 

if i = 2 then reshen[x, i] := piecewise(a[i] <= x and x <= b[i], (4487/10800)*x^2, 0); p[x, i] := piecewise(a[i] <= x and x <= b[i], 1, 0); Q[x, i] := 1 fi;

 if i = 3 then reshen[x, i] := piecewise(a[i] <= x and x <= b[i], -(13461*(1/10800))*x^2+(4487/2700)*x, 0); p[x, i] := piecewise(a[i] <= x and x <= b[i], 1, 0); Q[x, i] := 1 fi;

 for r to MMM do fs[r, i] := limit(diff(reshen[x, i], x), x = t[i, r], left) end do: end do:

Maple неправильно считает предел справа!!! Где я ошибся?
PM MAIL   Вверх
  
 Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman
 
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »

Подписаться на тему | Подписка на этот форум | Скачать/Распечатать тему


 




[ Время генерации скрипта: 0.0788 ]   [ Использовано запросов: 21 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
 По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.