THMMY.gr

% creates two matrices one for all the x-values on the grid, and one for
% all the y-values on the grid.
x1 = linspace(-10,10,20);
x2 = linspace(-5,5,20);
[x,y] = meshgrid(x1,x2);

odefun = @(t,x) [ οι διαφορικές σου];
t=0; % we want the derivatives at each point at t=0
for i = 1:numel(x)
  xprime = odefun(t,[x(i); y(i)]);
   u(i) = xprime(1);
   v(i) = xprime(2);
end
    
figure();
quiver(x,y,u,v,'r'); 
xlabel('x1')
ylabel('x2')
axis tight equal;