def __init__(self, params): Mesh.__init__(self, params) self.n_cell = params.get("n_cell") self.L = params.get("length") self.x_min = params.get("x_min") h = self.L / self.n_cell self.h = [h for e in xrange(self.n_cell)] self.x = np.zeros(self.n_cell + 1) self.x_max = self.x_min + self.L self.x[0] = self.x_min for e in xrange(self.n_cell): self.x[e + 1] = self.x[e] + self.h[e]