function [X,Y,F]=ros_plot()

xx=-2:0.01:2;
yy=-2:0.01:3;

nx=length(xx);
ny=length(yy);

X=zeros(nx,ny);
Y=zeros(nx,ny);
F=zeros(nx,ny);

for i=1:nx
    for j=1:ny
        x_val=xx(i);
        y_val=yy(j);
        [f_val]=fun([x_val;y_val]);
        
        F(i,j)=f_val;
        X(i,j)=x_val;
        Y(i,j)=y_val;
    end
end