Вычисление значения функции 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, задан массив: