Рассчитать N значений y(x) на интервале от а до b

[Юлия67]

Рассчитать N значений y(x) на интервале х принадлежит [а;b]. Каждое значение функции рассчитывать с точностью эпсилон (е).

[Юлия67]

Program Lena_1;
uses crt;
function Fact(N:word):real;
begin
  if N<=1 then Fact:=1
          else Fact := N*Fact(N-1);
  end;
function Exp1(x, e:real):real;
  var s, t:real;
  i:word;
  begin
   i:=7;
   s:=2*x/Fact(2*i+1);
   t:=1;
   while abs(t)>=e do
  begin
   inc(i);
   t:=2*x/Fact(2*i+1);
   s:=s+t;
  end;
Exp1:=S;
end;

var
  a,b,N,c,j:integer;
  x,y,d:real;
  e:real;

BEGIN
clrscr;
repeat
write('vvedite tochnost 0<E<1 E:=');
readln(e);
until (e>0)and(e>1);
writeln('Ukaghite otrezok glya x');
repeat
write('a:=');
readln(a);
write('b:=');
readln(b);
Writeln ('N:=');
Readln(N);
until a<>b;
if a>b then
   begin
   c:=a; a:=b; b:=c
   end;
writeln('x prinadleghit otrezky [',a,';',b,']');
x:=a; j:=1;

d:=(b-a)/(N-1);
writeln;
while x<=b+0.1 do
   begin
   if x<>0 then
   writeln('x',j,'=',x:5:2,' y(x)=', Exp1(x,e):0:(N-1),'s tochnostyu E=',e:0:(N-1));
    inc(j);
    x:=x+d
    end;
readln
END.

помогите пожалуйста подкорректировать программу, где то ошибка не могу понять где

[Stilet]

Каждое значение функции рассчитывать с точностью эпсилон (е).

Что означает эта фраза?
Как можно функцию рассчитывать с какой-то точностью? И зачем?

[Юлия67]

сама не понимаю, такое вот задание дали...

[BDA]

Как можно функцию рассчитывать с какой-то точностью? И зачем?

Вот есть функция, равная сумме бесконечного ряда.
Поскольку факториал растет гораздо быстрее линейной функции от х, то очень быстро очередной член ряда оказывается меньше наперед заданного эпсилон.

Чуть-чуть исправил код:

Program Lena_1;
uses crt;
function Fact(N: word): real;
begin
  if N <= 1 then
    Fact := 1
  else
    Fact := N * Fact(N - 1);
end;

function Exp1(x, e: real): real;
var
  s, t: real;
  i: word;
begin
  i := 7;
  s := 0;
  t := 2 * x / Fact(2 * i + 1);
  while abs(t) >= e do
  begin
    s := s + t;
    inc(i);
    t := t / (2 * i + 1);
  end;
  Exp1 := s;
end;

var
  N, j: integer;
  a, b, c, x, y, e: real;

BEGIN
  clrscr;
  repeat
    write('vvedite tochnost 0<E<1 E:=');
    readln(e);
  until (e > 0) and (e < 1);
  writeln('Ukaghite otrezok glya x');
  repeat
    write('a:=');
    readln(a);
    write('b:=');
    readln(b);
    write('N:=');
    readln(N);
  until a <> b;
  if a > b then
  begin
    c := a;
    a := b;
    b := c;
  end;
  writeln('x prinadleghit otrezky [', a:0:0, ';', b:0:0, ']');
  x := a;
  j := 1;
  writeln;
  while x <= b do
  begin
    writeln('x', j, '=', x:5:2, ' y(x)=', Exp1(x, e):0:10, ' s tochnostyu E=',
      e:0:10);
    inc(j);
    x := a + (b - a) * (j - 1) / (N - 1);
  end;
  readln;
END.

[Юлия67]

отлично, спасибо, а можно еще вопрос? почему у меня y(x) всегда 0 равен??

[BDA]

Пожалуйста.
Потому что факториал в первом члене равен 15! = 1307674368000.
Для теста я ввел a = 1307674368000, b - любое число, чуть большее a, и получил y, равный 2.11 (по-моему, вполне реальный результат).

[Юлия67]

все ясно)) спасибо огромное еще раз)