A solver for 1D systems of parabolic/elliptic PDE's in Julia