[Relativistic velocity addition] #demo #sample Matlab #script file used in #astrophysics #coursework

rvel.m

% Script : rvel.m
% Graph demonstrating the velocity addition formula 
% Easily convertible to a Live Script

% Create basic plot
beta=0:0.02:1;
rv1=2*beta;
rv2=rv1./(1+beta.*beta);
rvel_limit=ones(size(beta));
Traces=plot(beta,rv1,beta,rv2,beta,rvel_limit);

% Annotate the plot
xlabel('\beta');
ylabel('\beta_{relative}');
legend('Classical','STR');
title('Velocity Addition Formula');


% Modify the Traces
set(Traces(1),'linewidth',2)
set(Traces(2),'linewidth',2)
set(Traces(3),'linestyle','--')

% Set button for printing in color
uicontrol('string','Print in Color','pos',[10 20 100 20], ...
          'call','print -dwinc');