Here's a way to do this in python with cruzdb using the TRIM13 from hg19 gene as an example:

from cruzdb import Genome
g = Genome('hg19')

for transcript in g.refGene.filter_by(name2='TRIM13'):
    print transcript.exons

the output is a list of start, end tuples for each transcript.

Similar questions like these pop up quite often here, so please search the archives.

That having been said:

• Genomic Cordinates From Ucsc to one such question that you can adapt to get what you want.
• If you are familiar with R/Bioconductor, you can easily address such questions with facilities provided by the GenomicFeatures, and you can read through its vignette to become more familiar with it.