%%
% Divergence
% https://jp.mathworks.com/help/matlab/ref/divergence.html
clear; clc;     % 変数をクリア、コマンドウィンドウの表示を消去
[x,y] = meshgrid(-3:0.3:3,-3:0.3:3);    % 2変数グラフ描画用の変数を準備
Fx = x.^2;
Fy = y;
figure;                     % 新たにグラフ描画ウィンドウを作成
quiver(x,y,Fx,Fy);          % ベクトル場の成分(Fx,Fy)をプロット
%%
D = divergence(x,y,Fx,Fy);
hold on;
contour(x,y,D,'ShowText','on');
%%
hold off;
figure;
surf(x,y,D);