Love What You Do, Do What You Love

Tampilkan postingan dengan label Metode Numerik. Tampilkan semua postingan
Tampilkan postingan dengan label Metode Numerik. Tampilkan semua postingan

Program Mencari Akar Hampiran dengan Metode Iterasi Titik Tetap Menggunakan Delphi 7

Posted by : Unknown
Selasa, 26 November 2013 0 comments

Pseudo Code:

var
x0,x1,e,p,q:real;
i,n:integer;
begin
  x0:=strtofloat(edit1.Text);
  e:=0.0001;
  i:=0;
  n:=strtoint(edit2.Text);
  repeat
    i:=i+1;
    listbox1.Items.add(inttostr(i));
    listbox2.Items.add(floattostr(x0));
    x1:=sqrt(cos(x0));
    listbox3.Items.add(floattostr(x1));
    p:=(x1-x0)/x1;
    q:=abs(p);
    listbox4.Items.add(floattostr(q));
    x0:=x1;
  until
    q<e;
  edit3.Text:=floattostr(x1);
end;

end.

Program Mencari Akar Hampiran dengan Metode Secant Menggunakan Delphi 7

Posted by : Unknown
0 comments

Pseudo Code:

var
  x0,x1,x2,y0,y1,y2,e,d,m,j:real;
  i,n:integer;
begin
  x0:=strtofloat(edit1.text);
  x1:=strtofloat(edit2.Text);
  e:=0.000001;
  d:=0.000001;
  i:=0;
  n:=strtoint(edit3.Text);
repeat
  i:=i+1;
  listbox1.Items.Add(inttostr(i));
  listbox2.Items.add(format('%8.4f',[x0]));
  listbox3.items.Add(format('%8.4f',[x1]));
  y0:=(cos(x0)-(x0*x0));
  y1:=(cos(x1)-(x1*x1));
  listbox5.Items.add(format('%8.4f',[y0]));
  listbox6.Items.add(format('%8.4f',[y1]));
  j:=abs(y1-y0);
  if j<d then showmessage ('terlalukecil');
  x2:=((x0*y1)-(x1*y0))/(y1-y0);
  listbox4.Items.add(format('%8.4f',[x2]));
  y2:=(cos(x2)-(x2*x2));
  m:=abs(y2);
  listbox7.Items.add(format('%8.6f',[y2]));
  x0:=x1;
  x1:=x2;
  if i=n then m:=e;
until
  m<e;
  edit4.Text:=format('%8.4f',[x2]);
end;
end.

Program Mencari Akar Hampiran dengan Metode Posisi Palsu Menggunakan Delphi 7

Posted by : Unknown
1 comments

Pseudo Code:

var
a,b,e,x, c, fa, fb, fc : real;
i:integer;

begin
  a:= strtofloat (edit1.text);
  b:= strtofloat (edit2.Text);

  fa := cos(a) - (a*a);
  fb := cos(b) - (b*b);
    if fa*fb > 0 then  label1.Caption:='selang a dan b tidak memuat akar'
    else
    begin
    x:= strtofloat(edit3.Text);
    i:=-1;
    repeat
    i:= i+1;
    c:= (b- (((fb*b)-(fb*a))/(fb-fa)));
    listbox1.Items.add(floattostr(a));
    listbox2.Items.add(floattostr(b));
    listbox3.Items.add(floattostr(c));
    fa := cos(a) - (a*a);
    fb := cos(b) - (b*b);
    fc := cos(c) - (c*c);
    listbox4.Items.add(floattostr(fa));
    listbox5.Items.add(floattostr(fb));
    listbox6.Items.add(floattostr(fc));
    if fa*fc <= 0 then b:= c else a:=c;
    e:= abs ((b-a)/b);
    until e<= x;
    edit4.text:= floattostr(c);
    edit5.Text:= floattostr(i);
    end;
    end;

end.

Program Mencari Akar Hampiran dengan Metode Newton Raphson Menggunakan Delphi 7

Posted by : Unknown
0 comments

Pseudo Code:

var
  x0,x1,y0,y1,d,e,m:Real;
  i,n:integer;
begin
  x0:=StrToFloat(edit1.text);
  e:=0.000001 ;
  d:=0.000001 ;
  n:=strtoint(edit2.text);
  i:=0;
repeat
  i:=i+1;
  listbox1.Items.add(inttostr(i));
  listbox2.Items.add(format('%8.4f',[x0]));
  y0:=(cos(x0)-(x0*x0));
  y1:=(-(sin(x0))-2*x0);
  listbox4.Items.add(format('%8.4f',[y0]));
  listbox5.Items.add(format('%8.4f',[y1]));
  if (abs(y1))<d then showmessage('kemiringan terlalu kecil');
  x1:=(x0)-(y0/y1);
  listbox3.Items.add(format('%8.4f',[x1]));
  m:=abs((x1-x0)/x1);
  listbox6.Items.add(format('%8.4f',[m]));
  if m>e then x0:=x1;
until
  i=n;
  edit4.Text:=format('%8.4f',[x1]);
end;

end.

Program Mencari Akar Hampiran dengan Metode Bagi Dua Menggunakan Delphi 7

Posted by : Unknown
0 comments

Pseudo Code:

r
a,b,e,x, c, fa, fb, fc : real;
i:integer;

begin
  a:= strtofloat (edit1.text);
  b:= strtofloat (edit2.Text);

  fa := cos(a) - (a*a);
  fb := cos(b) - (b*b);
    if fa*fb > 0 then  label1.Caption:='selang a dan b tidak memuat akar'
    else
    begin
    x:= strtofloat(edit3.Text);
    i:=-1;
    repeat
    i:= i+1;
    c:=(a+b)/2;
    listbox1.Items.add(floattostr(a));
    listbox2.Items.add(floattostr(b));
    listbox3.Items.add(floattostr(c));
    fa := cos(a) - (a*a);
    fb := cos(b) - (b*b);
    fc := cos(c) - (c*c);
    listbox4.Items.add(floattostr(fa));
    listbox5.Items.add(floattostr(fb));
    listbox6.Items.add(floattostr(fc));
    if fa*fc <= 0 then b:= c else a:=c;
    e:= abs ((b-a)/b);
    until e<= x;
    edit4.text:= floattostr(c);
    edit5.Text:= floattostr(i);
    end;
    end;

end.

Program Newton-Gregory Mundur Menggunakan Delphi 7

Posted by : Unknown
0 comments

Pseudo Code:

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;
  x,s,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(edit15.Text);

  //Del f
  d0:=f1-f0;
  d1:=f2-f1;
  d2:=f3-f2;
  d3:=f4-f3;
  d4:=f5-f4;
  d5:=f6-f5;
  edit21.Text:=(format('%8.4f',[d0]));
  edit20.Text:=(format('%8.4f',[d1]));
  edit19.Text:=(format('%8.4f',[d2]));
  edit18.Text:=(format('%8.4f',[d3]));
  edit17.Text:=(format('%8.4f',[d4]));
  edit16.Text:=(format('%8.4f',[d5]));

  //Del^2 f
  dd0:=d1-d0;
  dd1:=d2-d1;
  dd2:=d3-d2;
  dd3:=d4-d3;
  dd4:=d5-d4;
  edit26.Text:=(format('%8.4f',[dd0]));
  edit25.Text:=(format('%8.4f',[dd1]));
  edit24.Text:=(format('%8.4f',[dd2]));
  edit23.Text:=(format('%8.4f',[dd3]));
  edit22.Text:=(format('%8.4f',[dd4]));

  //Del^3 f
  ddd0:=dd1-dd0;
  ddd1:=dd2-dd1;
  ddd2:=dd3-dd2;
  ddd3:=dd4-dd3;
  edit30.Text:=(format('%8.4f',[ddd0]));
  edit29.Text:=(format('%8.4f',[ddd1]));
  edit28.Text:=(format('%8.4f',[ddd2]));
  edit27.Text:=(format('%8.4f',[ddd3]));

  //Del^4 f
  dddd0:=ddd1-ddd0;
  dddd1:=ddd2-ddd1;
  dddd2:=ddd3-ddd2;
  edit33.Text:=(format('%8.4f',[dddd0]));
  edit32.Text:=(format('%8.4f',[dddd1]));
  edit31.Text:=(format('%8.4f',[dddd2]));

  //Del^5 f
  ddddd0:=dddd1-dddd0;
  ddddd1:=dddd2-dddd1;
  edit35.Text:=(format('%8.4f',[ddddd0]));
  edit34.Text:=(format('%8.4f',[ddddd1]));

  //Del^6 f
  dddddd0:=ddddd1-ddddd0;
  edit36.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]));
end;

end.

Program Interpolasi Newton-Gregory Maju Menggunakan Delphi 7

Posted by : Unknown
1 comments

Pseudo Code:

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;
  x,s,p6x :real;

begin
  x0:=StrToFloat(edit1.Text);
  x1:=StrToFloat(edit2.Text);
  x2:=StrToFloat(edit3.Text);
  x3:=StrToFloat(edit4.Text);
  x4:=StrToFloat(edit5.Text);
  x5:=StrToFloat(edit6.Text);
  x6:=StrToFloat(edit7.Text);

  f0:=StrToFloat(edit8.Text);
  f1:=StrToFloat(edit9.Text);
  f2:=StrToFloat(edit10.Text);
  f3:=StrToFloat(edit11.Text);
  f4:=StrToFloat(edit12.Text);
  f5:=StrToFloat(edit13.Text);
  f6:=StrToFloat(edit14.Text);

  x:=StrToFloat(edit36.Text);

  //Del f
  d0:=f1-f0;
  d1:=f2-f1;
  d2:=f3-f2;
  d3:=f4-f3;
  d4:=f5-f4;
  d5:=f6-f5;
  edit15.Text:=(format('%8.4f',[d0]));
  edit16.Text:=(format('%8.4f',[d1]));
  edit17.Text:=(format('%8.4f',[d2]));
  edit18.Text:=(format('%8.4f',[d3]));
  edit19.Text:=(format('%8.4f',[d4]));
  edit20.Text:=(format('%8.4f',[d5]));

  //Del^2 f
  dd0:=d1-d0;
  dd1:=d2-d1;
  dd2:=d3-d2;
  dd3:=d4-d3;
  dd4:=d5-d4;
  edit21.Text:=(format('%8.4f',[dd0]));
  edit22.Text:=(format('%8.4f',[dd1]));
  edit23.Text:=(format('%8.4f',[dd2]));
  edit24.Text:=(format('%8.4f',[dd3]));
  edit25.Text:=(format('%8.4f',[dd4]));

  //Del^3 f
  ddd0:=dd1-dd0;
  ddd1:=dd2-dd1;
  ddd2:=dd3-dd2;
  ddd3:=dd4-dd3;
  edit26.Text:=(format('%8.4f',[ddd0]));
  edit27.Text:=(format('%8.4f',[ddd1]));
  edit28.Text:=(format('%8.4f',[ddd2]));
  edit29.Text:=(format('%8.4f',[ddd3]));

  //Del^4 f
  dddd0:=ddd1-ddd0;
  dddd1:=ddd2-ddd1;
  dddd2:=ddd3-ddd2;
  edit30.Text:=(format('%8.4f',[dddd0]));
  edit31.Text:=(format('%8.4f',[dddd1]));
  edit32.Text:=(format('%8.4f',[dddd2]));

  //Del^5 f
  ddddd0:=dddd1-dddd0;
  ddddd1:=dddd2-dddd1;
  edit33.Text:=(format('%8.4f',[ddddd0]));
  edit34.Text:=(format('%8.4f',[ddddd1]));

  //Del^6 f
  dddddd0:=ddddd1-ddddd0;
  edit35.Text:=(format('%8.4f',[dddddd0]));

  //mencari S
  s:=(x-x0)/(x1-x0);
  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]));
end;

end.

Program Interpolasi Beda Terbagi Newton Menggunakan Delphi 7

Posted by : Unknown
0 comments

Pseudo Code:

var
  x0,x1,x2,x3,x4,x5,x6 :real;
  f0,f1,f2,f3,f4,f5,f6 :real;
  f01,f12,f23,f34,f45,f56 :real;
  f012,f123,f234,f345,f456 :real;
  f0123, f1234, f2345, f3456 :real;
  f01234, f12345, f23456 :real;
  f012345, f123456 :real;
  f0123456 :real;
  x,p6x :real;

begin
  x0:=StrToFloat(edit1.Text);
  x1:=StrToFloat(edit2.Text);
  x2:=StrToFloat(edit3.Text);
  x3:=StrToFloat(edit4.Text);
  x4:=StrToFloat(edit5.Text);
  x5:=StrToFloat(edit6.Text);
  x6:=StrToFloat(edit7.Text);

  f0:=StrToFloat(edit8.Text);
  f1:=StrToFloat(edit9.Text);
  f2:=StrToFloat(edit10.Text);
  f3:=StrToFloat(edit11.Text);
  f4:=StrToFloat(edit12.Text);
  f5:=StrToFloat(edit13.Text);
  f6:=StrToFloat(edit14.Text);

  //f[Xi,Xi+1]
  f01:=(f1-f0)/(x1-x0);
  f12:=(f2-f1)/(x2-x1);
  f23:=(f3-f2)/(x3-x2);
  f34:=(f4-f3)/(x4-x3);
  f45:=(f5-f4)/(x5-x4);
  f56:=(f6-f5)/(x6-x5);
  edit15.Text:=(format('%8.4f',[f01]));
  edit16.Text:=(format('%8.4f',[f12]));
  edit17.Text:=(format('%8.4f',[f23]));
  edit18.Text:=(format('%8.4f',[f34]));
  edit19.Text:=(format('%8.4f',[f45]));
  edit20.Text:=(format('%8.4f',[f56]));

  //f[X1,...,Xi+2]
  f012:=(f12-f01)/(x2-x0);
  f123:=(f23-f12)/(x3-x1);
  f234:=(f34-f23)/(x4-x2);
  f345:=(f45-f34)/(x5-x3);
  f456:=(f56-f45)/(x6-x4);
  edit21.Text:=(format('%8.4f',[f012]));
  edit22.Text:=(format('%8.4f',[f123]));
  edit23.Text:=(format('%8.4f',[f234]));
  edit24.Text:=(format('%8.4f',[f345]));
  edit25.Text:=(format('%8.4f',[f456]));

  //f[Xi,...,Xi+3]
  f0123:=(f123-f012)/(x3-x0);
  f1234:=(f234-f123)/(x4-x1);
  f2345:=(f345-f234)/(x5-x2);
  f3456:=(f456-f345)/(x6-x3);
  edit26.Text:=(format('%8.4f',[f0123]));
  edit27.Text:=(format('%8.4f',[f1234]));
  edit28.Text:=(format('%8.4f',[f2345]));
  edit29.Text:=(format('%8.4f',[f3456]));

  //f[Xi,...,Xi+4]
  f01234:=(f1234-f0123)/(x4-x0);
  f12345:=(f2345-f1234)/(x5-x1);
  f23456:=(f3456-f2345)/(x6-x2);
  edit30.Text:=(format('%8.4f',[f01234]));
  edit31.Text:=(format('%8.4f',[f12345]));
  edit32.Text:=(format('%8.4f',[f23456]));

  //f[X1,...,Xi+5]
  f012345:=(f12345-f01234)/(x5-x0);
  f123456:=(f23456-f12345)/(x6-x1);
  edit33.Text:=(format('%8.4f',[f012345]));
  edit34.Text:=(format('%8.4f',[f123456]));

  //f[Xi,...,Xi+6]
  f0123456:=(f12345-f01234)/(x6-x0);
  edit35.Text:=(format('%8.4f',[f0123456]));

  //mencari p6(x)
  x:=StrToFloat(edit36.Text);
  p6x:=(f0)+(f01*(x-x0))+(f012*(x-x0)*(x-x1))+(f0123*(x-x0)*(x-x1)*(x-x2))+(f01234*(x-x0)*(x-x1)*(x-x2)*(x-x3))+(f012345*(x-x0)*(x-x1)*(x-x2)*(x-x3)*(x-x4))+(f0123456*(x-x0)*(x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x5));
  edit37.Text:=(format('%8.4f',[p6x]));

end;

end.

Program Interpolasi lagrange Menggunakan Delphi 7

Posted by : Unknown
2 comments


Pseudo Code:
var
  x0,x1,x2,x3,x4,x5,x6 :real;
  y0,y1,y2,y3,y4,y5,y6 :real;
  L0x,L1x,L2x,L3x,L4x,L5x,L6x :real;
  x,p6x :real;

begin
  x0:=StrToFloat(edit1.Text);
  x1:=StrToFloat(edit2.Text);
  x2:=StrToFloat(edit3.Text);
  x3:=StrToFloat(edit4.Text);
  x4:=StrToFloat(edit5.Text);
  x5:=StrToFloat(edit6.Text);
  x6:=StrToFloat(edit7.Text);

  y0:=StrToFloat(edit8.Text);
  y1:=StrToFloat(edit9.Text);
  y2:=StrToFloat(edit10.Text);
  y3:=StrToFloat(edit11.Text);
  y4:=StrToFloat(edit12.Text);
  y5:=StrToFloat(edit13.Text);
  y6:=StrToFloat(edit14.Text);

  x:=STrToFloat(edit15.Text);

  //mencari Li(x)
  L0x:=((x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x5)*(x-x6))/((x0-x1)*(x0-x2)*(x0-x3)*(x0-x4)*(x0-x5)*(x0-x6));
  L1x:=((x-x0)*(x-x2)*(x-x3)*(x-x4)*(x-x5)*(x-x6))/((x1-x0)*(x1-x2)*(x1-x3)*(x1-x4)*(x1-x5)*(x1-x6));
  L2x:=((x-x0)*(x-x1)*(x-x3)*(x-x4)*(x-x5)*(x-x6))/((x2-x0)*(x2-x1)*(x2-x3)*(x2-x4)*(x2-x5)*(x2-x6));
  L3x:=((x-x0)*(x-x1)*(x-x2)*(x-x4)*(x-x5)*(x-x6))/((x3-x0)*(x3-x1)*(x3-x2)*(x3-x4)*(x3-x5)*(x3-x6));
  L4x:=((x-x0)*(x-x1)*(x-x2)*(x-x3)*(x-x5)*(x-x6))/((x4-x0)*(x4-x1)*(x4-x2)*(x4-x3)*(x4-x5)*(x4-x6));
  L5x:=((x-x0)*(x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x6))/((x5-x0)*(x5-x1)*(x5-x2)*(x5-x3)*(x5-x4)*(x5-x6));
  L6x:=((x-x0)*(x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x5))/((x6-x0)*(x6-x1)*(x6-x2)*(x6-x3)*(x6-x4)*(x6-x5));

  //mencari P6(x)
  p6x:=(y0*L0x)+(y1*L1x)+(y2*L2x)+(y3*L3x)+(y4*L4x)+(y5*L5x)+(y6*L6x);
  edit16.Text:=(format('%8.4f',[p6x]));

end;

Copyright © 2012 Zuka Zuka World | Hatsune Miku Theme | Designed by Johanes DJ