Many biomolecular condensates act as viscoelastic complex fluids with distinct cellular functions. Deciphering the viscoelastic behavior of biomolecular condensates can provide insights into their spatiotemporal organization and physiological roles within cells. Although there is significant interest in defining the role of condensate dynamics and rheology in physiological functions, the quantification of their time-dependent viscoelastic properties is limited and is mostly done through experimental rheological methods. Here, we demonstrate that a computational passive probe microrheology technique, coupled with continuum mechanics, can accurately characterize the linear viscoelasticity of condensates formed by intrinsically disordered proteins (IDPs). Using a transferable coarsegrained protein model, we first provide a physical basis for choosing optimal values that define the attributes of the probe particle, namely, its size and interaction strength with the residues in an IDP chain. We show that the technique captures the sequencedependent viscoelasticity of heteropolymeric IDPs that differ in either sequence charge patterning or sequence hydrophobicity. We also illustrate the technique's potential in quantifying the spatial dependence of viscoelasticity in heterogeneous IDP condensates. The computational microrheology technique has important implications for investigating the time-dependent rheology of complex biomolecular architectures, resulting in the sequence−rheology−function relationship for condensates.