ПРИЛОЖЕНИЕ 1. Исходный текст главной программы

Исходный текст главной программы

use msflib! подключение модуля стандартных процедур

implicit none! указание обязательного объявления переменых

! -----------------------------------------------------------

integer, parameter:: fi=1! устройство ввода

integer, parameter:: fpw=2! устройство вывода

integer, parameter:: MAX=100! максимальное количество точек

integer, parameter:: NT=150! количество точек кривой полинома

real(8), parameter:: D_X=0.05! припуск на диапазон рисования

real(8), parameter:: D_Y=0.05! припуск на диапазон

real(8) XN, XK, YN, YK, HX! границы и шаг изменения Х графика

integer i! счетчик точек полинома

real(8) Xmin, Xmax! минимальный и максимальный элементы X

real(8) Ymin, Ymax! минимальный и максимальный элементы Y

real(8):: X(0: MAX), Y(0: MAX)! таблица данных X(N), Y(N)

real(8):: C(0: MAX)! массив коэффициентов аппрокс.полинома Pm(x)

real(8):: A(0: MAX, 0: MAX)! матрица коэффициентов системы уравнений

integer N! количество точек табл. функции (i=0..N-1)

! -----------------------------------------------------------

logical res! логический результат графических функций

integer(2) XE, YE! XE, YE - размеры экрана в пикселях

INTEGER(2) status! целый результат графических функций

logical(2) finv /.true./! Ось у направлена снизу вверх

type(windowconfig) wc! структура параметров окна

! -----------------------------------------------------------

open(unit=fi, file='DZ_V6.txt')! открытие файла данных

read(fi, *) n! чтение количества узлов аппроксимации

do i=0, n-1

read(fi, *) X(i), Y(i)! чтение координат узлов

Enddo

close(fi)! закрытие файла

! -----------------------------------------------------------

! Автоматическая настройка конфигурации окна

data wc.numxpixels, wc.numypixels, wc.numtextcols, &

wc.numtextrows, wc.numcolors, wc.fontsize / 6*-1 /

wc.title = " Интерполяционный многочлен Лагранжа" C

res = getwindowconfig(wc)! Читаем параметры видеоокна

XE = wc.numxpixels! numxpixels - число пикселей по оси х

YE = wc.numypixels! numypixels - число пикселей по оси у

status = RECTANGLE($GFILLINTERIOR, INT2(0), INT2(0), XE, YE)

status = setcolor(1_2)! цвет рамки

status = RECTANGLE($GBORDER, XE/4_2-20, YE/4_2-20, &

3_2*XE/4_2+20, 3_2*YE/4_2+20)

call axis()! Рисуем оси координат

! Задание видового порта размером XE/2 * YE/2 в центре видеоокна

call setviewport(XE/4_2, YE/4_2, 3_2*XE/4_2, 3_2*YE/4_2)

! Оконная система координат (ОСК)

! -----------------------------------------------------------

Xmin=MinVal(X); Xmax=MaxVal(X)! поиск минимума и максимума Х

Ymin=MinVal(Y); Ymax=MaxVal(Y)! поиск минимума и максимума Y

XN=Xmin-D_X! левая граница Х графика

XK=Xmax+D_X! правая граница Х графика

YN=Ymin-D_Y! левая граница Y графика

YK=Ymax+D_Y! правая граница Y графика

status = setwindow(finv, XN, YN, XK, YK)

call curve()! Рисуем график

! -----------------------------------------------------------

ПРИЛОЖЕНИЕ 2

Секция CONTAINS с определением внутренних подпрограмм

! ******************************************************************

subroutine axis()! Рисуем оси координат

Type(xycoord) xy

status = setcolor(2_2)! цвет осей координат

call moveto(int2(XE/4 - 10), int2(YE/2), xy)

status = lineto(3_2*XE/4_2 + 10_2, YE/2_2)! Ось х

call moveto(int2(XE/2), int2(YE/4 - 10), xy)

status = lineto(XE/2_2, 3_2*YE/4_2 + 10_2)! Ось y