McGill University
School of Computer Science

Computer Science COMP 199 (Winter term)
Excursions in Computing Science

%MATLABpak03: StanTrav

% StanTrav.m	THM	070129
% Stanley standing and Travis travelling for Galileo and Lorentz
hold off
clear Ug, clear Vg, clear Ul, clear Vl
v = 0.1;
G = [1,0;-v,1];
L = (1/sqrt(1-v^2))*[1,-v;-v,1];
T = [0,2,0,-2,-1.5,-.5];
X = [0,0.2,0.8,4/3,0,0];
for j = 1:6
    TX = [T(j);X(j)];            %TX is Stan's frame
    UV = G*TX;                   %UV is Trav's frame: this is Galileo
    Ug(j) = UV(1);
    Vg(j) = UV(2);
    UV = L*TX;                   %UV is Trav's frame: this is Lorentz
    Ul(j) = UV(1);
    Vl(j) = UV(2);
end
plot(X,Y,'o')
hold on
plot(Ug,Vg,'+')                   %put this in for Galileo
%plot(Ul,Vl,'+')                   %put this in for Lorentz
xlabel('hours')
ylabel('light-hours')
axis([-3 3 -1 2])
legend('in Stan''s frame','in Trav''s frame')