function g1(v1:integer):real;
var fak:real;
begin
fak:=1;
for i:=1 to (v1-1) do
fak:=fak*i;
g1:=fak;
end;
function
g2(v2:integer):real;
var fak,k:real;
begin
if v2=0 then
g2:=sqrt(3.14)
else
k:=-0.5;
fak:=1;
for i:=1 to v2 do
fak:=fak*(k+i);
g2:=fak*sqrt(3.14);
end;
function
pangkat(x,p:real):real;
begin
pangkat:=exp(p*ln(x));
end;
function
fx1(v:integer):real;
var
h,hi,p,x,xo,xn,fungsi,fungsi1,fungsi2,fsigma:real;
begin
v1:=(v+1) div 2;
v2:=v div 2;
h:=(t)/1000;
xo:=(1);
xn:=(1+(t*t/v));
p:=(v+1)/2;
fungsi:=g1(v1)/(g2(v2)*sqrt(v*3.14));
fungsi1:=fungsi/(pangkat(xo,p));
fungsi2:=fungsi/(pangkat(xn,p));
fsigma:=0;
for i:=1 to 999 do
begin
hi:=i*h;
x:=(1+(sqr(hi)/v));
fsigma:=fsigma+(fungsi/(pangkat(x,p)));
end;
fx1:=1-(0.5+(h/2)*(fungsi1+(2*fsigma)+fungsi2));
end;
function
fx2(v:integer):real;
var
h,hi,p,x,xo,xn,fungsi,fungsi1,fungsi2,fsigma:real;
begin
v1:=v div 2;
v2:=(v+1) div 2;
h:=(t)/1000;
xo:=(1);
xn:=(1+(t*t/v));
p:=(v+1)/2;
fungsi:=g2(v1)/(g1(v2)*sqrt(v*3.14));
fungsi1:=fungsi/(pangkat(xo,p));
fungsi2:=fungsi/(pangkat(xn,p));
fsigma:=0;
for i:=1 to 999 do
begin
hi:=i*h;
x:=(1+(sqr(hi)/v));
fsigma:=fsigma+(fungsi/(pangkat(x,p)));
end;
fx2:=1-(0.5+(h/2)*(fungsi1+(2*fsigma)+fungsi2));
end;
{$R *.dfm}
procedure
TForm1.Button2Click(Sender: TObject);
begin
n:=StrToInt(Edit1.Text);
x:=StrToFloat(Edit2.Text);
s:=StrToFloat(Edit3.Text);
x0:=StrToFloat(Edit4.Text);
alpha:=StrToFloat(Edit5.Text);
t:=(x-x0)/(s/sqrt(n));
v:=n-1;
d:=ComboBox1.Text;
sisa:=v mod 2;
if sisa=1 then
pvalue:=fx1(v)
else
pvalue:=fx2(v);
if (d='greater than')
then
pvalue:=pvalue
else if (d='less than')
then
pvalue:=1-pvalue
else
begin
if t<0 then t:=-t;
sisa:=v mod 2;
if sisa=1 then
pvalue:=fx1(v)
else
pvalue:=fx2(v);
pvalue:=2*pvalue;
end;
label7.caption:=FloatToStrf(pvalue,ffnumber,4,3);
if pvalue<alpha then
label8.caption:='tolak
H0. mu not equals '+FloatToStr(x0)
else
label8.caption:='gagal
tolak H0. mu equals '+FloatToStr(x0);
end;
