Matlab_Numerik

f=inline(‘x-2*sin(x)’,’x’);
disp(‘ ===============================================’)
disp(‘ METODE BISECTION ‘);
disp(‘ ===============================================’)
a = input(‘ Masukkan batas kiri interval : ‘ );
b = input(‘ Masukkan batas kanan interval : ‘ );
eps = input (‘ Nilai ketelitian (Toleransi) : ‘ );
itemax = input (‘ Masukkan jumlah iterasi maksimum : ‘ );
disp(‘ ‘)
disp(‘ ‘)
disp(‘ Memulai Proses Iterasi ‘)
disp(‘ ===============================================’)
disp(‘ k a b c selisih’)
disp(‘ ===============================================’)
selisih = abs(b-a);
D=a;E=b;M=(a + b)/2;F=selisih;
format short
tic
for k = 1:1:itemax
c = (a + b)/2;
disp ([k-1,a,b,c,selisih])
if (f(b)*f(c) <= 0)
a = c ;
else
b = c;
end
selisih = abs(b-a);
D=[D;a];E=[E;b];M=[M;c];F=[F;selisih];
if (abs(a-b)) <= eps
break
end
end
disp ([k,a,b,c,selisih])
disp(' ===============================================')
waktu = toc;
plot([D,E,M,F])
legend('a','b','c','selisih')
xlabel('Iterasi ke :')
ylabel('Iterasi metode Biseksi')
disp (' ')
disp ([' Akarnya adalah : ',num2str(c)])
disp ([' Jumlah iterasi : ',num2str(k)])
disp ([' Selang waktu konvergensi : ',num2str(waktu)])

ganti mawal dengan ini
f=inline('x-2*sin(x)','x'); df=inline('1-2*cos(x)','x');

3 Tanggapan

  1. halo,
    senang bertemu Anda melalui blog ini sy Agus Suhanto, tulisan yg menarik🙂 … salam kenal yee

  2. Mohom mangap mas Rahman, blogku sing bener iki http://www.wandisukoharjo.com , soale aku wis ura ngapdet sing wordpress kui tapi aktif neng dot come.😆

  3. salam balik, sori yo mas tulisan ini ntuk mahasiswa stmik yang kesulitan ngarap tugas, sebab kurang titik aja tidak jadi tuh, programnya dari dosen… dosennya nyomot dari dosen stmik lain… dst.
    jadi tidak menariiiiiiiik mas.
    salam kenal. tanx u

Komentar ditutup.

%d blogger menyukai ini: