Список литературы

1. Г.Г. Рапаков, С.Ю. Ржеуцкая «Программирование на языке Pascal», СПб, «БХВ-Петербург», 2004 г.-467с.

2. В.М. Пестриков, А.Н. Маслобоев «Turbo PASCAL 7.0 изучаем на примерах» «Наука и техника», СПб, 2004 г.-356с.

3. Н.Культин «Turbo PASCAL в задачах и примерах», СПб, «БХВ-Петербург», 2003 г.-254с.

4. А. Адаменко «Turbo PASCAL на примерах из математики»,СПб, «БХВ-Петербург», 2005 г.-397с.

5. Г.Ч. Шушкевич, С.В. Шушкевич «Введение а MathCAD 2000». Учебное пособие./Гродно, ГрГУ,2001. – 140с.

Приложение 1

Программа решения уравнения методом касательных на языке Pascal.

program kasat; {уточнение корней уравнения методом касательных}

uses crt;

var

xk, xk1: real; {k-e и (k+1)-е приближение корня}

eps: real; {точность расчета}

m: integer; {кол-во корней}

n: integer; {счетчик итераций}

i: integer; {параметр цикла}

{подпрограмма вычисления значения функции f(x)}

function f (x : real) : real;

begin

f:=5*x*sqr(x)+4*sqr(x)+4*х+3;

end;

{подпрограмма вычисления значения производной функции }

function fi (x : real) : real;

begin

fi:=15*sqr(x)+8*x+4;

end;

{основная программа}

begin

clrScr;

write('Введите число корней = ');

readln(m);

write('Введите точность расчета  = ');

readln(eps);

For i:=1 to n do {Цикл по корням}

begin

writeln;

writeln(' Уточнение’, i , 'корня = ');

writeln ('----------------------');

write('Введите начальное приближение корня = ');

readln(xk1); {Ввод начального приближения}

n:=0; {Счетчик числа итераций}

repeat

xk:=xk1;

xk1:=xk - f(xk) / fi(xk);

n:=n+1;

until (abs(xk1-xk)<eps) or (n>100);

writeln('Число итераций n = ',n); {Вывод результатов расчета}

if n<101 then

begin

writeln('Корень уравнения = ', xk1:12:9);

writeln('Значение функции f(', xk1:12:9, ') = ', f(xk1):10);

end

else

writeln('Решение не найдено');

end;

readln;

end.

Приложение 2

Программа расчета определенного интеграла методом трапеций на языке Pascal.

program trapecia; {Вычисление интеграла по формуле трапеций}

uses crt;

var

h : real; {Элементарный отрезок интегрирования}

Sum : real; {Переменная для подсчета суммы}

a, b : real; {Пределы интегрирования}

integ : real; {Значение интеграла}

i : integer;

n : integer; {Число разбиений области интегрирования на

элементарные отрезки}

function f (x : real) : real; {Интегрируемая функция f(х)}

begin

f:=x*sqr(x)*sqr(sin(x));

end;

begin

clrScr;

{Ввод исходных данных}

write ('Нижний предел интегрирования  = ');

readln (a);

write ('Верхний предел интегрирования = ');

readln (b);

write ('Число разбиений отрезка интегрирования n = ');

readln (n);

{Расчет интеграла}

h:=(b-a)/n; {Длина элементарного отрезка}

Sum:=0; {Начальное значение суммы}

for i:=0 to n-1 do

Sum:=Sum+f(a+(i*h); {Подсчет суммы}

integ:=h*(f(a)/2+Sum+f(b)/2);

{Вывод результата расчета}

writeln('Интеграл от функции f = ', integ:12:9);

readln;

end.

	Страница 19