(Назад)
(Cкачать работу)графика функции. Модуль использует стандартные модули «crt» и «graph», а также библиотечный модуль inout, в котором содержатся описания массива для хранения значений функции и файла для хранения данных. unit count;graph,inout;m=array[1..1] of real;=^m;=array[1..2] of pm;schet(var f:fm; var y:real; x:real);result (var f,n:fm);grafik(var f,n:fm);
Процедура schet используется для вычисления значения у(х) для одного конкретного значения х, которое, наряду с данными для вычисления, передается этой процедуре. Вычисленное значение у(х) передается при помощи механизма передачи параметров. При вычислении используется полином Лагранжа, при этом возникает необходимость обращения по номерам к данным для вычисления. Поэтому все значения до начала вычислений заносятся в динамический массив. procedure schet;i,j:integer; m:z; c,u,d:real; w:mas;(f);(m[1],6*filesize(f)); getmem(m[2],6*filesize(f));i:= 1 to filesize(f) do begin(f,w);[1]^[i]:=w[1];[2]^[i]:=w[2];;:=0;i:= 1 to filesize(f) do begin-вычисление у(х)по формуле полинома Лагранжа:=1; d:=1; c:=1;j:=1 to filesize(f) doij then begin:=u*(x-m[1]^[j]);:=d*(m[1]^[i]-m[1]^[j]);;:=(u/d);:=y+c*m[2]^[i];;(m[1],6*filesize(f)); freemem(m[2],6*filesize(f));(f);
end;
Процедура result имеет своими параметрами оба файла - для хранения данных для вычисления и файл со значениями х, для которых и производится вычисление. Эта процедура считает значение у(х) и записывает его в соответствующий файл. Именно эта процедура используется в главной программе. Предыдущая является вспомогательной.
procedure result;i:integer; m:mas;(n);i:=1 to filesize(n) do begin(n,m);(f,m[2],m[1]);(n,filepos(n)-1);(n,m);;
close(n);
end;
Процедура grafik использует два файла с данными. При помощи значений, содержащихся в этих файлах, она строит график функции, оси координат в графическом режиме. Также отмечаются точки, заданные изначально и вычисленные в процессе программы, и некоторые значения по осям.
программа функция полином лагранж
procedure grafik;m:z; i,j,k,gd,gm,cx,cy,xo,yo:integer; koefy,koefx:integer;,buf2,q,p:real; w,e:mas; s:string[7];(f); reset(n);:=filesize(f)+filesize(n);(m[1],6*k); getmem(m[2],6*k);i:= 1 to filesize(f) do begin(f,w);[1]^[i]:=w[1];[2]^[i]:=w[2];;i:= filesize(f)+1 to k do begin(n,w);[1]^[i]:=w[1];[2]^[i]:=w[2];;i:=1 to k-1 doj:=i+1 to k doabs(m[1]^[i])>abs(m[1]^[j]) then begin:=m[1]^[i];:=m[2]^[i];[1]^[i]:=m[1]^[j];[2]^[i]:=m[2]^[j];[1]^[j]:=buf1;[2]^[j]:=buf2;;:=m[2]^[1]; p:=m[2]^[1]; xo:=0; yo:=0;i:=1 to k do begin -
определение количества отрицательных чисел для сохранения правильности расстановки осей if m[1]^[i]
Похожие работы
Интересная статья: Быстрое написание курсовой работы