diff --git a/prism/src/explicit/ConstructModel.java b/prism/src/explicit/ConstructModel.java
index 5af18f9e..8f720e9f 100644
--- a/prism/src/explicit/ConstructModel.java
+++ b/prism/src/explicit/ConstructModel.java
@@ -35,6 +35,7 @@ import java.util.List;
 import parser.State;
 import parser.Values;
 import parser.VarList;
+import parser.ast.Expression;
 import parser.ast.ModulesFile;
 import prism.ModelType;
 import prism.Prism;
@@ -137,11 +138,6 @@ public class ConstructModel extends PrismComponent
 		int i, j, nc, nt, src, dest;
 		long timer;
 
-		// Don't support multiple initial states
-		if (modulesFile.getInitialStates() != null) {
-			throw new PrismException("Cannot do explicit-state reachability if there are multiple initial states");
-		}
-
 		// Display a warning if there are unbounded vars
 		VarList varList = modulesFile.createVarList();
 		if (varList.containsUnboundedVariables())
@@ -180,16 +176,30 @@ public class ConstructModel extends PrismComponent
 		// Initialise states storage
 		states = new IndexedSet<State>(true);
 		explore = new LinkedList<State>();
-		// Add initial state to lists/model
-		if (modulesFile.getInitialStates() != null) {
-			throw new PrismException("Explicit model construction does not support multiple initial states");
+		// Add initial state(s) to 'explore'
+		// Easy (normal) case: just one initial state
+		if (modulesFile.getInitialStates() == null) {
+			state = modulesFile.getDefaultInitialState();
+			explore.add(state);
 		}
-		state = modulesFile.getDefaultInitialState();
-		states.add(state);
-		explore.add(state);
-		if (!justReach) {
-			modelSimple.addState();
-			modelSimple.addInitialState(0);
+		// Otherwise, there may be multiple initial states
+		// For now, we handle this is in a very inefficient way
+		else {
+			Expression init = modulesFile.getInitialStates();
+			List<State> allPossStates = varList.getAllStates();
+			for (State possState : allPossStates) {
+				if (init.evaluateBoolean(modulesFile.getConstantValues(), possState)) {
+					explore.add(possState);
+				}
+			}
+		}
+		// Copy initial state(s) to 'states' and to the model
+		for (State initState : explore) {
+			states.add(initState);
+			if (!justReach) {
+				modelSimple.addState();
+				modelSimple.addInitialState(modelSimple.getNumStates() - 1);
+			}
 		}
 		// Explore...
 		src = -1;
diff --git a/prism/src/parser/VarList.java b/prism/src/parser/VarList.java
index cd64b3c7..08d32e91 100644
--- a/prism/src/parser/VarList.java
+++ b/prism/src/parser/VarList.java
@@ -420,6 +420,48 @@ public class VarList
 		return allValues;
 	}
 
+	/**
+	 * Get a list of all possible states over the variables in this list. Use with care!
+	 */
+	public List<State> getAllStates() throws PrismLangException
+	{
+		List<State> allStates;
+		State state, stateNew;
+
+		int numVars = getNumVars();
+		allStates = new ArrayList<State>();
+		allStates.add(new State(numVars));
+		for (int i = 0; i < numVars; i++) {
+			if (getType(i) instanceof TypeBool) {
+				int n = allStates.size();
+				for (int j = 0; j < n; j++) {
+					state = allStates.get(j);
+					stateNew = new State(state);
+					stateNew.setValue(i, true);
+					state.setValue(i, false);
+					allStates.add(stateNew);
+				}
+			} else if (getType(i) instanceof TypeInt) {
+				int lo = getLow(i);
+				int hi = getHigh(i);
+				int n = allStates.size();
+				for (int j = 0; j < n; j++) {
+					state = allStates.get(j);
+					for (int k = lo + 1; k < hi + 1; k++) {
+						stateNew = new State(state);
+						stateNew.setValue(i, k);
+						allStates.add(stateNew);
+					}
+					state.setValue(i, lo);
+				}
+			} else {
+				throw new PrismLangException("Cannot determine all values for a variable of type " + getType(i));
+			}
+		}
+
+		return allStates;
+	}
+
 	/**
 	 * Convert a bit vector representing a single state to a State object. 
 	 */