Как построить в Maple график решения системы ур-ий движения?

[Divergence]

Интересует построение графиков, используя Maple для функции $u_n=u_n(t)$ при фиксированном $n$, то есть по оси OX: $x=t$ по оси OY: $y=u_n$.
При этом функции $u_n(t)$ описываются линейными уравнениями

$$\frac{d^2 u_n(t)}{dt^2} = \sum^N_{m=-N} c(m) u_{n-m}(t) ,$$

где $n \in \mathbb{Z}$, а $c(m)$ - некоторая элементарная функция (например, $c(m)=m^{-2}$ для ненулевых $m$).

К сожалению, не знаю, как моделировать даже в простейшем случае

$$ \frac{d^2 u_n(t)}{dt^2} = u_{n+1}(t)- 2u_n(t)+ u_{n-1}(t) . $$

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

[GAA]

Тема перемещена из форума «Околонаучный софт» в форум «Карантин»
по следующим причинам:

Условие, даже в простейшем случае, не полное.
Если это система уравнений, в которой $n$ изменяется от 1 до $n_{max}$, то какой вид имеют уравнения системы при $n=1$ и $n=n_{max}$?
Решение должно находиться символьно или численно? Есть ли начальные условия?
В общем, надо уточнить условие.


Исправьте все Ваши ошибки и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума.

-- Fri 31.07.2015 08:51:40 --

Если я угадываю задачу правильно, то график строиться стандартным образом. Вот пример построения графиков для оригинальной системы Вольтерра в приведенной форме (после замены переменных, позволяющей свести задачу к системе с одним параметром) при фиксированном значении параметра 1.
Maple 7 (но будет работать и в более поздних версиях)

Оригинальная система Вольтерра в приведенном виде. Значение параметра положено равным 1
> eqs:= diff(u(t), t) = u(t)-u(t)*v(t), diff(v(t), t)= u(t)*v(t)-v(t):
Задаём начальные условия 
> ics:= u(0)=0.5, v(0)=0.5:
Получаем модуль численного решения 
> sol:= dsolve({eqs, ics}, numeric, range=0..100, output=listprocedure):
Строим графики решений u(t), v(t)
> plot([rhs(sol[2]), rhs(sol[3])], 0..10, color=[red, blue]);