Компьютерные технологии в металлургии и литейном производстве. Дембовский В.В. - 28 стр.

UptoLike

Составители:

Рубрика:

Это можно сделать по оценке исходного дифференциаль-

ного уравнения. Так, заметим, что в нашем примере после окон-

чания переходного процесса изменения y от его начального y0

до конечного y1 = const состояний имеем dy/dt = 0. Отсюда

f(t,y) = 0 и y1 := 1.

Определяем вектор начальных условий:

ic := y0

и функцию производной:

D(t,Y) := f(t,Y)

Здесь второй аргумент Y является вектором вычисляемых

значений функции y(t).

Формируем матрицу решений

S:=rkfixed(ic, t0, t1, n, D).

Определяем её столбцы:

первый (фактически нулевой)

t:=S

<0>

и второй

y:=S

<1>

Если нас интересуют численные значения аргумента t и

функции y в отдельных точках решения, то после последнего

выражения можно ввести

S =

Это вызовет появление на экране начальной области мат-

рицы, содержащей в первом столбце значения

t, а во втором –

y. Её можно расширить путем действия [1Л], после чего растя-

нуть вниз по экрану с помощью мыши.