Функция решения системы уравнений:

{программу подготовил студент группы 106410}

{Мойсеня Е. }

Program var_18 (moisenja);

uses crt;

type massiv = array [1.3,1.2] of real;

var

i: integer; {счётчик}

vyb: byte; {Для диалога}

syst: massiv; {для cсистемы уравнений}

g: real;

p, eps, x0: real; {для нелинейного уравнения}

x, w: real; {коэффициенты временной функции}

fi: byte; {угл для коэф-а v}

t, t0, tkon, tkv: real; {для функции}

j: byte; {для вывода таблицы (коорд x) }

{*************** решение системы ***************}

function sist (A: massiv): real;

var i,j: byte; {счётчики}

r1, r2: real; {решения системы}

begin

for j: = 1 to 2 do

for i: = 3 downto 1 do

A [i,j]: = A [i,j] /A [1,j];

for i: = 1 to 3 do A [i,2]: = A [i,1] - A [i,2];

r2: = A [3,2] /A [2,2];

r1: = (A [3,1] - r2*A [2,1]);

if abs (r1) >abs (r2) then sist: = r1 else sist: = r2;

end;

{**************************************************}

{** * решение нелинейного уравнения **************}

function nelin (fx0, feps: real): real;

var x, xpred: real;

function f (r: real): real; {Заданная функция}

begin

f: = exp (r) + ln (r) - 10*r;

end;

function fpr (r: real): real; {Её производная}

begin

fpr: = exp (r) - 10 + 1/r;

end;

begin

x: = fx0;

repeat

begin

xpred: =x;

x: = xpred - f (xpred) /fpr (xpred);

end;

until (f (x) < feps) and (abs (x-xpred) < feps);

nelin: = x;

end;

{*************************************************}

{************* функция **************************}

function ft (xf,wf,pf,gf,tf: real): real;

begin

ft: =abs (pf + gf + tf* (wf + tf*xf));

end;

{*************************************************}

BEGIN

{********** Решаем систему уравнений ************}

clrscr;

Writeln ('Решаем систему уравнений + a1z + b1v = d1 ');

Writeln (' + a2z + b2v = d2 ');

Writeln ('Введите коэффициенты a1, b1, d1, a2,. ');

for j: = 1 to 2 do

begin

for i: = 1 to 3 do

begin

gotoXY (25 + i*6, 4 + j);

read (syst [i,j]);

end;

end;

g: = sist (syst);

Writeln ('Наибольшее по модулю решение системы g= ',g: 4: 2);

Write ('Для продолжения любую клавишу'); readkey;

{ readln; }

{****** Решаем нелинейное уравнение ***************}

clrscr;

Writeln ('Решаем нелинейное уравнение e^ (x) +ln (x) - 10x=0 методом Ньютона ');

Write ('Вводим начальное значение x0 = '); readln (x0);

Write ('Вводим точность расчёта E = '); readln (eps);

p: = nelin (x0,eps);

Writeln ('Корень уравнения e^ (x) +ln (x) - 10x=0 x= ',p: 6: 4);

Write ('Для продолжения нажмите любую клавишу'); readkey;

{**** подготовка к выводу таблицы **************}

clrscr;

Writeln ('Сейчас происходит подготовка к выводу значений ');

Writeln ('временной функции xt^2 + wt + p + g, где ');

Writeln (' p - корень уранения e^ (x) +ln (x) - 10x=0 ');

Writeln (' m - больлший по модулю корень системы. ');

Writeln ('Ввод недостающих коэффициентов: ');

Write ('Введите x; x = '); readln (x);

Write ('w = tg (fi); введите fi (в градусах); fi = '); readln (fi);

w: = (sin (fi)) / (cos (fi));

Writeln ('');

Writeln ('Ввод параметров функции: ');

Write ('Введите начальное время t0; t0 = '); readln (t0);

Write ('Введите конечное время tkon; tkon = '); readln (tkon);

Write ('Введите время квантования tkv; tkv = '); readln (tkv);

Writeln ('');

Writeln ('В каком времени рассчитывать функцию? ');

Writeln (' 1 - если в машинном ');

Writeln (' 2 - если в реальном ');

readln (vyb);

Write ('Для вывода результатов всё готово. Нажмите любую клавишу'); readkey;

{****** в ы в о д р е з у л ь т а т о в *************}

clrscr;

j: =0;

gotoXY (1+j, 1);

Writeln ('+-----------------------------+ ');

gotoXY (1+j,

2);

Writeln ('¦ Значение t ¦ Значение f (t) ¦ ');

gotoXY (1+j,

3);

Writeln ('+------------+----------------¦ ');

t: = t0; i: =1;

while t <= tkon do

begin

gotoXY (1+j, i+3); write ('¦');

gotoXY (14+j, i+3); write ('¦');

gotoXY (31+j, i+3); write ('¦');

gotoXY (4+j, i+3);

writeln (t: 5: 3);

gotoXY (18+j, i+3);

writeln (ft (x,w,p,g,t): 6: 4);

if vyb = 2 then delay (trunc (tkv*58000));

if i > 19 then

begin

gotoXY (j+1, i+4);

Writeln ('+-----------------------------+ ');

i: = - 2;

j: = j+40;

end;

t: = t + tkv; i: = i + 1;

end;

gotoXY (j+1, i+3);

Writeln ('+-----------------------------+ ');

readkey;