SC2Mapster

% Program Desc.: The program will create a virtual simulation of cue balls on a pool table.

clear all; clf;                     % Resets parameters

% [Region] Define global variables
table = [2.54, 1.27];               % Size of virtual table [meters]
dt = 0.01;                          % Time interval between calculations [seconds]
time = 10;                          % Time to run simulation [seconds]
% [Region] Define global variables

% [Region] Initial ball parameters
balls(1).pos = [0.4, 0.6];          % Ball no. 1 initial transform from bottom-left [meters]
balls(1).v = [0, 0];                % Ball no. 1 initital velocity [m/s]
balls(1).r = 0.05;                  % Ball no. 1 initial radius for impact calculations [meters]
balls(1).color = 'green';           % Ball no. 1 initial color.

balls(2).pos = [table(1)-0.4 0.6];  % Ball no. 2 initial transform from bottom-left [meters]
balls(2).v = [0, 0];                % Ball no. 2 initital velocity [m/s]
balls(2).r = 0.05;                  % Ball no. 2 initial radius for impact calculations [meters]
balls(2).color = 'red';             % Ball no. 2 initial color.

balls(3).pos = [0.25, 1];           % Ball no. 3 initial transform from bottom-left [meters]
balls(3).v = [6, -1.8];             % Ball no. 3 initital velocity [m/s]
balls(3).r = 0.05;                  % Ball no. 3 initial radius for impact calculations [meters]
balls(3).color = 'blue';            % Ball no. 3 initial color.
% [Region] Initial ball parameters



for t = 0:dt:time
    for a = 1:length(balls)
       
        % [Region] Update position
        % balls(a).pos = balls(a).pos + balls(a).v .* dt;     % Update transform from old transform + current speed * passed time interval.
        balls(a).pos = balls(a).pos + balls(a).v .* dt;

NEW CODE HERE?????

        % [Region] Update position
        
        
        % [Region] Check for collision with wall
        if balls(a).pos(1)-balls(a).r <= 0 ...              % Compare (1) coordinate of transform to radius of current transform.
                || balls(a).pos(1)+balls(a).r >= table(1)   % Compare (1) coordinate of transform to table variables (1) coordinate.
            balls(a).v(1) = -balls(a).v(1);                 % If True: Change velocity to reverse direction for (1) coordinate.
        end
        if balls(a).pos(2)-balls(a).r <= 0 ...              % Compare (2) coordinate of transform to radius of current transform.
                || balls(a).pos(2)+balls(a).r >= table(2)   % Compare (2) coordinate of transform to table variables (2) coordinate.
            balls(a).v(2) = -balls(a).v(2);                 % If True: Change velocity to reverse direction for (2) coordinate.
        end
        % [Region] Check for collision with wall
        
        % [Region] Check for collision with balls
        for b = a+1:length(balls)
            if norm(balls(b).pos - balls(a).pos) <= balls(b).r + balls(a).r
 
                % Calculate unit normal
                un = (balls(b).pos - balls(a).pos) ./ norm(balls(b).pos - balls(a).pos);
                ut = [-un(2), un(1)]; % Unit tangent
                % Project velocity vectors into (normal,tangent)-coordinats
                v_an = dot(un,balls(a).v);
                v_at = dot(ut,balls(a).v);
                v_bn = dot(un,balls(b).v);
                v_bt = dot(ut,balls(b).v);
                % Find velocities after collision
                v_anp = v_bn.*un;
                v_atp = v_at.*ut;
                v_bnp = v_an.*un;
                v_btp = v_bt.*ut;
                % Add up the components for the resulting velocity
                balls(a).v = v_anp + v_atp;
                balls(b).v = v_bnp + v_btp;
            end
        end
        % [Region] Check for collision with balls
    end
    
    show(table,balls)
    pause(dt)
end