Mencari nilai dari PI mengunakan "Nilakantha Series".

Pi (π). merupakan salah satu bilangan unik yang dihasilkan dari perbandingan antara keliling lingkaran dengan diameter yang angkanya selalu konstan yaitu sekitar(mendekati) 3,14. Nah dengan bantuan program Scilab (saudaranya Matlab ). kita bisa mencari nilai PI dengan "Nilakantha Series" tadi. Berikut persamaan dari "Nilakantha Series" sendiri:

sumber : https://www.mathsisfun.com/numbers/pi.html

 nah, untuk ilustrasi dari pengambilan nilai PI  sendiri bisa dilihat di bawah :

sumber : https://www.mathsisfun.com/numbers/pi.html

Oke, langsung saja ke Codingan Scilab untuk Nilakantha Series. Saya sendiri masih terinspirasi dari codingan dua artikel ini mencari nilai sin(x) dan mencari nilai cos(x), dalam pembuatan codingan kali ini dan berikut codingannya :

clear;
f0=3;
f1=f0+(4/(2*3*4));
Galat = abs(f1-f0)/f1*100;

j=4;
n=1;
y=3;
disp([n,f1,Galat]);

for i=1:9
    n=n+1;
    f0=f1;
    c=modulo(y,2);
    if (c==1) then
        f1 = f0 - (4/(j*(j+1)*(j+2)));
    else
        f1 = f0 + (4/(j*(j+1)*(j+2)));
    end
    
    Galat = abs(f1-f0)/f1*100;
    disp([n,f1,Galat,c,j]);
    
    j=j+2;
    y=y+1;
        
end

 

nah untuk penjelasannya sebagai berikut :

clear;

f0=3;

f1=f0+(4/(2*3*4));

Galat = abs(f1-f0)/f1*100;

clear; untuk membersihkan serta mengosongkan memory.

f0=3; menginialisasi nilai 3 ke variabel f0.

f1=f0+(4/(2*3*4)); f1 diberi nilai sesuai dengan rumus "Nilakantha Series".

Galat = abs(f1-f0)/f1*100; variabel Galat digunakan untuk proses seberapa error kah nilai yang dihasilkan.

j=4;

n=1;

y=3;

disp([n,f1,Galat]);

j=4; memberi nilai 4 ke variabel j. variabel J sendiri digunakan untuk penggantian angka 4,6,8;

n=1; memberi nilai 1 ke variabel n. variabel n sendiri digunakan sebagai petunjuk proses keberapa.

y=3; memberi nilai 3 ke variabel y. variabel y sendiri digunakan untuk petunjuk apakah penjumlahan atau pengurangan untuk proses selanjutnya.

disp([n,f1,Galat]); digunakan untuk menampilkan n,f1,Galat.

for i=1:9
    n=n+1;
    f0=f1;
    c=modulo(y,2);
    if (c==1) then
        f1 = f0 - (4/(j*(j+1)*(j+2)));
    else
        f1 = f0 + (4/(j*(j+1)*(j+2)));
    end
    
    Galat = abs(f1-f0)/f1*100;
    disp([n,f1,Galat,c,j]);
    
    j=j+2;
    y=y+1;
        
end

nah. untuk bagian ini sendiri sebagai looping yang dilakukan sebanyak 9 kali tanpa memperhatikan faktor Bena.

dan berikut hasilnya :

 Semoga bermanfaat dan membantu.

 

Lokasi:

Berbagi :