#!/usr/bin/gnuplot # concordia_u_pb.plt # Krivka zvana concordia pro radiometrickou metodu U-Pb. # Miroslav Broz (miroslav.broz@email.cz), May 31st 2011 yr = 365.25*86400. Gyr = 1.e9*yr T12_238 = 4.468e9*yr T12_235 = 7.038e8*yr lambda_238 = log(2)/T12_238 lambda_235 = log(2)/T12_235 print "lambda_238 = ", lambda_238, " s^-1" print "lambda_235 = ", lambda_235, " s^-1" ######################################################################## set tr [0:4.5*Gyr] set parametric set samples 10 x(t) = exp(lambda_235*t)-1 y(t) = exp(lambda_238*t)-1 set term table set out "concordia_u_pb.dat" p x(t),y(t) ######################################################################## set term x11 set out set xl "207 Pb/235 U" set yl "206 Pb/238 U" set samples 100 set style line 1 lt 1 set style line 2 lt 2 p x(t),y(t) w l ls 1,\ "concordia_u_pb.dat" w p ls 2 pa -1 set term png small set out "concordia_u_pb.png" rep