Контрольная работа: Вычисление значения функции y(x)

Задание 1.1

Вычислить значение функции y(x) на промежутке [1;2.8] с шагом h=0.3.

 

В Pascal:

 

program z1;

uses crt;

var Xn,Xk,X,Y,H,Z:REAL;

begin

clrscr;

write('‚Enter Xn,Xk,H=');

readln(Xn,Xk,H);

X:=Xn;

repeat

z:= 1-exp((0.5)*x);

if z=0 then writeln('NO') else

Y:=ln(6.8*sqrt(x))*(exp((1/3)*(ln(abs(Z))))*(abs(Z)/Z));

writeln('X=',X:6:1,' Y=',Y:8:3);

X:=X+H;

until X>=Xk+H/2;

readkey;

end.

 

Результаты вычислений:

 

В Excel:

 

 

В MathCAD:


В Delphi:

 

 

var

i:integer;

s,Xn,z,F,h,Xk:real;

procedure TForm1.BitBtn1Click(Sender: TObject);

begin

Xn:=strtoFloat(edit1.text);

h:=strtoFloat(edit2.text);

Xk:=strtoFloat(edit3.text);

repeat

z:= 1-exp((0.5)*Xn);

F:=ln(6.8*sqrt(Xn))*(exp((1/3)*(ln(abs(Z))))*(abs(Z)/Z));

Memo1.Lines.Add('x='+FormatFloat('0.0',xn)+' y='+FormatFloat('0.000',F));

Xn:=Xn+h;

until Xn>Xk;

end;

procedure TForm1.BitBtn2Click(Sender: TObject);

begin

form1.close;

end;


end.

 

Вывод: т.к. ответы в Pascal,Excel,MathCAD и Delphi сошлись – решение верно.

 

Блок-схема алгоритма решения задания №1.1

 


Задание 1.2

 

Вычислить значение функции y(x) с условием на промежутке [-2;3] с шагом h=0.5

В Pascal:

program z2;

uses crt;

label 20;

var x,y,Xn,Xk,h:real;

begin clrscr;

writeln('Please ENTER Xn,Xk,h=');

readln(Xn,Xk,h);

x:=Xn;

while x<=Xk+h/2 do

begin

if x<0 then y:=1/sin(x) else

if (x>=0) and (x<2) then y:=(2*x)/((abs(x-2)/(x-2))*(exp(1/3*(ln(abs(x-2)))))) else

if x>2 then y:=sqr(cos(x))/(x-2)

else

begin

writeln('NO answer');

goto 20;

end;

writeln('x=',x:3:1,' y=',y:6:3);

20: x:=x+h;

end;

readkey;

end.

 

Результаты вычислений:

компьютерный вычисление функция алгоритм


В Excel:

 

В MathCAD:

 


В Delphi:

 

var

x,Xn,Xk,h,y: Real;

begin

Xn:=strtoFloat(edit1.text);

Xk:=strtoFloat(edit2.text);

h:=strtoFloat(edit3.text);

begin

x:=Xn;

while x<=Xk+h/2 do

begin

If (x=2) then

begin

memo1.Lines.Add('No answer');

end

else

begin

if x<0 then y:=1/sin(x) else

if (x>=0) and (x<2) then y:=(2*x)/((abs(x-2)/(x-2))*(exp(1/3*(ln(abs(x-2)))))) else

if x>2 then y:=sqr(cos(x))/(x-2);

memo1.Lines.Add('y='+FormatFloat('0.000',y));

Series1.AddXY(x,y,'',clred);

end;

 x:=x+h;

end;

end;

end;

end.

 

Вывод: т.к. ответы в Pascal, Excel, MathCAD и Delphi сошлись – решение верно.

Блок-схема алгоритма решения задания: 1.2

Нахождение функции, заданной условием:


Вычислить функцию y(x) по заданному массиву x

  

В Pascal:

program zadanie3;

uses crt;

const n=4;

var x,y:array[1..n] of real;

i:integer;

z,j,d:real;

begin

clrscr;

for i:=1 to n do

begin

write('Enter x[i]=');

readln(x[i]);

end;

for i:=1 to n do

begin

j:=exp(x[i])-8.35*x[i]*x[i];

z:=exp(1/3*ln(abs(j)))*(abs(j)/j);

d:=exp(x[i]*ln(37))+ln(abs(sin(2*x[i])))*(abs(2*x[i])/2*x[i]);

y[i]:=z/d;

writeln('x[i]=',x[i]:5:1,' y[i]=',y[i]:5:3);

end;

readkey;

end.

Результаты вычислений:

В Excel:

 

 


В MathCAD:

 

 


В Delphi:

var x,y:array[0..3]of real;

i:integer;

z,j,d:real;

begin

randomize;

for i:=0 to 3 do

begin

x[i]:=strtofloat(stringGrid1.cells[i,0]);

j:=exp(x[i])-8.35*x[i]*x[i];

z:=exp(1/3*ln(abs(j)))*(abs(j)/j);

d:=exp(x[i]*ln(37))+ln(abs(sin(2*x[i])))*(abs(2*x[i])/2*x[i]);

y[i]:=z/d;

stringgrid1.Cells[i,1]:=FormatFloat('0.000',y[i]);

end;

end;

procedure TForm1.BitBtn2Click(Sender: TObject);

begin

form1.close;

end;

end.

 

Вывод: т.к. ответы в Pascal,Excel,MathCAD и Delphi сошлись – решение верно.

Блок-схема алгоритма решения задания №1.3

Нахождение функции заданном массиве:

Задание 1.4

 

Найти сумму, если задана функция y(x), задан массив а, даны Xn, n, h.


 

В Pascal:

Program Summa;

uses crt;

const n=6;

Var x,j,z,d,Xn,h,F,S:Real;

i:Integer;

a:array [1..n] of Real;

BEGIN

clrscr;

randomize;

Write('enter please Xn,h=');

ReadLN(Xn,h);

for i:=1 to n Do

begin

Write('enter please a[i]=');

ReadLN(a[i]);

end;

x:=Xn;

S:=0;

for i:=1 to n Do

begin

j:=(ln(abs(3+x*x*x))*(abs(3+x*x*x)/(3+x*x*x)))-(exp((cos(x)*cos(x))*ln(2)));

z:=exp(1/7*ln(abs(j)))*abs(j)/j;

d:=sin(2*x)*cos(3*x)+sqrt(1+exp((x)*ln(2)));

F:=z/d;

S:=S+F*a[i];

x:=x+h;

end;

WriteLN('S=',S:10:3);

END.

Результаты вычислений:

 

В Excel:

 

 


В MathCAD:

 

В Delphi:

var a:array[0..5]of real;

i:integer;

x,s,Xn,z,F,j,d,h:real;

begin

Xn:=strtoFloat(edit1.text);

h:=strtoFloat(edit2.text);

x:=Xn;

s:=0;

for i:=0 to 5 do

begin

j:=(ln(abs(3+x*x*x))*(abs(3+x*x*x)/(3+x*x*x)))-(exp((cos(x)*cos(x))*ln(2)));  j:=(ln(abs(3+x*x*x))*(abs(3+x*x*x)/(3+x*x*x)))-(exp((cos(x)*cos(x))*ln(2)));

z:=exp(1/7*ln(abs(j)))*abs(j)/j;

d:=sin(2*x)*cos(3*x)+sqrt(1+exp((x)*ln(2)));

F:=z/d;

a[i]:=strtofloat(stringGrid1.cells[i,0]);

s:=s+F*a[i];

x:=x+h;

end;

memo1.Lines.Add(' Ñóììà= '+FormatFloat('0.000',s));

end;

procedure TForm1.BitBtn2Click(Sender: TObject);

begin

form1.close;

end;

end.

 

Вывод: т.к. ответы в Pascal, Excel, MathCAD и Delphi сошлись – решение верно.


Блок-схема алгоритма решения задания №1.4

Нахождение функции, если дан x, h, n, задан массив:

Основы программирования на языке Паскаль
Как работать с книгой Внимательно прочитайте соответствующий раздел теории (одну главу), разберите все примеры, чтобы вам все было понятно, при этом ...
Begin s:=s+(abs(cos(x))+abs(cos(x+h)))/2*h;
Begin f1:=exp(x/10)+sqrt(x/(x+3)); end;
Раздел: Рефераты по информатике, программированию
Тип: учебное пособие
Разностные схемы для уравнения переноса на неравномерных сетках
Министерство образования и науки Российской Федерации Федеральное агентство по образованию Государственное образовательное учреждение высшего ...
begin ut:= A*exp(p+r);end;
x[i]:=i*h[i];t[j]:=j*tau[j];
Раздел: Рефераты по информатике, программированию
Тип: дипломная работа
Object Pascal
... Основными символами языка Object Pascal являются: - символы _ - 26 больших и 26 малых латинских букв A,B, .Y,Z, a,b, ., y,z - 10 арабских цифр 0, 1, 2 ...
Function Exp (X: Real):
Function Ln (X: Real):
Раздел: Рефераты по информатике, программированию
Тип: реферат
Система математических расчетов MATLAB
ГОСУДАРСТВЕННЫЙ ИНЖЕНЕРНЫЙ УНИВЕРСИТЕТ АРМЕНИИ MATLAB УЧЕБНОЕ ПОСОБИЕ Гаспарян Олег Николаевич Д.т.н, с.н.с 2005 СОДЕРЖАНИЕ Система математических ...
ZI = interp2(X, Y, Z, XI, YI, method)
Если Y содержит комплексные числа, то запись plot(Y) эквивалентна записи plot(real(Y), imag(Y)). Во всех остальных случаях, при использовании plot мнимые компо-ненты игнорируются.
Раздел: Рефераты по информатике, программированию
Тип: учебное пособие
Проектирование трансляторов
ЛЕКЦИЯ 1 СУЩНОСТЬ ПРЕДМЕТА. СОДЕРЖАНИЕ КП. СРОКИ. ОРГАНИЗАЦИЯ РАБОТ. МАТЕМАТИЧЕСКИЙ АППАРАТ. СТРУКТУРНАЯ СХЕМА ТРАНСЛЯТОРА. ПРОХОДЫ ТРАНСЛЯТОРА ...
x=x0, y=xN, xI-1 -> xI (I=1,2,...,N).
(1) begin real x,y
Раздел: Рефераты по информатике, программированию
Тип: реферат