We exhibit a method to numerically compute power series expansions of modular forms on a cocompact Fuchsian group, using the explicit computation a fundamental domain and linear algebra. As applications, we compute Shimura curve parametrizations of elliptic curves over a totally real field, including the image of CM points, and equations for Shimura curves.