metode biseksi menggunakan scilab

 

function []=biseksi()

clc;

a = input("Masukkan Nilai a=");

b = input("Masukkan nilai b=");

epsilon = input("Masukkan nilai Epsilon=");

E = abs(a-b);

i = 0;

disp("-----------------------------------------------------------------------")

disp ("   i        a            b         m          f(a)         f(m)         E ")

disp ("------------------------------------------------------------------------")

while mtlb_logic(E,">",mtlb_double(epsilon))

  i = i+1;

   mtlb_fprintf("%5.0f%12.7f%12.7f",i,a,b);

    fa = a^3+2*(a^2)+10*a-20;

  fb = b^3+2*(b^2)+10*b-20;

  m = b-(fb*(b-a))/(fb-fa)

  fm = m^3+2*(m^2)+10*m-20;

  if (fa*fm > 0) then

    a = m;

  else

    b = m;

  end;

  E = abs(b-a);

  mtlb_fprintf("%12.7f%12.7f%12.7f%12.7f\n",m,fa,fm,E);

end;

endfunction

simpan program ini dengan nama biseksi lalu 

ketik pada consol scilab biseksi


maka outputnya adalah