diff --git a/.gitignore b/.gitignore
index cd041b2c6..812fd6b63 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,3 +1,4 @@
+*.pyd
 *.pyc
 *~
 *build*
diff --git a/.travis.yml b/.travis.yml
index f454589ca..f9aa44ec4 100644
--- a/.travis.yml
+++ b/.travis.yml
@@ -1,7 +1,6 @@
 branches:
   only:
     - master
-    - ft-mixture-of-experts
 
 sudo: required
 
@@ -15,6 +14,7 @@ language: generic
 
 env:
 - PY=2.7
+- PY=3.6
 
 addons:
   apt:
diff --git a/README.md b/README.md
index 9bb471cb2..e2eb77ed5 100644
--- a/README.md
+++ b/README.md
@@ -13,17 +13,7 @@ SMT is distributed under the New BSD license.
 Version 0.2
 
 # Required packages
-Numpy    >= Version 1.9.2
-
-Scikit-learn  >= Version 0.13.1
-
-pyDOE >= Version 0.3.8
-
-six
-
-Cython
-
-Scipy
+SMT depends on the following modules: numpy, scipy, sk-learn, pyDOE2 and Cython. 
 
 # Installation
 Clone the repository from github then run:
diff --git a/appveyor.yml b/appveyor.yml
index 482e1aece..552ed35f9 100644
--- a/appveyor.yml
+++ b/appveyor.yml
@@ -7,11 +7,12 @@ platform:
 
 environment:
   matrix:
+    - PYTHON: 3.6
     - PYTHON: 2.7
 
 init:
   - set CONDA=Miniconda
-  - if "%PYTHON%" == "3.4" (set CONDA=%CONDA%3)
+  - if "%PYTHON%" == "3.6" (set CONDA=%CONDA%3)
   - if "%PLATFORM%" == "x64" (set CONDA=%CONDA%-x64)
   - echo %CONDA%
   - set PATH=C:\%CONDA%;C:\%CONDA%\Scripts;%PATH%
@@ -24,7 +25,7 @@ init:
   - IF "%ARCH%"=="64" (call "C:\Program Files (x86)\Microsoft Visual Studio 14.0\VC\vcvarsall.bat" x86_amd64)
 
 install:
-  - pip install -i https://pypi.anaconda.org/carlkl/simple mingwpy
+  #pip install -i https://pypi.anaconda.org/carlkl/simple mingwpy
   #Grab testflo for running test, to ensure examples are run
   - pip install git+https://github.com/OpenMDAO/testflo.git
   - pip install git+https://github.com/hwangjt/sphinx_auto_embed.git
diff --git a/setup.py b/setup.py
index 9b7a94f58..a21c510b7 100644
--- a/setup.py
+++ b/setup.py
@@ -67,7 +67,7 @@
     ],
     install_requires=[
         'scikit-learn',
-        'pyDOE',
+        'pyDOE2',
         'matplotlib',
         'numpydoc',
         'six>=1.10',
diff --git a/smt/extensions/tests/test_moe.py b/smt/extensions/tests/test_moe.py
index 70469066c..a7324b161 100644
--- a/smt/extensions/tests/test_moe.py
+++ b/smt/extensions/tests/test_moe.py
@@ -204,7 +204,7 @@ def run_moe_example():
         x0 = np.linspace(xlim[0, 0], xlim[0, 1], 20)
         x1 = np.linspace(xlim[1, 0], xlim[1, 1], 20)
         xv, yv = np.meshgrid(x0, x1)
-        x = np.array(zip(xv.reshape((-1,)), yv.reshape((-1,))))
+        x = np.array(list(zip(xv.reshape((-1,)), yv.reshape((-1,)))))
         prob = moe._proba_cluster(x)
 
         plt.subplot(221, projection='3d')
diff --git a/smt/sampling_methods/lhs.py b/smt/sampling_methods/lhs.py
index a7a717fa1..c2544dcfe 100644
--- a/smt/sampling_methods/lhs.py
+++ b/smt/sampling_methods/lhs.py
@@ -3,10 +3,10 @@
 
 This package is distributed under New BSD license.
 
-LHS sampling; uses the pyDOE package.
+LHS sampling; uses the pyDOE2 package.
 """
 from __future__ import division
-import pyDOE
+from pyDOE2 import lhs
 from six.moves import range
 from scipy.spatial.distance import pdist, cdist
 import numpy as np
@@ -38,7 +38,7 @@ def _compute(self, n):
         xlimits = self.options['xlimits']
         nx = xlimits.shape[0]
         if self.options['criterion'] != 'ese':
-            return pyDOE.lhs(nx, samples=n, criterion=self.options['criterion'])
+            return lhs(nx, samples=n, criterion=self.options['criterion'])
         elif self.options['criterion'] == 'ese':
             return self._ese(nx,n)
 
diff --git a/smt/surrogate_models/krg_based.py b/smt/surrogate_models/krg_based.py
index eb93fc8f0..45e1f4482 100644
--- a/smt/surrogate_models/krg_based.py
+++ b/smt/surrogate_models/krg_based.py
@@ -193,7 +193,7 @@ def _reduced_likelihood_function(self, theta):
         try:            
             C = linalg.cholesky(R, lower=True)
         except (linalg.LinAlgError, ValueError) as e:
-            print "exception : ", e
+            print("exception : ", e)
             return reduced_likelihood_function_value, par
         
         # Get generalized least squares solution
diff --git a/smt/utils/kriging_utils.py b/smt/utils/kriging_utils.py
index 50d34339c..96ac924d6 100644
--- a/smt/utils/kriging_utils.py
+++ b/smt/utils/kriging_utils.py
@@ -6,7 +6,7 @@
 
 import numpy as np
 from sklearn.cross_decomposition.pls_ import PLSRegression as pls
-from pyDOE import *
+from pyDOE2 import bbdesign
 
 def standardization(X,y,copy=False):