Contoh source code matlab untuk metode posisi palsu

Program Matlab dengan menggunakan Metode Posisi Palsu

Dengan soal sebagai berikut :

F(x) = sinx^2 + cosx^2

a=0, b=3, epsilon=0.01

Perintah program :

syms x;

f=sin(x^2)+cos(x^2);

a=0;

b=3;

et=0.01;

e=abs(b-a);

i=1;

disp(' i      a       b      c     f(a)    f(b)   f(c)    E');

disp('----------------------------------------------------------');

clama=a;

cbaru=b;

while (e > et ) & (clama ~= cbaru);

    fa=subs(f,x,a);

    fb=subs(f,x,b);

    %c=(a+b)/2;

    clama=cbaru;

    c=(fb*a-fa*b)/(fb-fa);

    cbaru=c;

    fc= subs(f,x,c);

    fprintf('%3.0f %6.4f %6.4f %12.10f %7.4f %7.4f %7.4f %7.4f \n', i, a, b, c, fa, fb, fc, e);

    if fa*fc < 0

    b=c; %geser kiri

    else

      a=c; %geser kanan

    end

    e=abs(b-a);  % menghitung error

    i=i+1;

end

akar = c;

fprintf('\n\nAkarnya adalah = %8.7f\n', akar);

fprintf('Jumlah Iterasi = %g\n', i);