Версия для печати темы

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

Форум «Всё о Паскале» _ Задачи _ Вычисление факториала методом Стирлинга

Автор: Гость_Andrew 20.10.2004 18:13

Собсно сабж.

Автор: Гость_Andrew 20.10.2004 18:19

Формула стерлинга выглядит так

n!=квадратный корень_из_2*Пи*n*n^n*e^n


у меня получилось так (сразу предупрежу что я корявый и с паскалём лет 5 не сталкивался smile.gif

Цитата
program vova;
var n,o:real;
begin
Writeln('Enter the N number');
Readln(n);
o:=sqrt(2)*3.14*n*exp(n*ln(n))*exp(n*ln(2.718));
Writeln('factorial n!=',o);
Readln;
end.


Вроде что то считает (не уверен что правильно)
ввожу число 99 (и больше) - пишет ошибка 205: Переполнение плавающей запятой.
То есть вся суть пользования формулой стирлина пропадает. А пользуют её для приближённого вычисления факториалов больших чисел.

Ткните где что поправить можно.

Автор: Atos 20.10.2004 18:31

По-моему эта формула должна выглядеть так:
Квадратный корень(2*pi*n) *((n/e)^n)

Автор: Amro 20.10.2004 18:57

А по моему так, без всяких тама корней, хотя может я и ошибаюсь
n!= (n/e)^n*sqrt(2*pi*n)
или лучше так n! := Exp(n*Ln(n/Exp(1))) * Sqrt(2*Pi*n);
Ну это ладно!!!
Andrew А чтоб ошибка не возникала, напиши начало так

Код
{$N+} {для подклучения всех вещественных типов}
type Real = Extended; {берём более большой вещест тип}
var n,o:Real;
......................

Автор: Atos 20.10.2004 19:00

Amro Sqrt - это и есть квадратный корень smile.gif

Автор: Amro 20.10.2004 19:04

Цитата
Квадратный корень(2*pi*n) *((n/e)^n)

Atos В смысле не всё выражение под корнем!!!!

Автор: Atos 20.10.2004 19:07

А, понятно... Я так и имел в виду - под корнем всё, что в первых скобках

Автор: Amro 20.10.2004 19:09

ИМХО Значит я не доглядел
smile.gif Пардон

Автор: Altair 20.10.2004 22:14

так, вы с формулой не ошиблись?
или мне разобраться? smile.gif ;)

Автор: viagra online without prescripti 1.09.2021 3:46

Buy Std Medication Online Canada

Автор: how much does viagra cost at wal 5.11.2021 19:03

puedo comprar cialis sin receta

Автор: ivermectin tablet online purchas 11.11.2021 18:47

Viagra Paypal Accepted