Спасибо
myhand, это то что нужно.
У меня есть еще один вопрос по gnuplot. Я пытаюсь нарисовать линии электрического поля, делаю по аналогии с
http://gnuplot.sourceforge.net/demo/vector.html, но векторы рисуются не правильно, почему то их начало не перпендикулярно к эквипотенциалям.
В примере все нормально, почему есть разница мне не понятно.
Код:
# harmonic number
m=2.
r(x,y)=sqrt(x*x+y*y)
t(x,y)=atan2(y,x)
# potential
# r<1
v1(x,y)=r(x,y)**m/m*cos(m*t(x,y))
# r>1
v2(x,y)=r(x,y)**(-m)/m*cos(m*t(x,y))
# for all r
v(x,y)=(r(x,y) < 1 ? v1(x,y) : v2(x,y))
#electric field
# inside r<1
e1x(x,y)= -r(x,y)**(m-1)*cos((m-1)*t(x,y))
e1y(x,y)= r(x,y)**(m-1)*sin((m-1)*t(x,y))
# r>1
e2x(x,y)= -r(x,y)**(-m-1)*cos((m-1)*t(x,y))
e2y(x,y)= r(x,y)**(-m-1)*sin((m-1)*t(x,y))
# field abы val
enorm1(x,y)=sqrt(e1x(x,y)*e1x(x,y)+e1y(x,y)*e1y(x,y))
enorm2(x,y)=sqrt(e2x(x,y)*e2x(x,y)+e2y(x,y)*e2y(x,y))
dx1(x,y)=e1x(x,y)/enorm1(x,y)
dy1(x,y)=e1y(x,y)/enorm1(x,y)
dx2(x,y)=e2x(x,y)/enorm2(x,y)
dy2(x,y)=e2y(x,y)/enorm2(x,y)
# form for vec ploting x,y,dx,dy
dx(x,y)=(r(x,y) < 1 ? dx1(x,y) : dx2(x,y))
dy(x,y)=(r(x,y) < 1 ? dy1(x,y) : dy2(x,y))
xmin=-2.
xmax=2.
ymin=-2.
ymax=2.
unset autoscale
set xr [xmin:xmax]
set yr [ymin:ymax]
set isosam 31,31
set view map
unset surface
set contour base
set cntrparam order 5
set cntrparam levels auto 30
set cntrparam points 20
unset border
unset xtics
unset ytics
splot v(x,y) w l
set table "equipo2.tmp"
replot
unset table
set term fig
set output "m3.fig"
set size ratio -1
plot "equipo2.tmp" w l
set term x11
set isosam 31,31
set table "field2xy.tmp"
splot v(x,y) w l
unset table
unset autoscale
set xr [xmin:xmax]
set yr [ymin:ymax]
set isosam 31,31
set key under Left reverse
plot "field2xy.tmp" u 1:2:(dx($1,$2)):(dy($1,$2)) w vec, \
"equipo2.tmp" w l
pause -1
