2020/03/30

トラクトリックス(tractrix)の書き方

トラクトリックスとは


出典: フリー百科事典『ウィキペディア(Wikipedia)

 トラクトリックス (tractrix) は直交座標の方程式 (1) によって表される曲線である。


数式表現


x=loga±a2y2ya2y2(1)=acosh1(ay)a2y2

書き方


 R言語での書き方は以下のとおり。

 条件: a=1 の時

1
2
3
4
5
6
7
8
9
10
11
y=seq(0, 1, length=300)
a=1
x1=a*log((a+sqrt(a^2-y^2))/y)-sqrt(a^2-y^2)
x2=a*log((a-sqrt(a^2-y^2))/y)+sqrt(a^2-y^2)
svg("tractrix.svg", width=5, height=5)
plot(x1, y, type="l", asp=1, lwd=2, col=2, xlim=c(-2, 2), ylim=c(-2, 2), ann=F)
par(new=T)
plot(x2, y, type="l", asp=1, lwd=2, col=2, xlim=c(-2, 2), ylim=c(-2, 2), xlab="", ylab="")
abline(h=0)
abline(v=0)
dev.off()

0 件のコメント:

コメントを投稿

マンデルブロ集合の彩色方法(5)

06.knit 1 発散判定式を変更する mandelbrot() 内の発散判定式 |zn|>2 を変更する...