Script Matlab Metode titik tengah

clear all;

clc;

disp('*************************************************************************************************        ')

disp('*                                     METODE NUMERIK                                              *')

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

disp('*                    Menghitungakar - akarpersamaandengan c                 *')

disp('*                   JurusanPendidikanMatematika  fakultasTarbiyahdanKeguruan               *')

disp('*                                IAIN Antasari Banjarmasin                                      *')

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

disp('* Oleh  :                                                                                       *')

disp('* Kelompok      6                                                                               *')

disp('* Hairianor(1201250929)                                                                  *')

disp('* Muhammad Arif Rahman(1201250938)                                                                          *')

disp('* Kartini  (1201250859)                                                                        *')

disp('* Hardiati Rizki Amelia  (1201250850)                                                                        *')

disp('* Lika Wartita  (1201250864)                                                                        *')

disp('*************************************************************************************************        ')

disp('*                                                                                               *')

syms x;

 f=input('masukkan persamaan f(x): ');

    a=input('masukkan nilai a : ');

    b=input('masukkan nilai b : ');

    et=input('masukkan Error Toleransi : ');

    e=abs(b-a)/b;

    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(cbaru-clama)/clama;  % menghitung error

        i=i+1;

    end

Semoga Bermanfaat......