function [ x y ] = randomwalk2d( xinit, yinit, n )
%RANDOMWALK2D Summary of this function goes here
%   Detailed explanation goes here

    x = zeros(1,n+1);
    y = zeros(1,n+1);
    
    x(1) = xinit;
    y(1) = yinit;
    
    for j=2:n+1
        r = rand;
        if r <= 0.25
            x(j) = x(j-1);
            y(j) = y(j-1) + 1;
        elseif r <= 0.50
            x(j) = x(j-1) + 1;
            y(j) = y(j-1);
        elseif r <= 0.75
            x(j) = x(j-1);
            y(j) = y(j-1) - 1;
        else
            x(j) = x(j-1) - 1;
            y(j) = y(j-1);
        end            
    end

end