time_of_transfer.m

function time = time_of_transfer(R1, R2, mu)
    %% Calculate the time of interplanetary transfer
    %
    % Jeremy Penn
    % 06/11/17
    %
    % Input:    o R1    - The orbital semimajor axis of the departing planet
    %           o R2    - The orbital semimajor axis of the arrival planet
    %           o mu    - grav parameter [OPTIONAL]
    %
    % Output:   o time - transfer flight time [s]
    %
    clc;
    
    %% constants
    if nargin == 2
        mu = 132.71e9; % [km^3/s^2]
    end
    
    %% calculate the transfer time
    time = (pi/sqrt(mu)) * ( (R1 + R2)/2 )^1.5; %[s]
    
end