diff --git a/examples/injection_examples/basic_tutorial.py b/examples/injection_examples/basic_tutorial.py
index 90202b5f492755e2f4b23f18e0b7760e857e1807..1e90bcddcda848fb378792bb6ff7f788ba10a41f 100644
--- a/examples/injection_examples/basic_tutorial.py
+++ b/examples/injection_examples/basic_tutorial.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a reduced parameter space for an injected signal.
 
diff --git a/examples/injection_examples/binary_neutron_star_example.py b/examples/injection_examples/binary_neutron_star_example.py
index cec8fc5e098182d8da972a4b6adc0b1089a675ad..4f826a060037129c891286c28e2a022d626a6edd 100644
--- a/examples/injection_examples/binary_neutron_star_example.py
+++ b/examples/injection_examples/binary_neutron_star_example.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a binary neutron star
 system taking into account tidal deformabilities.
diff --git a/examples/injection_examples/calibration_example.py b/examples/injection_examples/calibration_example.py
index d111d3896bbdc4d2839e79b84c759a9f003254c8..9f2052ffb55d019aa1509ae251aa6f57d4bc9c14 100644
--- a/examples/injection_examples/calibration_example.py
+++ b/examples/injection_examples/calibration_example.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation with calibration
 uncertainties included.
diff --git a/examples/injection_examples/create_your_own_source_model.py b/examples/injection_examples/create_your_own_source_model.py
index f304b97abb0d3e418d8b93ff3a41b4d56e32f3a5..254fdfc52f4dd76aed5df74fbb9b31c79b134787 100644
--- a/examples/injection_examples/create_your_own_source_model.py
+++ b/examples/injection_examples/create_your_own_source_model.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 A script to demonstrate how to use your own source model
 """
diff --git a/examples/injection_examples/create_your_own_time_domain_source_model.py b/examples/injection_examples/create_your_own_time_domain_source_model.py
index ac45e1ef987ec27868d643ce590b74debd481c53..5ea57fb58897f5e0cace5f9a9ad22c40c0ba3afb 100644
--- a/examples/injection_examples/create_your_own_time_domain_source_model.py
+++ b/examples/injection_examples/create_your_own_time_domain_source_model.py
@@ -1,3 +1,4 @@
+#!/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
diff --git a/examples/injection_examples/eccentric_inspiral.py b/examples/injection_examples/eccentric_inspiral.py
index 235536ca48482e4bf39d6e58f7533f9822f22239..75ff93b572e85f9e2abfae62b341057d9c51b2f3 100644
--- a/examples/injection_examples/eccentric_inspiral.py
+++ b/examples/injection_examples/eccentric_inspiral.py
@@ -1,4 +1,4 @@
-#!/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
diff --git a/examples/injection_examples/how_to_specify_the_prior.py b/examples/injection_examples/how_to_specify_the_prior.py
index 706bdaedc1da97fde24db50b10ded0a56ec0c086..b92288d90237f57fceb84b48704ed296e67bf9bf 100644
--- a/examples/injection_examples/how_to_specify_the_prior.py
+++ b/examples/injection_examples/how_to_specify_the_prior.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate how to specify the prior distributions used for parameter estimation.
 """
diff --git a/examples/injection_examples/marginalized_likelihood.py b/examples/injection_examples/marginalized_likelihood.py
index 300649719a3519999ed78eab6256a18f5c98367f..83259d627e80909796d646be386ced0fbd2e2e99 100644
--- a/examples/injection_examples/marginalized_likelihood.py
+++ b/examples/injection_examples/marginalized_likelihood.py
@@ -1,4 +1,4 @@
-#!/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.
diff --git a/examples/injection_examples/plot_time_domain_data.py b/examples/injection_examples/plot_time_domain_data.py
index d03c3200a4f5303c3ecad0f045701dd0e1521efd..f3d056b7a9a594a328a200c84ac334ec50775093 100644
--- a/examples/injection_examples/plot_time_domain_data.py
+++ b/examples/injection_examples/plot_time_domain_data.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 """
 from __future__ import division, print_function
diff --git a/examples/injection_examples/sine_gaussian_example.py b/examples/injection_examples/sine_gaussian_example.py
index c63691af107ab8b9ac680fe034da0f30b4f5f581..4bbf7c6b0181c28212c26731acbff17539c29f28 100644
--- a/examples/injection_examples/sine_gaussian_example.py
+++ b/examples/injection_examples/sine_gaussian_example.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on a sine gaussian injected signal.
 
diff --git a/examples/injection_examples/using_gwin.py b/examples/injection_examples/using_gwin.py
index f6f0e28d640c45a9fbbc7fe9c47286e9568df0da..a213c490c6a30bd51dda582baa1675b6c300f3c3 100644
--- a/examples/injection_examples/using_gwin.py
+++ b/examples/injection_examples/using_gwin.py
@@ -1,4 +1,4 @@
-#!/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.
diff --git a/examples/open_data_examples/GW150914.py b/examples/open_data_examples/GW150914.py
index b0876d243a6634ff76368a72bdcf6cdcf465a625..aa5ca44af19f6575feb8ed7e8b4e373b001bec0a 100644
--- a/examples/open_data_examples/GW150914.py
+++ b/examples/open_data_examples/GW150914.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate running parameter estimation on GW150914 using open
 data.
diff --git a/examples/open_data_examples/GW150914_minimal.py b/examples/open_data_examples/GW150914_minimal.py
index 339eb3b13558c9b5f05e31916e150c6db6da3b30..520fe6ae9aa3c7428207e6880660fc62dbc9f827 100644
--- a/examples/open_data_examples/GW150914_minimal.py
+++ b/examples/open_data_examples/GW150914_minimal.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 Tutorial to demonstrate the minimum number of steps required to run parameter
 stimation on GW150914 using open data.
diff --git a/examples/other_examples/add_multiple_results.py b/examples/other_examples/add_multiple_results.py
index b4f99ba9669ae95ab7adb3372b007f41fd840086..f4c565555ad6250390622cd71fd1f433a7599c64 100644
--- a/examples/other_examples/add_multiple_results.py
+++ b/examples/other_examples/add_multiple_results.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 An example of running two sets of posterior sample estimations and adding them
 """
diff --git a/examples/other_examples/gaussian_example.py b/examples/other_examples/gaussian_example.py
index 318222e2d1fe62b003bc543ceb5c8c1aedc95d9d..bf7730286f12ddc70ac9bfe647e6e637fd4004ac 100644
--- a/examples/other_examples/gaussian_example.py
+++ b/examples/other_examples/gaussian_example.py
@@ -1,4 +1,4 @@
-#!/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
diff --git a/examples/other_examples/get_LOSC_event_data.py b/examples/other_examples/get_LOSC_event_data.py
index 7a8836ba110f881e687860a718b0c4092c734a35..6197426eb6e02869fb73e18b349d8dde92ca1fc5 100644
--- a/examples/other_examples/get_LOSC_event_data.py
+++ b/examples/other_examples/get_LOSC_event_data.py
@@ -1,3 +1,4 @@
+#!/usr/bin/env python
 """ Helper script to faciliate downloading data from LOSC
 
 Usage: To download the GW150914 data from https://losc.ligo.org/events/ 
diff --git a/examples/other_examples/linear_regression_pymc3.py b/examples/other_examples/linear_regression_pymc3.py
index 3328d8cf5722fcb5f3db5047f6d423ff645285f7..c10d5d3c93480254bc6ca2fa4d05ea128f0e2718 100644
--- a/examples/other_examples/linear_regression_pymc3.py
+++ b/examples/other_examples/linear_regression_pymc3.py
@@ -1,4 +1,4 @@
-#!/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
diff --git a/examples/other_examples/linear_regression_pymc3_custom_likelihood.py b/examples/other_examples/linear_regression_pymc3_custom_likelihood.py
index 2f36967947d1ff0d4087758bb3070b8aa37569f9..9321ea4c82621250201b3af8b97484845cd2e9a7 100644
--- a/examples/other_examples/linear_regression_pymc3_custom_likelihood.py
+++ b/examples/other_examples/linear_regression_pymc3_custom_likelihood.py
@@ -1,4 +1,4 @@
-#!/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
diff --git a/examples/other_examples/occam_factor_example.py b/examples/other_examples/occam_factor_example.py
index b164a3298b55fc54ff6225a19136480c63521b73..b977eaeadf53086326ca4047a8a9aa64bb607b0c 100644
--- a/examples/other_examples/occam_factor_example.py
+++ b/examples/other_examples/occam_factor_example.py
@@ -1,4 +1,4 @@
-#!/bin/python
+#!/usr/bin/env python
 """
 
 As part of the :code:`bilby.result.Result` object, we provide a method to