Hubble Profile Script

Here’s a useful little script I wrote in IDL. It allows you to send it 2 arrays (x and y data) and then you can plot the data, and fit it with a Hubble profile:

y = I * (1 + x/r_c)^(-a)

where I, r_c and a are constants representing the vertical scaling, the core radius and the exponent of decline, respectively.

It’s interactive, and allows you to adjust parameters and re-fit on the fly.

Click here to download it.

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

w

Connecting to %s