Как построить гиперболическое векторное поле в R?

Question

Я пытаюсь построить векторное поле следующим образом:

http://www.aanda.org/articles/aa/full/2006/11/aa4228-05/img58.gif

который имел бы уравнение вида (y2 /b)-(x2 /k)=1,

в сюжете что-то вроде этого ( http://math.haifa.ac.il/ROVENSKI/rovenski/images/ROV05_628.gif)

с милыми стрелочками и тому подобным.

Кто-нибудь знает, как это можно сделать? Спасибо!

Answer 1

df <- expand.grid(x=seq(-2, 2, .1), y=seq(-2, 2, .1)); k=2; b=4
df$z <- with(df, (y^2)-(x^2) )
library(raster)
library(rasterVis)
r <- rasterFromXYZ(df)
projection(r) <- CRS("+proj=longlat +datum=WGS84")
vectorplot(r, par.settings=RdBuTheme)

Ypu может играть с размерами стрелок. Эта настройка выглядела немного более приятной на моем экране: lwd.arrows=0.3,length=unit(.04, 'npc'))