27 Desember 2008

Komputasi Saintifik

Kemampuan R untuk melakukan Komputasi Saintifik adalah sangat bagus. Berikut ini contoh peragaan grafik trigonometri.

> h <- 0.1 # beri nilai h sama dengan 0.1


> x <- seq( 0, 2*pi, by=h)# x merupakan sekuens dari 0 samapi 2 pi

> y <- sin(x)# merupakan vektor yang nilainya adalah sin dari vektor x

> x1 <- x + h

> y1 <- sin(x1)

> dy <-(y1 - y)/h # bila h mendekati 0 maka dy adalah turunan y

> plot( x, y, type="l" )

Hasilnya adalah:
Plot01


























Gambar 1: Plot x dan y = sin x


> plot( x, y, type="n", axes=FALSE ) # menyiapkan bidang plot

> lines(x, y) # menggambar plot x dan y

> axis(side=1, at=seq(0, 2*pi, by=pi/2),
+ labels= expression(0,paste(pi,"/2"),
pi, paste("3",pi,"/2"),
+ paste("2",pi)) )


> axis(side=2, at=seq(-1, 1, by=0.25), labels= seq(-1, 1,
+ by=0.25) )


> box()

> lines(x, dy, lty=2, col=2 )

> abline(h=0, lwd=2)

> legend(4, 1, legend=c("f(x) = sin(x)","f’(x) = cos(x)"),
+ lty=c(1,2), col=c(1,2))

Hasilnya adalah:

Plot02b


























Gambar 2: Plot x terhadap y = sin x dan dy/dx = cos x

Sekian dulu ya, nanti dilanjukan lagi dengan topik lain ......;-)

1 komentar:

Riyan adriyan mengatakan...

misi pak... saya dah tes pake template bapa dan hasilnya bisa...

contoh hasilnya akan seperti ini :
CONTOH


ada perubahan sedikit cara membuatnya..

tips.doc