こちらで、２次元平面における拡散が、ベクトル場で起きている ここで言うベクトル場では、以下で周回条件が入っている f<-function(x,y,Ox,Oy,phi){ #r<-sqrt((x-Ox)^2+(y-Oy)^2) r<-1/3 th<-atan2(y-Oy,x-Ox) v<-c(-r*sin(th+phi),r*cos(th+phi)) v } for(…

こちらで、２次元正円を２次元周回曲線に対応付ける関数について考えた 対応関数は解けないことも多いので、RのrootSolveパッケージを使うことも書いた ロトカ=ヴォルテラのx,y対称な曲線を描いてみる( うまくいっているように見える(こちらで紹介したロトカ…