from visual.graph import *

## Heat capacity data for Problem 10.4

nA = 6.02e23    ## Avogadro's number

## This converts heat capacity data from
## a per mole to a per atom basis and
## plots C vs. T (as dots) for Al and for Pb.
## Data are in the form (T,C).

Aldata =[(20,0.23/nA),(40,2.09/nA), (60,5.77/nA), (80,9.65/nA), (100,13.04/nA),
         (150,18.52/nA), (200, 21.58/nA), (250,23.25/nA), (300, 24.32/nA),
         (400,25.61/nA)]
Pbdata = [(20,11.01/nA),(40,19.57/nA), (60,22.43/nA), (80,23.69/nA), (100,24.43/nA),
         (150,25.27/nA), (200, 25.87/nA), (250,26.36/nA), (300, 26.82/nA),
         (400,27.45/nA)]

Alg = gdots(color=color.magenta,pos=Aldata)
Pbg = gdots(color=color.cyan, pos=Pbdata)