Computing coefficients of modular forms