classdef beads2d
    
    properties
    end
    
    methods
        % Samples the vector field of the Beads problem
        function [u, v] = sample(obj, x, y, t)
            u = -2 * y + (2 * sin(t)) / 3;
            v = +2 * x - (2 * cos(t)) / 3;
        end
        
        % Samples the x-partial of the vector field of the Beads problem
        function [u, v] = sample_dx(obj, x, y, t)
            u = 0;
            v = 2;
        end
        
        % Samples the y-partial of the vector field of the Beads problem
        function [u, v] = sample_dy(obj, x, y, t)
            u = -2;
            v = 0;
        end
        
        % Samples the t-partial of the vector field of the Beads problem
        function [u, v] = sample_dt(obj, x, y, t)
            u = (2 * cos(t)) / 3;
            v = (2 * sin(t)) / 3;
        end
    end
    
end

