Форум программистов [Powered by Invision Power Board]

Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате

Форум программистов > Центр помощи > [Delphi] Вычислить: Sin(x)/x

Автор: Etrom 5.10.2012, 00:33

Помогите пожалуйста доделать задачу, если можно вместе с шапкой задачи

Вычислить: Sin(x)/x=1-x2/3!+x4/5!-x6/7!+...
1) Для N слагаемых
2) Для E>0 до тех пор пока слагаемое >E
3) Сравнить истинное значение.

Код


Procedure SinXdivX;
Var
    n:Word;
    F1, F2, P, k, Xn, Xk, dX, eps, fact : Real;
    sign : Boolean;
Begin
  Read(Xn,Xk);
  If Xn>Xk then
    Begin
      P:=Xn; Xn:=Xk; Xk:=P;
    end;
  eps:=0.001;
  dx:=(Xk-Xn)/10;
  While Xn<=Xk do
    Begin
      F2:=Sin(Xn)/Xn;
      n:=0;
      P:=1;
      F1:=1;
      fact:=1;
      k:=1;
      sign:=false;
      Repeat
        F2:=F1;
        fact:=fact*k;
        If sign then
          F1:=F1+P/fact
        Else
          F1:=F1-P/fact;
        sign:=not sign;
        inc(n);
      Until Abs(F1-F2)<eps;
      F2:=Sin(Xn)/Xn;
      WriteLn('|', Xn:12:6,'|', F1:12:6,'|',n:5,'|', F2:12:6,'|');
      Xn:=Xn+dx;
    end;

Автор: Pawl 6.10.2012, 11:43

2. Для заданой точности eps:

Код


program Project2;

{$APPTYPE CONSOLE}

uses
  SysUtils, Windows;

const eps = 0.001;

var a, s, x, u, p, y, z, dx : real;
  f : longint; n, k : integer;

begin
  setconsoleoutputcp(1251);
  write('введите начало и конец отрезка ');
  readln(y, z);
  dx := (z - y) / 10;
  x := y;
  writeln('x       сумма ряда      sin(x)');
  while x <= z do
  begin
    a := 1;
    f := 1;
    k := 1;
    s := 1;
    n := 1;
    repeat
      a := a * x * x;
      f := f * (2 * n) * ( 2 * n + 1);
      k := -k;
      u := k * a / f;
      s := s + u;
      inc(n);
    until abs(u) < eps;
    p := sin(x) / x;
    write(x:3:5, '    ');
    write(s:3:5, '    ');
    writeln(p:3:5);
    x := x + dx;
  end;
  readln;
end.

На отрезке [1..2] дает погрешность в 0,00001
1. Для заданного количества слагаемых:

Код


program Project2;

{$APPTYPE CONSOLE}

uses
  SysUtils, Windows;

var a, s, x, u, p, y, z, dx : real;
  f : longint; n, count, k : integer;

begin
  setconsoleoutputcp(1251);
  write('введите начало и конец отрезка ');
  readln(y, z);
  write('введите число слагаемых ');
  readln(count);
  dx := (z - y) / 10;
  x := y;
  writeln('x        сумма ряда      sin(x)');
  while x <= z do
  begin
    a := 1;
    f := 1;
    k := 1;
    s := 1;
    for n := 1 to count do
    begin
      a := a * x * x;
      f := f * (2 * n) * ( 2 * n + 1);
      k := -k;
      u := k * a / f;
      s := s + u;
    end;
    p := sin(x) / x;
    write(x:3:5, '    ');
    write(s:3:5, '    ');
    writeln(p:3:5);
    x := x + dx;
  end;
  readln;
end.

На отрезке [1..2] при count = 5 вообще не дает погрешности