場の成分の計算

• ポテンシャルから場を計算する
  • こちら の2.7 にもポテンシャルがでてくる
• つまりは一階偏微分係数を求める

Nx<-10
Ny<-10

V<-matrix((1:(Nx*Ny))^2,Nx,Ny,byrow=T)
V_x<-matrix(0,Nx,Ny)

V_x[c(-1,-Nx),]<-(V[c(-1,-2),]-V[c(-(Nx-1),-Nx),])/2
V_x[1,]<-V[2,]-V[1,]
V_x[Nx,]<-V[Nx,]-V[Nx-1,]

persp(V,col=4)
persp(V_x,col=4)

# gx<-gxy<-1
# Ix<-matrix(0,Nx-1,Ny)
# Ix<-(-1)*gx*(V[-1,]-V[-Nx,])+(-1)*gxy*(V_x[-1,]+V[-Nx,])/2