Результати співпадають з одержаними у розділі 3 в межах точності.

Висновки

Текст модуля бібліотеки програм

UNIT biblunit; {Бiблiотека пiдпрограм}

INTERFACE

uses crt;

const nom=50;

TYPE ar=array[1..50] of integer;

chyslo=record

a:array[1..50] of integer; {знаки до коми}

b:array[1..50] of integer; {знаки пiсля коми}

n:integer; {число знакiв до коми}

m:integer; {число знакiв пiсля коми}

end;

var zero,od:chyslo;

fi:text;

function max(n,m:integer):integer;

procedure vvid(var x:chyslo);

procedure vyvid(x:chyslo);

function poriadok(x:integer):integer;

function por(x:chyslo):integer;

function mpor(x:chyslo):integer;

function riv(x,y:chyslo):boolean;

function stdes(x:chyslo):integer;

function bilse(x,y:chyslo):boolean;

procedure dodavan(x,y:chyslo;var z:chyslo);

procedure vidnim(x,y:chyslo;var z:chyslo);

procedure pidvst(l:integer;var x:chyslo);

procedure ponst(l:integer;var x:chyslo);

procedure mnoint(x:chyslo;l:integer;var y:chyslo);

procedure zaokrug(l:integer;var x:chyslo);

procedure mnojen(x,y:chyslo;var z:chyslo);

procedure stepin(x:chyslo;l:integer;var y:chyslo);

procedure dilen(x,y:chyslo;var z:chyslo);

procedure dilenna(x,y:chyslo;var z:chyslo);

procedure intchyslo(k:integer;var t:chyslo);

procedure zerod;

implementation

function max(n,m:integer):integer;

begin

{бiльше з 2 чисел}

if (n>=m) then max:=n else max:=m

end;

procedure vvid(var x:chyslo);

{ввiд числа}

var i:integer;

begin

clrscr;

x:=zero;

with x do begin

write('Введiть число знакiв до коми ');

readln(n);

write('Введiть число знакiв пiсля коми ');

readln(m);

for i:=n downto 1 do begin gotoxy(20-i,10);read(a[i]);end;

gotoxy(20,10);

write('.');

for i:=1 to m do begin gotoxy(20+i,10);read(b[i]);end;

end;

end;

procedure vyvid(x:chyslo);

var i:integer;

begin

{вивiд числа}

with x do begin

if (n>0) then

for i:=n downto 1 do write(fi,a[i])

else

write(fi,'0');

write(fi,'.');

if (m>0) then

for i:=1 to m do write(fi,b[i]);

end;

end;

function poriadok(x:integer):integer;

var i,des:integer;

begin

{перша значуща цифра}

i:=1;des:=1;

while x>=des do begin

des:=des*10;

i:=i+1;

end;

poriadok:=i-1;

end;

function por(x:chyslo):integer;

var i:integer;

begin

{перша значуща цифра}

i:=nom;

while ((x.a[i]=0)and(i>0)) do i:=i-1;

por:=i;

end;

function mpor(x:chyslo):integer;

var i:integer;

begin

{число знакiв пiсля коми}

i:=nom;

while ((x.b[i]=0) and (i>0)) do

begin

i:=i-1;

end;

mpor:=i;

end;

function riv(x,y:chyslo):boolean;

var i,ni:integer;

begin

ni:=0;

x.n:=por(x);

x.m:=mpor(x);

y.n:=por(y);

y.m:=mpor(y);

if((x.n<>y.n)or(x.m<>y.m)) then ni:=1

else

begin

for i:=1 to nom do if (x.a[i]<>y.a[i]) then ni:=1;

for i:=1 to nom do if (x.b[i]<>y.b[i]) then ni:=1

end;

if (ni=1) then riv:=false else riv:=true;

end;

function stdes(x:chyslo):integer;

{розряд числа}

var i,n1:integer;

begin

x.n:=por(x);

x.m:=mpor(x);

if(riv(x,zero)) then n1:=-nom-1

else

if (por(x)>0) then n1:=por(x)-1

else

begin

i:=1;

while (x.b[i]=0) do i:=i+1;

n1:=-i;

end;

stdes:=n1;

end;

function bilse(x,y:chyslo):boolean;

var ni:boolean;

i,l:integer;

begin

{перевiрка яке число бiльше}

if (riv(x,y)) then ni:= false

else

i:=nom;

while ((x.a[i]=y.a[i])and(i>0)) do i:=i-1;

if (i>0) then

if(x.a[i]>y.a[i]) then ni:=true else ni:=false;

if (i=0) then

begin

i:=1;

while ((x.b[i]=y.b[i])and(i<nom)) do i:=i+1;

if(x.b[i]>y.b[i]) then ni:=true else ni:=false;

end;

bilse:=ni;

end;

procedure intchyslo(k:integer;var t:chyslo);

begin

t:=od;

if (k<10) then

begin

t.n:=1;

t.a[1]:=k;

end;

if ((k>=10) and (k<100)) then

begin

t.n:=2;

t.a[2]:=k div 10;

t.a[1]:=k-t.a[2]*10;

end;

end;

procedure dodavan(x,y:chyslo;var z:chyslo);

var i,perenos:integer;

begin

{додавання двох чисел}

z:=zero;

z.m:=max(x.m,y.m);

z.n:=max(x.n,y.n);

perenos:=0;

for i:=z.m downto 1 do begin

z.b[i]:=x.b[i]+y.b[i]+perenos;

if z.b[i]>=10 then

begin

z.b[i]:=z.b[i]-10;

perenos:=1;

end

else perenos:=0;

end;

for i:=1 to z.n do begin

z.a[i]:=x.a[i]+y.a[i]+perenos;

if z.a[i]>=10 then

begin

z.a[i]:=z.a[i]-10;

perenos:=1;

end

else perenos:=0;

end;

if (perenos=1) then begin z.n:=z.n+1;z.a[z.n]:=1 end;

z.m:=mpor(z);

end;

procedure vidnim(x,y:chyslo;var z:chyslo);

var i,j:integer;

begin

{вiднiмання двох чисел}

z:=zero;

z.m:=max(x.m,y.m);

z.n:=max(x.n,y.n);

for i:=z.m downto 1 do

if (x.b[i]>=y.b[i]) then z.b[i]:=x.b[i]-y.b[i]

else

begin

x.b[i]:=x.b[i]+10;

z.b[i]:=x.b[i]-y.b[i];

j:=i-1;

while ((x.b[j]=0)and(j>=1)) do

begin {4}

x.b[j]:=x.b[j]+9;

j:=j-1;

end;

if (j>=1) then x.b[j]:=x.b[j]-1

else

begin

j:=1;

while(x.a[j]=0) do

begin

x.a[j]:=9;

j:=j+1;

end;

x.a[j]:=x.a[j]-1;

end;

end;

i:=1;

for i:=1 to z.n do

if(x.a[i]>=y.a[i]) then z.a[i]:=x.a[i]-y.a[i]

else

begin

x.a[i]:=x.a[i]+10;

z.a[i]:=x.a[i]-y.a[i];

j:=i+1;

while (x.a[j]=0) do

begin

x.a[j]:=x.a[j]+9;

j:=j+1;

end;

x.a[j]:=x.a[j]-1;

end;

z.n:=por(z);

z.m:=mpor(z);

end;

procedure pidvst(l:integer;var x:chyslo);

procedure pidv(var x:chyslo);

var i:integer;

begin

with x do begin

for i:=n+1 downto 2 do a[i]:=a[i-1];

a[1]:=b[1];

for i:=1 to m-1 do b[i]:=b[i+1];

b[m]:=0;

x.n:=por(x);

x.m:=mpor(x);

end;

end;

var i:integer;

begin

{домноження на 10^l}

for i:=1 to l do pidv(x);

x.n:=por(x);

x.m:=mpor(x);

end;

procedure ponst(l:integer;var x:chyslo);

procedure pon(var x:chyslo);

var i:integer;

begin

with x do begin

for i:=m+1 downto 2 do b[i]:=b[i-1];

b[1]:=a[1];

for i:=1 to n-1 do a[i]:=a[i+1];

a[n]:=0;

x.n:=por(x);

x.m:=mpor(x);

end;

end;

var i:integer;

begin

{дiлення на 10^l}

for i:=1 to l do pon(x);

x.n:=por(x);

x.m:=mpor(x);

end;

procedure mnoint(x:chyslo;l:integer;var y:chyslo);

var i,perenos,dob:integer;

begin

{домноження на цiле число}

y:=zero;

perenos:=0;

for i:=x.m downto 1 do

begin

dob:=x.b[i]*l+perenos;

if (dob>=10) then begin

perenos:=dob div 10;

dob:=dob-10*perenos;

end

else perenos:=0;

y.b[i]:=dob;

end;

for i:=1 to x.n do

begin

dob:=x.a[i]*l+perenos;

if (dob>=10) then begin

perenos:=dob div 10;

dob:=dob-10*perenos;

end

else perenos:=0;

y.a[i]:=dob;

end;

if (perenos>0) then begin

y.a[x.n+1]:=perenos;

end;

y.n:=por(y);

y.m:=mpor(y);

end;

procedure zaokrug(l:integer;var x:chyslo);

{заокруглення з точнiстю}

var i:integer;odst:chyslo;

begin

odst:=od;

ponst(l,odst);

if x.b[l+1]<5 then for i:=l+1 to nom do x.b[i]:=0

else

begin

for i:=l+1 to nom do x.b[i]:=0;

dodavan(x,odst,x);

end;

end;

procedure mnojen(x,y:chyslo;var z:chyslo);

var i:integer;

s,t:chyslo;

begin

{множення двох чисел}

s:=zero;

for i:=y.m downto 1 do

begin

mnoint(x,y.b[i],t);

ponst(i,t);

dodavan(s,t,s);

end;

for i:=1 to y.n do

begin

mnoint(x,y.a[i],t);

pidvst(i-1,t);

dodavan(s,t,s);

end;

z:=s;

zaokrug(nom-1,z);

z.n:=por(z);

z.m:=mpor(z);

zaokrug(z.m,z);

z.n:=por(z);

z.m:=mpor(z);

end;

procedure stepin(x:chyslo;l:integer;var y:chyslo);

var i:integer;

begin

y:=od;

for i:=1 to l do mnojen(y,x,y);

end;

procedure dilen(x,y:chyslo;var z:chyslo);

label 1;

var u,v:chyslo;i,j,k,l,t1,t2:integer;

begin

{дiлення двох чисел}

t1:=0;t2:=0;

z:=zero;

if(riv(x,y)) then z:=od;

if (bilse(x,y)) then

begin

i:=0;

while (bilse(x,y)) do begin pidvst(1,y);i:=i+1 end;

t1:=i;

end;

if (bilse(y,x)) then

begin

i:=0;

while (bilse(y,x)) do begin pidvst(1,x);i:=i+1 end;

ponst(1,x);t2:=i-1;

end;

for i:=1 to nom-t2 do

begin

j:=0;

pidvst(1,x);

u:=zero;

1:mnoint(y,j,u);

if (riv(x,u)) then begin k:=j end

else

if (bilse(u,x)) then begin k:=j-1;end

else

begin j:=j+1;goto 1; end;

z.b[i]:=k;

mnoint(y,k,v);

vidnim(x,v,x);

end;

pidvst(t1,z);

ponst(t2,z);

zaokrug(nom-1,z);

z.n:=por(z);

z.m:=mpor(z);

end;

procedure dilenna(x,y:chyslo;var z:chyslo);

var i:integer;t1,t2,t3,t4,t5:chyslo;

{уточнення дiлення}

begin

t1:=x;t2:=y;

dilen(x,y,z);

for i:=1 to 20 do

begin

mnojen(y,z,t3);

if (bilse(t3,x)) then

begin

vidnim(t3,x,t4);

dilen(t4,y,t5);

vidnim(z,t5,z)

end

else

if (bilse(x,t3)) then

begin

vidnim(x,t3,t4);

dilen(t4,y,t5);

dodavan(z,t5,z);

end;

zaokrug(nom-1,z);

z.n:=por(z);

z.m:=mpor(z);

zaokrug(mpor(z),z);

z.n:=por(z);

z.m:=mpor(z);

end;

end;

procedure zerod;

var i:integer;

begin

with zero do

begin

n:=0;

m:=20;

for i:=1 to 50 do a[i]:=0;

for i:=1 to 50 do b[i]:=0;

end;

od:=zero;

od.n:=1;

od.a[1]:=1;

end;

end.

Текст програми kurs0.pas

program kurs0;

{$M 65520,0,655360}

uses biblunit,Crt;

{iлюстрацiя роботи програми}

var t1,t2,t3,t4,t5,t6:chyslo;

i,m:integer;

fi:text;

begin

zerod;

assign(fi,'tor\\res.txt');

rewrite(fi);

vvid(t1);

vvid(t2);

clrscr;

writeln(fi,'аргументи');

vyvid(t1);

writeln(fi);

vyvid(t2);

writeln(fi);

writeln(fi,'Сума');

dodavan(t1,t2,t3);

vyvid(t3);

writeln(fi);

writeln('Рiзниця');

vidnim(t1,t2,t4);

vyvid(t4);

writeln(fi);

writeln(fi,'Добуток');

mnojen(t1,t2,t5);

vyvid(t5);

writeln(fi);

writeln(fi,'Частка');

dilenna(t1,t2,t6);

vyvid(t6);

writeln(fi);

close(fi);

end.

Текст програми kurs1.pas

program kurs1;

{$M 65520,0,655360}

{знаходження числа e}

uses biblunit,Crt;

var factor,od1,suma,t1,t2,t3:chyslo;

i,m:integer;

begin

assign(fi,'tor\\rese.txt');

rewrite(fi);

clrscr;

zerod;

od1:=od;

ponst(47,od1);

factor:=od;

mnoint(od,2,suma);i:=2;

repeat

writeln('i=',i);

intchyslo(i,t1);

mnojen(factor,t1,factor);

dilen(od,factor,t2);

dodavan(suma,t2,suma);

i:=i+1;

until(bilse(od1,t2));

writeln(fi,'Значення числа e');

vyvid(suma);

close(fi);

end.

Текст програми kurs2.pas

program kurs2;

{$M 65520,0,655360}

{знаходження числа pi}

uses biblunit,Crt;

var ko1,ko2,ko3,suma,suma1,suma2,spi,t0,t1,t2,t3,od1,dwa,try,kt:chyslo;

i,m,k:integer;

begin

assign(fi,'tor\\respi.txt');

rewrite(fi);

clrscr;

zerod;

{блок знаходження кореня з 3}

od1:=od;

ponst(45,od1);

try:=od;

try.a[1]:=3;

dwa:=od;

dwa.m:=1;

dwa.a[1]:=0;

dwa.b[1]:=5;

ko1:=od;

ko2:=od;

ko1.m:=1;

ko2.m:=1;

ko1.b[1]:=7;

ko2.b[1]:=8;

repeat

dodavan(ko1,ko2,ko3);

mnojen(ko3,dwa,ko3);

mnojen(ko3,ko3,t2);

if (bilse(t2,try)) then ko2:=ko3 else ko1:=ko3;

if (bilse(t2,try)) then vidnim(t2,try,t1)

else vidnim(try,t2,t1);

until (bilse(od1,t1));

dilen(od,ko3,ko3);

writeln('корiнь пiдрахували');

{блок пiдсумовування ряду}

suma:=od;

m:=1;

t3:=od;

i:=1;

repeat

m:=-m;

mnoint(t3,3,t3);

k:=2*i+1;

intchyslo(k,kt);

mnojen(t3,kt,t1);

dilen(od,t1,t2);

if (m>0) then dodavan(suma,t2,suma)

else vidnim(suma,t2,suma);

i:=i+1;

until (bilse(od1,t2));

mnoint(suma,6,spi);

mnojen(spi,ko3,spi);

writeln(fi,'Значення число П');

writeln;

vyvid(spi);

writeln(fi);

close(fi);

end.

Текст програми kurs3.pas

program kurs3;

{$M 65520,0,655360}

{знаходження числа pi методом Ейлера}

uses biblunit,Crt;

var suma,suma1,suma2,spi,t0,t1,t2,t3,od1,dwa,try,kt,chot,dev:chyslo;

i,m,k:integer;

begin

assign(fi,'tor\\respi1.txt');

rewrite(fi);

clrscr;

zerod;

od1:=od;

ponst(45,od1);

dwa:=od;

dwa.m:=1;

dwa.a[1]:=0;

dwa.b[1]:=5;

mnojen(dwa,dwa,chot);

try:=od;

try.a[1]:=3;

dilen(od,try,try);

mnojen(try,try,dev);

{dwa=1/2 try=1/3 chot=1/4 dev=1/9}

{блок пiдсумовування ряду}

suma1:=dwa;

m:=1;

t3:=od;

i:=1;

t0:=dwa;

repeat

m:=-m;

mnojen(t0,chot,t0);

k:=2*i+1;

intchyslo(k,kt);

dilen(od,kt,t2);

mnojen(t2,t0,t1);

if (m>0) then dodavan(suma1,t1,suma1)

else vidnim(suma1,t1,suma1);

i:=i+1;

until (bilse(od1,t1));

writeln('Перший крок пройшов');

suma2:=try;

m:=1;

t0:=try;

i:=1;

repeat

m:=-m;

mnojen(t0,dev,t0);

k:=2*i+1;

intchyslo(k,kt);

dilen(od,kt,t2);

mnojen(t2,t0,t1);

if (m>0) then dodavan(suma2,t1,suma2)

else vidnim(suma2,t1,suma2);

i:=i+1;

until (bilse(od1,t1));

dodavan(suma1,suma2,suma);

mnoint(suma,4,spi);

zaokrug(45,spi);

writeln(fi,'Значення число П');

writeln;

vyvid(spi);

writeln(fi);

close(fi);

end.

Результати роботи програми

Значення числа e

2.7182818284590452353602874713526624977572470936998

Значення число П (метод Шарпа)

3.14159 26535 89793 23846 26433 83279 50288 41971 69399

Значення число П (метод Ейлера)

3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510