Mirza Nur Hidayat

Home
Scilab
About
Visualisasi Sains - Animated Plot


Muatan q didistribusikan dalam cakram tipis yang jari-jari dalamnya a dan jari-jari luarnya b. Besar potensial listrik V di pusat cakram setelah melalui perhitungan matematis adalah V = kq (a+b) / 2ab dengan 0 < a < b. Guna menghasilkan visualisasi contour (gambar atas) dan model 3d (gambar bawah), dapat dilakukan dengan langkah-langkah sebagai berikut.

Buat data V sebagai fungsi a dan b. Data ini dapat dihasilkan dengan komputasi (misal dengan bahasa Fortran) dengan kode sebagai berikut.

program potential
    implicit none

    real :: a, b, V
    integer :: i, j

    real, parameter :: k = 9.0e+9
    real, parameter :: q = 1.602e-19

    open(1, file = "../potential.dat")

    a = 0.0
	
    do i = 1, 100
        a = a + 1.0
        
        b = 1.0

        do j = 1, 100
            b = b + 1.0

            V = k * q * ((a + b) / (2.0 * (a * b))) * 1.0e+3

            if (b > a) then
                write(1, *) a, b, V * 1.0e+6
            end if
            
        end do

        write(1, *) "", "", ""
        write(1, *) "", "", ""

    end do

end program potential
Jika kode di atas dieksekusi (dengan gfortran compiler) akan menghasilkan file potential.dat.

Berikutnya, plot data di atas (potential.dat) dengan menggunakan aplikasi gnuplot (http://gnuplot.info) dengan kode program sebagai berikut.
    set term gif animate delay 10 size 400,300
    set output "potential-3d.gif"
    set title "Grafik potensial di pusat cakram"
    set grid
    set xrange [0:100]
    set yrange [0:100]
    set zrange [0:1.1]
    set xyplane at 0
    set xtics
    set ytics
    set ztics
    set label "jari-jari dalam (mm)" at 20,-10,-0.3 rotate by -9
    set label "jari-jari luar (mm)" at 120,20,-0.1  rotate by 30
    set label "potensial (µV)" at -25,0,0.2 rotate by 90
    set pm3d
    do for [i=0:99] {
    	splot "potential.dat" index i notitle
    }
Selamat mencoba!