программа ПЗЗ. Решение системы из двух дифференциальных: уравнений методом Рунге - Кутта. Ввод: .vo = R00, уо = RQ\ = = RQ2, Zo = /?03 = /?04, А/2 = /?07, данные у{х, у, г) и г(х,у,г) при x-RQQ, y-RQ2 и z->-ROA.

и Zo = 0,49235 получим при h = 0,1:

Программа П34

уравнений методом =. /?02, Zo = R03 = у(х, у, Z, w), z(x, z-ROA и w-R06.

SBR SBR STO

x:t STO 08

10 STO 07

11 SUM 00 RCL 05 = Ч- RCL 03 SUM 07 + SUM 09 SUM STO 06 x:.ti

. Решение системы нз трех дифференциальных Рунге - Кутта. Ввод; х = R00, г/о = WI = = R04, Wo = R05 .= RQ6, h/2 = RU, данные у, г, w) и ui{x, у, Z, w) при x~RQO, y-RQi2, Вывод результатов: х, w, z м у ъ регистр X.

09 + +

SBR STO

RCL RCL SBR

STO 04 01 =

X 2

03 =

01 =

SUM 09

= STO 06

STO 04 RCL

STO 02 RCL

SUM 09 -h

SUM 08 SUM 08

SUM 08

RCL = STO 2 +

SUM 08

RCL X

SUM 07

SBR SBR

03 =

2 4-

Для системы dx

dz dx

= - 6г - \ \y -&w;

y

при xo = 0, «0 = -II, 0 = 3, wo = 0 и h

0,1 получим:

Точное решение w {x) = 2e- - e- - e-.

Время счета по программам ПЗЗ и П34 составляет около 20 и 25 с на каждый шаг. Программы П32-34 наглядно иллюстрируют возможности программируемых микрокалькуляторов с повышенной емкостью запоминающих устройств.

Приложение 3

ПАКЕТ ПРОГРАММ С ОСОБЫМИ ПРИЕМАМИ ПРОГРАММИРОВАНИЯ ДЛЯ МИКРОКАЛЬКУЛЯТОРОВ ЭЛЕКТРОНИКА БЗ-34

Программа П3.1. Вычисление факториала я1 прн п = РХ, включая случай я = О (см. программу 3.16).

ВП ПО 1 ИПО X FLO 03 С/П

Программа П3.2. Вычисление значений многочлена степени п 12 по схеме Горнера. Ввод: оо = РД, а%.±а„ = Р1-j-Рл и да-

лее. для каждого х п - РО и х = РХ (см. программу 4.10).

t t С;с КИП + X FLO 03 ИПД + С/П БП 00

Программа ПЗ.З. Решение системы из п < 13 нелинейных уравнений Xi = Fi{Xi} модифицированным методом Зейделя [29]. Подпрограмма должна заносить лг1-*-РД, л;2-+-РС, Xin-t-Pl, Ввод: начальное значение х, В/О С/П, начальное значение Xi С/П и т. д., затем число итераций от = РХ БП 09 С/П получаем г,) С/П (получаем Хг) и Т. Д.

ПД 1 3 по XY С/П КПО БП 05 ПО ПП 21 FLO 10 1 4 ПО КИПО С/П БП 17 ........................ В/О

Для системы Xj = \/lXi (xi + b) - 1]/2; Х2 = л/хГ+ПтГхГ подпрограмма имеет вид

ИПС 5 + ИПД X I - 2 -v- FV ПД t Fin 3 X + fV~ ПС

При начальных х, - 4, хг =2,5 и m = 6 получим xi = 3,7607299 (точное значение 3,7568) и 2 = 2,7811093 (точное значение 2,7798) при tc « 70 с.

Программа П3.4. Поиск экстремума функции F{x) одной переменной X методом ее квадратичной аппроксимации в интервале Xi±:h. Ввод: данные F(x) при х заносимом в начало подпрограммы, h = РВ, г = PC, начальное значение х х, = РХ. Результат: Хм-»-РХ, F{Xh)PI- Заняты регистры О, 1, А, В и С. Программа вычисляет Xi = Хк для аппроксимирующей F(x) параболы до тех пор, пока (Xi+i - Xi) не станет меньше el

ПА ИПА ИПВ - ПП 58 ПО ИПА ПП 58

П1 ИПА ИПВ + ПП 58 t ИПА 2 X

ИПВ - X ИПА 2 X ИПВ + ИПО X

-f ИП1 4 X ИПА X - XY ИПО -Ь

ИП1 2 X - 2-4- ИПА XY ПА

- Fx ИПС - Fx < О 01 ИПА С/П ... В/О

Для F{x) = \0х - 2x2 + 0,1л: при х, = 3; А = 0,1; = МО-» получим хк = 3,33375 и F(x>) = 14,814815 при tc « 90 с.

Программа П3.5. Минимизация гладких функций F{Xi) с и <9 переменными методом координатного спуска, Ввод: данные F{Xi)