DuckieDai Clausius-Clapeyron Equation Python Program

← Back to Calculator Program Library

Independent compatibility resource: Hidden Word Scanner is independent and is not affiliated with or endorsed by Texas Instruments. TI product names describe compatibility only.

Action-count privacy: Successful downloads, preview starts, and preview completions increment same-origin aggregate counters. No entered values, source text, IP address, account, cookie, device identifier, or fingerprint is stored with these counts.

chemistry · MIT License

Overview

Use the integrated Clausius-Clapeyron relation to solve either pressure, either temperature, or enthalpy of vaporization.

Usage

Inputs

Two vapor pressures in the same unit
Absolute temperatures in kelvin
Enthalpy of vaporization in kJ/mol when supplied

Outputs

Requested pressure, temperature, or enthalpy of vaporization

Units: Pressure units may vary but must match; temperature is kelvin; enthalpy is kJ/mol; R is 8.314 J/(mol*K).

Assumptions and limitations

Assumptions

Delta Hvap is constant over the temperature interval.
The vapor behaves ideally and liquid molar volume is negligible relative to vapor.

Constraints

Pressures, temperatures, and supplied Delta Hvap must be positive; paired temperatures must differ when solving Delta Hvap.

Known failures

Wide temperature ranges, near-critical conditions, or phase changes can violate the approximation.
Extreme pressure ratios may overflow or underflow floating-point display.

Compatibility and review

Review status: Not calculator tested
Tested on: The owner reports personally testing the supplied archive programs on a calculator; exact model, OS, and Python App versions were not recorded
What happened: Owner reported interactive testing is disclosed, but the published bytes are not labeled independently calculator tested because the environment is unrecorded and Arrhenius received a narrow intake correction
Dependencies: math
Suggested calculator name: CCEQ — you can give it another valid, unique name when you transfer it
Desktop test cases: Not available for this interactive-only source

This exact source auto-launches an interactive calculator session, so compatibility evidence comes from the recorded physical-device review rather than an importable desktop fixture.

Browser preview

Try the DuckieDai calculator screen

Follow the same opening, solve menu, prompts, validation, and result flow as the downloadable program.

Simulation boundary: This preview uses a reviewed browser adapter. It does not execute the downloaded Python and does not prove TI-Python or calculator compatibility.

JavaScript is required for the optional browser simulation. Source viewing and downloading remain available without it.

Exact reviewed bytes

Source code

Download .py

View source below, copy it when JavaScript is available, or download the exact .py bytes.

from math import log,exp

R=8.314

def duckiedai_intro(program_name,wait=True):
 title=program_name[:16]
 empty=16-len(title)
 left=empty//2
 right=empty-left
 print("+----------------------+")
 print("|   "+" "*left+title+" "*right+"   |")
 print("|                      |")
 print("|       __             |")
 print("|   ___(o )>  quack!   |")
 print("|   \\ <_. )           |")
 print("|    `---'             |")
 print("|     by DuckieDai     |")
 print("+----------------------+")
 print("Loading...")
 if wait:
  input("Press enter to start ")
  print("\n"*6)

def pos(s):
 while 1:
  try:
   x=float(input(s))
   if x>0:return x
  except:pass
  print("Enter a positive #")

def menu():
 while 1:
  try:
   x=int(input("Choose 1-5: "))
   if x>0 and x<6:return x
  except:pass
  print("Choose 1 to 5")

duckiedai_intro("CLAUSIUS CLAP")
print("CLAUSIUS-CLAPEYRON")
print("1 Find P2")
print("2 Find T2")
print("3 Find Hvap")
print("4 Find P1")
print("5 Find T1")
c=menu()

if c==1:
 p1=pos("P1: ")
 t1=pos("T1 K: ")
 t2=pos("T2 K: ")
 h=pos("Hvap kJ/mol: ")*1000
 try:
  p2=p1*exp((-h/R)*(1/t2-1/t1))
  print("P2 =",p2)
  print("P2/P1 =",p2/p1)
 except:print("P2 too large")

elif c==2:
 p1=pos("P1: ")
 p2=pos("P2: ")
 t1=pos("T1 K: ")
 h=pos("Hvap kJ/mol: ")*1000
 x=1/t1-(R/h)*log(p2/p1)
 if x<=0:print("T2 is not positive")
 else:
  t2=1/x
  print("T2 K =",t2)
  print("T2 C =",t2-273.15)

elif c==3:
 p1=pos("P1: ")
 p2=pos("P2: ")
 t1=pos("T1 K: ")
 t2=pos("T2 K: ")
 x=1/t2-1/t1
 if x==0:print("T1 and T2 differ")
 else:print("Hvap kJ/mol =",-R*log(p2/p1)/x/1000)

elif c==4:
 p2=pos("P2: ")
 t1=pos("T1 K: ")
 t2=pos("T2 K: ")
 h=pos("Hvap kJ/mol: ")*1000
 try:print("P1 =",p2/exp((-h/R)*(1/t2-1/t1)))
 except:print("P1 cannot display")

else:
 p1=pos("P1: ")
 p2=pos("P2: ")
 t2=pos("T2 K: ")
 h=pos("Hvap kJ/mol: ")*1000
 x=1/t2+(R/h)*log(p2/p1)
 if x<=0:print("T1 is not positive")
 else:
  t1=1/x
  print("T1 K =",t1)
  print("T1 C =",t1-273.15)

input("ENTER for guide: ")
print("QUICK GUIDE")
print("Use Kelvin only")
print("P units must match")
print("Higher T: higher P")
print("Hvap in kJ/mol")
print("DuckieDai: done!")
input("ENTER to exit: ")

from math import log,exp
R=8.314
def duckiedai_intro(program_name,wait=True):
title=program_name[:16]
empty=16-len(title)
left=empty//2
right=empty-left
print("+----------------------+")
print("| "+" "*left+title+" "*right+" |")
print("| |")
print("| __ |")
print("| ___(o )> quack! |")
print("| \\ <_. ) |")
print("| `---' |")
print("| by DuckieDai |")
print("+----------------------+")
print("Loading...")
if wait:
input("Press enter to start ")
print("\n"*6)
def pos(s):
while 1:
try:
x=float(input(s))
if x>0:return x
except:pass
print("Enter a positive #")
def menu():
while 1:
try:
x=int(input("Choose 1-5: "))
if x>0 and x<6:return x
except:pass
print("Choose 1 to 5")
duckiedai_intro("CLAUSIUS CLAP")
print("CLAUSIUS-CLAPEYRON")
print("1 Find P2")
print("2 Find T2")
print("3 Find Hvap")
print("4 Find P1")
print("5 Find T1")
c=menu()
if c==1:
p1=pos("P1: ")
t1=pos("T1 K: ")
t2=pos("T2 K: ")
h=pos("Hvap kJ/mol: ")*1000
try:
p2=p1*exp((-h/R)*(1/t2-1/t1))
print("P2 =",p2)
print("P2/P1 =",p2/p1)
except:print("P2 too large")
elif c==2:
p1=pos("P1: ")
p2=pos("P2: ")
t1=pos("T1 K: ")
h=pos("Hvap kJ/mol: ")*1000
x=1/t1-(R/h)*log(p2/p1)
if x<=0:print("T2 is not positive")
else:
t2=1/x
print("T2 K =",t2)
print("T2 C =",t2-273.15)
elif c==3:
p1=pos("P1: ")
p2=pos("P2: ")
t1=pos("T1 K: ")
t2=pos("T2 K: ")
x=1/t2-1/t1
if x==0:print("T1 and T2 differ")
else:print("Hvap kJ/mol =",-R*log(p2/p1)/x/1000)
elif c==4:
p2=pos("P2: ")
t1=pos("T1 K: ")
t2=pos("T2 K: ")
h=pos("Hvap kJ/mol: ")*1000
try:print("P1 =",p2/exp((-h/R)*(1/t2-1/t1)))
except:print("P1 cannot display")
else:
p1=pos("P1: ")
p2=pos("P2: ")
t2=pos("T2 K: ")
h=pos("Hvap kJ/mol: ")*1000
x=1/t2+(R/h)*log(p2/p1)
if x<=0:print("T1 is not positive")
else:
t1=1/x
print("T1 K =",t1)
print("T1 C =",t1-273.15)
input("ENTER for guide: ")
print("QUICK GUIDE")
print("Use Kelvin only")
print("P units must match")
print("Higher T: higher P")
print("Hvap in kJ/mol")
print("DuckieDai: done!")
input("ENTER to exit: ")

Filename: ccEQ.py
Size: 2167 bytes
SHA-256: 8cc5c176f0340878b9cc77812e2a04b48b406d072b87b10faab88b17d891f84f

Transfer and launch

Confirm that your calculator is the Python-capable model named above. This release has no verified calculator OS or Python App version; treat transfer and execution as unverified until independently tested.
Download the .py source above and verify its SHA-256 digest if your computer provides that option.
Use TI Connect CE to send the Python file to a compatible calculator. CCEQ is a suggested name; you may choose another valid, unique calculator name.
Open the Python App, select the program, and check sample inputs before relying on other results.

Read the complete installation, launch, troubleshooting, and removal guide.

Availability does not mean a teacher, school, or exam permits this program. Follow the applicable rules.

Version history

1.0.0 — current version, published 2026-09-23. Download this reviewed version

Any source change requires a new immutable version and digest. Historical downloads remain available only while their review records remain valid.