reimink function

Calculate upper and lower concordia intercepts from discordant detrital zircon data