Форум по Delphi программированию

Delphi Sources



Вернуться   Форум по Delphi программированию > Все о Delphi > [ "Начинающим" ]
Ник
Пароль
Регистрация <<         Правила форума         >> FAQ Пользователи Календарь Поиск Сообщения за сегодня Все разделы прочитаны

 
 
Опции темы Поиск в этой теме Опции просмотра
  #1  
Старый 15.12.2011, 01:20
Аватар для Kpyc
Kpyc Kpyc вне форума
Новичок
 
Регистрация: 05.05.2011
Сообщения: 75
Версия Delphi: Turbo Delphi
Репутация: 11
Вопрос Интеграл вероятности

Не хочет его считать, методом проб и ошибок выяснил, что вываливается он на n=34 с ошибкой деления на 0.
Полный листинг:
Код:
unit Unit2;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ExtCtrls, Math;

type
  TfrmMain = class(TForm)
    edtX: TLabeledEdit;
    lblResult: TLabel;
    btnCalc: TButton;
    function Summ(x: Extended): Extended;
    function Fact(n: integer): integer;
    procedure btnCalcClick(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  frmMain: TfrmMain;
  x: extended;

implementation

{$R *.dfm}

procedure TfrmMain.btnCalcClick(Sender: TObject);
begin
  try
    x := strtofloat(edtX.Text);
    lblResult.Caption := 'erf(x) = ' + floattostr(summ(x));
  except
    lblResult.Caption := 'Ошибка!';
  end;
end;

function TfrmMain.Fact(n: integer): Integer;
begin
  if n > 1 then
    Result := n * Fact(n - 1)
  else
    Result := 1;
end;

function TfrmMain.Summ(x: Extended): Extended;
var
  n: Integer;
begin
  Result := 0;
  for n := 0 to 33 do
      Result := Result + power( -1, n) * power(x, 2*n + 1)/(fact(n) * (2*n + 1));
  Result := Result * 2 / sqrt(pi);
end;

end.

Как лечить?
Ответить с цитированием
 


Delphi Sources

Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск
Опции просмотра

Ваши права в разделе
Вы не можете создавать темы
Вы не можете отвечать на сообщения
Вы не можете прикреплять файлы
Вы не можете редактировать сообщения

BB-коды Вкл.
Смайлы Вкл.
[IMG] код Вкл.
HTML код Выкл.
Быстрый переход


Часовой пояс GMT +3, время: 15:28.


 

Сайт

Форум

FAQ

Соглашения

Прочее

 

Copyright © Форум "Delphi Sources" by BrokenByte Software, 2004-2025