VPF::[Pascal] Вычислить значение функций

Kuvaldis

Дата 14.10.2006, 12:26 (ссылка)

(нет голосов)

Загрузка ...

механик-вредитель

Профиль
Группа: Участник Клуба
Сообщений: 1189
Регистрация: 16.6.2006
Где: Минск

Репутация: 32
Всего: 61

Koschei,

alexeis1 прав, даже не обязательно вычислять каждый раз x^n = exp(ln(x)) . Можно (и нужно) заметить, что x^n = x * x^(n-1)
Но если тебе уж очень приспичило, то вот вариант с "полноценным" вычислением факториала

Код


{$N+}

const
  err = 1e-4;
//******************************************************************************
function factorial(n : integer) : integer;
var
   f, i : integer;
begin
   f := 1;
   for i := 2 to n do
       f := f * i;
   factorial := f;    
end;
//******************************************************************************
function f(x : double) : double;
var
  n    : integer;
  s, c : double;
  sign : integer;
  fact : double;
Begin
  sign := 1;
  s := 0;
  n := 1;

  repeat
      fact := factorial(2 * n - 1);
      c := sign * exp(n * ln(x)) / ((2 * n - 1) * fact);
      s := s + c;
      sign := -sign;
      inc(n);
  until abs(c) < err;

  f := s;
end;
//******************************************************************************
var
  x : double;
begin
  x := 0.2;
  Repeat
    Writeln('x= ', x : 1 : 1, '  f= ', f(x) : 2 : 4);
    x := x + 0.1    
  until x >= 2.001;    
  readln;
end.
//******************************************************************************

--------------------

Помни - когда ты спишь, враг не дремлет
Спи чаще и дольше, изматывай врага бессоницей

Alexeis

Дата 14.10.2006, 13:40 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Да это сложнее и рекурсивный алгоритм как вегда, для тех кто не боится мозг свернуть.

Код


program Project2;

{$N+}

const
  err = 1e-4;

var
  sign : integer;
  fact : double;
  level: Integer;

function Recur(x : Double; n : integer) : Double;
var
  c : double;

begin
  if n = 1
  then
    begin
      sign := 1;
      fact := 3 * 2;
      n := 1;
      c := sign * exp(n * ln(x)) / ((2 * n - 1) * fact);

      if Abs(c) > Err
      then
        Recur := c + Recur(x, n + 1)
      else
        begin
          Recur := c;
          level := 1;
        end;
    end
  else
    begin
      c    := sign * exp(n * ln(x)) / ((2 * n - 1) * fact);
      sign := -sign;
      fact := fact * (2 * (n + 1)) * (2 * (n + 1) + 1);

      if Abs(c) > Err
      then
        Recur := c + Recur(x, n + 1)
      else
        begin
          Recur := c;
          level := n;
        end;
    end;
end;

var
  x : double;
begin
  x := 0.2;
  Repeat
    Writeln('x= ', x : 1 : 1, '  f= ', Recur(x, 1) : 2 : 4, ' n=', level);
    x := x + 0.1
  until x >= 2.001;
  readln;
end.

Модератор: На будщее читайте правила.

Правила (пункт 10): Если вы хотите, чтобы вашу проблему решили при помощи определенного

алгоритма, то не забудьте описать его!

Это сообщение отредактировал(а) alexeis1 - 14.10.2006, 13:41

--------------------

Vit вечная память.

Обсуждение действий администрации форума производятся только в этом форуме

гениальность идеи состоит в том, что ее невозможно придумать

Koschei

Дата 15.10.2006, 06:27 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

Спасибо,теперь надеюсь препода это устроит.А то виделите подрозумевалось таким образом решать именно...

Koschei

Дата 19.10.2006, 10:42 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

блин мне не рекурсия нужна,а рекуретный коофициент вычислить!вообшем вычислить значение функций с помошью рекуретного алгоритма!помогите,а?

Alexeis

Дата 19.10.2006, 13:34 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Вы противоречите сами себе

Цитата(Koschei @ 19.10.2006, 10:42

)

блин мне не рекурсия нужна

Цитата(Koschei @ 19.10.2006, 10:42

)

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

Это одно и тоже. Раберитесь наконец с формулировкой задачи, потому как узнать о чем думает ваш преподаватель я точно не в силах. smile

--------------------

Koschei

Дата 19.10.2006, 18:47 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

задача ток решена неправильно smile

этим рекуретным алгоритмом.С ответом прошлого решения не сходится.а там ответ правильный ибо преподаватель ответ проверял..

Alexeis

Дата 19.10.2006, 20:01 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Код


program Project1;

{$N+}

const
  err = 1e-4;

var
  sign : integer;
  fact : double;
  level: Integer;

function Recur(x : Double; n : integer) : Double;
var
  c : double;

begin
  if n = 1
  then
    begin
      sign := 1;
      fact := 3 * 2;
      n := 1;
      c := sign * exp(n * ln(x)) / ((2 * n - 1) * fact);
      sign := -sign;
      fact := fact * (2 * (n + 1)) * (2 * (n + 1) + 1);

      if Abs(c) > Err
      then
        Recur := c + Recur(x, n + 1)
      else
        begin
          Recur := c;
          level := 1;
        end;
    end
  else
    begin
      c    := sign * exp(n * ln(x)) / ((2 * n - 1) * fact);
      sign := -sign;
      fact := fact * (2 * (n + 1)) * (2 * (n + 1) + 1);

      if Abs(c) > Err
      then
        Recur := c + Recur(x, n + 1)
      else
        begin
          Recur := c;
          level := n;
        end;
    end;
end;

var
  x : double;
begin
  x := 0.2;
  Repeat
    Writeln('x= ', x : 1 : 1, '  f= ', Recur(x, 1) : 2 : 4, ' n=', level);
    x := x + 0.1
  until x >= 2.001;
  readln;
end.

--------------------

Koschei

Дата 27.10.2006, 12:39 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

А не могли бы вы написать текстовый алгоритм?как на той странице smile

буду очень признателен!

Alexeis

Дата 27.10.2006, 14:02 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Koschei, так ты бы сказал, что конкретно надо! Выводится и так все в текстовом виде. Побольше конкретики. Откуда же мне знать, в каком виде у вас там принято оформлять.

--------------------

Koschei

Дата 27.10.2006, 14:42 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

Не вы не поняли уважаемые,под текстовым алгоритмом я имел ввиду это:
Текстовый алгоритм:
Функция f:

1. Инициализируем локальные переменные
2. Вычисляем текущий член суммы
3. Добавляем его к суммы
4. Если текущий член суммы больше заданной точности, то переход на шаг 2.
5. Возвращаем в вызывающую функцию полученную сумму.

Главная функция

1. Инициализируем переменные
2. Во время выполнения цикла выводим аргумент и значение функции на экран для заданного интервала

и для полного счастья:
Не могли бы вы уважаемые написать что каждый шаг делает?
Например:
var
  n    : integer;
  s, c : double;         = здесь идет перечисление "переменных"
  sign : integer;
  fact : double;
А то как мне обьснить преподу что и как у меня получилось?
p.s.а точно прокатит рекурсия?"если препод точно имел ввиду рекуретность?"
А не могли бы вы как-нибудь по простому сделать.чтоб было видно что первокурсник сделал,а то как вы сделали у нас на курсе таких команд не знают даже smile

Alexeis

Дата 27.10.2006, 15:02 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Цитата(Koschei @ 27.10.2006, 14:42

)

А не могли бы вы как-нибудь по простому сделать.чтоб было видно что первокурсник сделал,а то как вы сделали у нас на курсе таких команд не знают даже

Каких команд не знают?

Вообще-то рекурсия вещь непростая сама по себе, и алгоритмы получаются хочешь не хочешь, закрученными. Я сам мозг чуть не свернул smile

.

Давайте вначале с командами разберемся.

--------------------

Koschei

Дата 27.10.2006, 15:07 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

Цитата

Каких команд не знают?

ну например эти
recur
sign

Alexeis

Дата 27.10.2006, 15:26 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

recur - это не команда, а просто имя функции которая вызывает сама себя т.е. рекурсивной. Его можно заменить хоть на Vasea_Pupkin smile

sign = это имя переменной, которая хранит значения знака члена ряда. (он же чередуется плюс / минус), вот она и меняет значение -1 +1.

Какие еще названия не ясны?

--------------------

Koschei

Дата 27.10.2006, 18:51 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

Цитата

Какие еще названия не ясны?

ну с остальными ясно.
Мне бы текстовой алгоритмик и что действия обозначают;)А то без этого мне серавно не сдать лабу...
p.s.задача не решена правильно ибо опять ответ с первым решением(тот который не рекуретным сделан) не сходится smile

(а там решено правильно),да еще какая то n появилось в ответе.это вообше чего такое?

Это сообщение отредактировал(а) Koschei - 27.10.2006, 19:11

Koschei

Дата 28.10.2006, 10:42 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 32
Регистрация: 22.9.2006

Репутация: нет
Всего: нет

Товариши не бросайте меня на пол пути.помогите сделать все правильно! smile

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »