Make all examples have a sensible shebang

examples/injection_examples/basic_tutorial.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a reduced parameter space for an injected signal.
 

examples/injection_examples/binary_neutron_star_example.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a binary neutron star
 system taking into account tidal deformabilities.

examples/injection_examples/calibration_example.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation with calibration
 uncertainties included.

examples/injection_examples/create_your_own_source_model.py

-#!/bin/python
+#!/usr/bin/env python
 """
 A script to demonstrate how to use your own source model
 """

examples/injection_examples/create_your_own_time_domain_source_model.py

+#!/usr/bin/env python
 """
 A script to show how to create your own time domain source model.
 A simple damped Gaussian signal is defined in the time domain, injected into noise in

examples/injection_examples/eccentric_inspiral.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a reduced parameter space
 for an injected eccentric binary black hole signal with masses & distnace similar

examples/injection_examples/how_to_specify_the_prior.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate how to specify the prior distributions used for parameter estimation.
 """

examples/injection_examples/marginalized_likelihood.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate how to improve the speed and efficiency of parameter estimation on an injected signal using
 phase and distance marginalisation.

examples/injection_examples/plot_time_domain_data.py

-#!/bin/python
+#!/usr/bin/env python
 """
 """
 from __future__ import division, print_function

examples/injection_examples/sine_gaussian_example.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a sine gaussian injected signal.
 

examples/injection_examples/using_gwin.py

-#!/bin/python
+#!/usr/bin/env python
 """
 An example of how to use bilby with `gwin` (https://github.com/gwastro/gwin) to
 perform CBC parameter estimation.

examples/open_data_examples/GW150914.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on GW150914 using open
 data.

examples/open_data_examples/GW150914_minimal.py

-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate the minimum number of steps required to run parameter
 stimation on GW150914 using open data.

examples/other_examples/add_multiple_results.py

-#!/bin/python
+#!/usr/bin/env python
 """
 An example of running two sets of posterior sample estimations and adding them
 """

examples/other_examples/gaussian_example.py

-#!/bin/python
+#!/usr/bin/env python
 """
 An example of how to use bilby to perform paramater estimation for
 non-gravitational wave data consisting of a Gaussian with a mean and variance

examples/other_examples/get_LOSC_event_data.py

+#!/usr/bin/env python
 """ Helper script to faciliate downloading data from LOSC
 
 Usage: To download the GW150914 data from https://losc.ligo.org/events/ 

examples/other_examples/linear_regression_pymc3.py

-#!/bin/python
+#!/usr/bin/env python
 """
 An example of how to use bilby to perform paramater estimation for
 non-gravitational wave data. In this case, fitting a linear function to

examples/other_examples/linear_regression_pymc3_custom_likelihood.py

-#!/bin/python
+#!/usr/bin/env python
 """
 An example of how to use bilby to perform paramater estimation for
 non-gravitational wave data. In this case, fitting a linear function to

examples/other_examples/occam_factor_example.py

-#!/bin/python
+#!/usr/bin/env python
 """
 
 As part of the :code:`bilby.result.Result` object, we provide a method to