Подпрограммы некоторых численных методов

12 13 14 15 16 17 18

100 REM ПОДПРОГРАММА РЕШЕНИЯ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИИ

110 REM МЕТОД ЖОРДАНА—ГАУССА

120 REM A(N, N+1) МАТРИЦА КОЭФФИЦИЕНТОВ, X(N) МАТРИЦА НЕИЗВЕСТНЫХ

130 REM N — ЧИСЛО УРАВНЕНИИ

140 REM СТОЛБЕЦ N+1 — ПРАВЫЕ ЧАСТИ УРАВНЕНИЙ

150 FOR I=1 ТО N

160 FOR J=l TO N

170 IF J=l THEN 190

180 M=A (J, I) /A (I, I)

190 FOR K=1 TO N+1

200 A(J, K)=A(J, K)—M*A (I, К)

210 NEXT К

220 NEXT J

230 NEXT I

240 FOR I=1 TO N

250 X(I)=A(1, N+1)/A(I, I)

260 NEXT I

270 RETURN

100 REM ПОДПРОГРАММА РЕШЕНИЯ ПО МЕТОДУ ЭЙЛЕРА

110 REM ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ

120 REM ШАГ РАСЧЕТА — Н

130 REM УРАВНЕНИЕ ЗАДАЕТСЯ ОПЕРАТОРОМ DEF FND...

140 REM ПОДПРОГРАММА ВЫЗЫВАЕТСЯ НА КАЖДОМ ШАГЕ

150 F=FND (X, Y)

160 Y=Y+H*F

170 Х=Х+Н

180 RETURN

100 REM ПОДПРОГРАММА ВЫЧИСЛЕНИЯ ПО МЕТОДУ РУНГЕ-КУТТА

110 RЕМ С ПОСТОЯННЫМ ШАГОМ

120 REM ШАГ — Н, НАЧ. УСЛОВИЯ ДЛЯ X, Y — R, Q

130 REM МАКСИМАЛЬНОЕ ЗНАЧЕНИЕ Х—М

140 REM РЕШАЕМОЕ УРАВНЕНИЕ СЛЕДУЕТ ЗАДАТЬ С ПОМОЩЬЮ

150 REM ОПЕРАТОРА DEF FND (X, Y)=....

160 L=INT ((M—R)/H)

170 DIM X(L+1), Y(L+1)

180 I=0

190 X(0)=R

200 Y(0)=Q

210 PRINT “ШАГ X Y=F(X)”

220 PRINT TAB (2); 1; TAB (10);

230 PRINT USING “*****,*****”:X(I);

240 PRINT TAB (23)

250 PRINT USING “*****,*****”: Y(l)

260 X(I+1)=X(1)+H

270 IF x(i+i)>|m THEN 350

280 K=FND (X(l), Y(I))

290 Kl=FND(X(I)+0,5*H,Y(I)+0,5*H*K)

300 K2= FND(X(I)+0,5*H,Y(I)+0,5*H*K1)

310 K3=FND(X(1)+H, Y(1)+H*K2)

320 Y(I+1)=Y(I)+(H/6)*(K+2*K1+2*K2+K3)

330 I = I + 1

340 GOTO 220

350 RETURN

СОДЕРЖАНИЕ

Раздел 1.	Общие методические рекомендации по изучению дисциплины ……………………………………………………..
Раздел 2.	Методические указания по изучению содержания тем и разделов курса ………………………………………………
Раздел 3.	Задания для контрольных работ и указания по их выполнению …………………………………………………….
Список использованных источников………………………..
Приложения…………………………………………………………………
Содержание ……………………………………………………………...