Tag: สมการ

เอาสักดอกไหม

จาก http://www.sakngoi.com/2015/10/30/สมการดอกกุหลาบ/ library(misc3d) rose <- function(x,theta,axis) { phi <- (pi/2)*exp(-theta/(8*pi)) X <- 1-(1/2)*((5/4)*(1- ((3.6 * theta) %% (2*pi))/pi)^2-1/4)^2 y <- 1.95653* x*x * (1.27689 *x-1)^2 * sin(phi) r <- X*(x *sin(phi)+y*cos(phi)) out<-switch(axis, x={ r*sin(theta) }, y={ r*cos(theta) }, z={ (X*(x* cos(phi)-y* sin(phi))) } ) return(out) } parametric3d( fx = function(u,v) rose(u,v,axis="x"), fy = function(u,v) rose(u,v,axis="y"), fz

maemod package

จากที่เคยเขียนไว้เกี่ยวกับ r package อันหนึ่งที่ผมเขียนชื่อ maemod สำหรับช่วยให้คนที่สนใจอยากคำนวณพวก ode ได้ง่ายขึ้น (ดูข้างล่าง) maemod (แม่มด) มีคนสนใจว่าถ้าพวกตัวแปร state ต่างๆนั้นเป็นแบบ array จะทำอย่างไร ผมเลยเขียนตัวอย่างพร้อมกับเพิ่มความสามารถด้านarrayนี้เข้าไป พอใช้ได้ไปคร่าวๆก่อน ดูตัวอย่างข้างล่างนี้ครับ # Example from Berkeley Madonna # for using Array # METHOD RK4 # # STARTTIME = 0 # STOPTIME = 20 # DT = 0.02 # # d/dt (A[1]) = -k[1]*A[1] # d/dt (A[2..n-1]) = k[i-1]*A[i-1]-k[i]*A[i]

HappyValentinesDay

R code สำหรับวันแห่งความรัก 🙂 t<-seq(-pi, pi,0.05) colors<-rainbow(length(t)) x<-16*sin(t)^3 y<-13*cos(t)-5*cos(2*t)-2*cos(3*t)-1*cos(4*t) p<-c(“H”, “a”, “p”,”p”,”y”,”V”,”a”,”l”,”e”,”n”,”t”,”i”,”n”,”e”,”s”,”D”,”a”,”y”) plot(x,y, type=”p”,pch=p,col=colors,xlab=”X”,ylab=”Y”)