3D版邪恶的曲线

我的公众号biobabble第一篇原创文章《邪恶的曲线》,因为太邪恶,当时只推送给男性关注者,由于不是推送给所有人,所以历史记录是翻不到的,之前有再文字推送一下链接,结果评论区全是差评,我也是蛮囧rz的…

无意中我又翻到了一张2012年2月份画的3D版本,尽管你们要差评,还是想分享一下。代码就在截屏之中,想玩的话,可以照着输入,支持鼠标旋转哦,可以慢慢把玩(太邪恶了)。

Continue reading

今天在微博上看到这坑爹的方程: $ (x^2+y^2-1)^3 = {x^2} {y^3}$ 画出来如下: 跟个屁股似的,sigh… 翻出azalea的老文:http://azaleasays.com/2008/06/18/fomula-of-love/ 用ggplot2画一下这个爱的方程:$17x^2-16|x|y+17y^2 = 225 $ x <- seq(-sqrt(17), sqrt(17), 0.001) y1 <- 8*abs(x)/17 + 15* sqrt(17 - x^2)/17 y2 <- 8*abs(x)/17 - 15* sqrt(17 - x^2)/17 d <- data.frame(x=c(x,x),y=c(y1,y2)) require(ggplot2) p <- ggplot(d, aes(x,y)) p + geom_point(color="red") 还可以画出红心,适合今天七夕的日子。 p + geom_line(color="red") 画函数图,这种事情,还是用CAS方便点。 通过maxima来画,一条指令就行,还不用解方程: contour_plot(17*x^2-16*abs(x)*y+17*y^2-225, [x,-sqrt(17),sqrt(17)],[y,-15*sqrt(17)/17, sqrt(17)])

Continue reading

下载imaxima   

http://members3.jcom.home.ne.jp/imaxima/Site/Download%20and%20Install_files/imaxima-imath-0.97a.tar.gz

./configure
make
sudo make install

进行安装。 默认将安装在/usr/local/share/emacs/site-lisp下,在~/.emacs下加入

(load "imaxima")

启动emacs,按M-x,键入imaxima运行maxima 但是进行运算时出错。

需要breqn094包,可从ftp://ftp.ams.org/pub/tex/breqn.zip下载,

cd /usr/share/texmf/tex/latex/
sudo unzip ~/Downloads/breqn.zip
sudo texhash

Continue reading

Author's picture

Guangchuang Yu

a senior-in-age-but-not-senior-in-knowledge bioinformatician

Postdoc researcher

Hong Kong