Scilab champ de vecteur 2D

fx(i,j) is the x component of the vector field at point (x(i),y(j)).

// using meshgrid to get fx and fy
clf()
x = linspace(-1,1,30);
y = linspace(-1,1,30);
[X,Y] = meshgrid(x,y);
fy = 2.*Y;
fx = 2.*X;
champ(x,y,fx',fy'); //here we have transposed fx and fy

\documentclass[12pt,tikz]{standalone}
\usepackage{pgfplots}
\begin{document}
\begin{tikzpicture}
  \begin{axis}[title=Quiver and plot table]
    \addplot 
      [
        blue,
        quiver={u=\thisrow{u},v=\thisrow{v}},
        -stealth,
      ] 
      table 
      {
        x y  u v
        1 0  0 1
        2 2 -2 2
        3 0  0 3
        0 1 -1 0
      };
  \end{axis}
\end{tikzpicture}
\end{document}