Metode Bisection dengan Matlab

Ini adalah contoh script metode bisection menggunakan matlab. Silahkan dianalisa dan diedit sendiri.


'Ujian Matlab 8.36';
'Metode Bisection';
'fungsi y=tan(Te)*x-g*x^2/(2*Vo^2*cos(Te)^2)+yo';
x=35;
g=9.81;
Vo=20;
yo=2;
Te1=input('masukkan tebakan kiri  (dalam derajat) : ');
Te2=input('masukkan tebakan kanan (dalam derajat) : ');
Te1=Te1/180*pi;
Te2=Te2/180*pi;
Te3=Te1;
z1=tan(Te1)*x-g*x^2/(2*Vo^2*cos(Te1)^2)+yo-1;
z2=tan(Te2)*x-g*x^2/(2*Vo^2*cos(Te2)^2)+yo-1;
if z1<0;
    Te1=Te1;
else Te1=Te2;
    Te2=Te3;
end
if z1*z2<0;
    for i=1:100;
    y=(Te1+Te2)/2;
    if tan(y)*x-g*x^2/(2*Vo^2*cos(y)^2)+yo-1<0;
        Te1=y;
    else Te2=y;
    end
    end
    Hasil=Te1/pi*180
else 'akar tidak diapit';
end

x=0:0.1:35;
y=tan(Te1).*x-g.*x.^2./(2*Vo^2*cos(Te1)^2)+yo;
plot(x,y);

Posted in: download, script program