Sabtu, 31 Januari 2015

Metode Newton Gregory Backward dengn Delphi7


procedure TForm3.Button1Click(Sender: TObject);
var
x0,x1,x2,x3,x4,x5,x6:real;
f0,f1,f2,f3,f4,f5,f6:real ;
d0,d1,d2,d3,d4,d5:real;
dd0,dd1,dd2,dd3,dd4:real;
ddd0,ddd1,ddd2,ddd3:real;
dddd0,dddd1,dddd2:real;
ddddd0,ddddd1:real;
dddddd0:real;
s,x,p6x:real;

begin
x6:=StrToFloat(Edit1.Text);
x5:=StrToFloat(Edit2.Text);
x4:=StrToFloat(Edit3.Text);
x3:=StrToFloat(Edit4.Text);
x2:=StrToFloat(Edit5.Text);
x1:=StrToFloat(Edit6.Text);
x0:=StrToFloat(Edit7.Text);

f6:=StrToFloat(Edit8.Text);
f5:=StrToFloat(Edit9.Text);
f4:=StrToFloat(Edit10.Text);
f3:=StrToFloat(Edit11.Text);
f2:=StrToFloat(Edit12.Text);
f1:=StrToFloat(Edit13.Text);
f0:=StrToFloat(Edit14.Text);

x:=StrToFloat(Edit36.Text);

//mencari DEL F

d0:=f1-f2;
d1:=f2-f1;
d2:=f3-f2;
d3:=f4-f3;
d4:=f5-f4;
d5:=f6-f5;
Edit20.Text:=(Format('%8.4f',[d0]));
Edit19.Text:=(Format('%8.4f',[d1]));
Edit18.Text:=(Format('%8.4f',[d2]));
Edit17.Text:=(Format('%8.4f',[d3]));
Edit16.Text:=(Format('%8.4f',[d4]));
Edit15.Text:=(Format('%8.4f',[d5]));

// Mencari DEL^2F

dd0:=d1-d0;
dd1:=d2-d1;
dd2:=d3-d2;
dd3:=d4-d3;
dd4:=d5-d4;
Edit25.Text:=(Format('%8.4f',[dd0]));
Edit24.Text:=(Format('%8.4f',[dd1]));
Edit23.Text:=(Format('%8.4f',[dd2]));
Edit22.Text:=(Format('%8.4f',[dd3]));
Edit21.Text:=(Format('%8.4f',[dd4]));

//mencari DEL^3F

ddd0:=dd1-dd0;
ddd1:=dd2-dd1;
ddd2:=dd3-dd2;
ddd3:=dd4-dd3;
Edit29.Text:=(Format('%8.4f',[ddd0]));
Edit28.Text:=(Format('%8.4f',[ddd1]));
Edit27.Text:=(Format('%8.4f',[ddd2]));
Edit26.Text:=(Format('%8.4f',[ddd3]));

//mencari DEL^4F

dddd0:=ddd1-ddd0;
dddd1:=ddd2-ddd1;
dddd2:=ddd3-ddd2;
Edit32.Text:=(Format('%8.4f',[dddd0]));
Edit31.Text:=(Format('%8.4f',[dddd1]));
Edit30.Text:=(Format('%8.4f',[dddd2]));

//mencari DEL^5F

ddddd0:=dddd1-dddd0;
ddddd1:=dddd2-dddd1;
Edit34.Text:=(Format('%8.4f',[ddddd0]));
Edit33.Text:=(Format('%8.4f',[ddddd1]));

//mencari DEL^6F

dddddd0:=ddddd1-ddddd0;
Edit35.Text:=(Format('%8.4f',[dddddd0]));

//mencari s

s:=(x-x0)/(x0-x1);
Edit37.Text:=(format('%8.4f',[s]));

//mencari P6X

p6x:=(f0)+(s/1)*d0+((s*(s+1))/2)*dd0+((s*(s+1)*(s+2))/6)*ddd0+((s*(s+1)*(s+2)*(s+3))/24)*dddd0+((s*(s+1)*(s+2)*(s+3)*(s+4))/120)*ddddd0+((s*(s+1)*(s+2)*(s+3)*(s+4)*(s+5))/720)*dddddd0;


Edit38.Text:=(format('%8.4f',[p6x]));


Selamat Mencoba, Semoga Bermanfaat. 

Tidak ada komentar:

Posting Komentar