share/ezunits/ezunits.dem

   1 if get ('ezunits, 'version) = false then load (ezunits) else 'already_loaded;
   2
   3 "... SI units ..."$
   4
   5 foo : 10 ` m;
   6
   7 qty (foo);
   8
   9 units (foo);
  10
  11 dimensions (foo);
  12
  13 "... \"Customary\" units ..."$
  14
  15 bar : x ` acre;
  16
  17 dimensions (bar);
  18
  19 fundamental_units (bar);
  20
  21 "... Unit conversions ..."$
  22
  23 1 ` ft `` m;
  24
  25 %, numer;
  26
  27 1 ` kg `` lbm;
  28
  29 %, numer;
  30
  31 1 ` W `` Btu/hour;
  32
  33 %, numer;
  34
  35 100 ` degC `` degF;
  36
  37 -40 ` degF `` degC;
  38
  39 1 ` acre*ft `` m^3;
  40
  41 %, numer;
  42
  43 "... Arithmetic operations on dimensional quantities ..."$
  44
  45 100 ` m + 100 ` ft;
  46
  47 (100 ` m + 100 ` ft) `` ft;
  48
  49 %, numer;
  50
  51 (100 ` m + 100 ` ft) `` m;
  52
  53 %, numer;
  54
  55 ((10 ` acre) * (10 ` ft)) `` m^3;
  56
  57 %, numer;
  58
  59 quux : (144 ` Btu/hour) / (48 ` ft^2);
  60
  61 dimensions (quux);
  62
  63 fundamental_units (quux);
  64
  65 quux `` fundamental_units (quux);
  66
  67 foo^2;
  68
  69 dimensions (foo^2);
  70
  71 foo^2 `` acre;
  72
  73 %, numer;
  74
  75 bar `` m^2;
  76
  77 %, numer;
  78
  79 sqrt (bar);
  80
  81 dimensions (sqrt (bar));
  82
  83 "... Units ad hoc ..."$
  84
  85 baz : 3 ` sheep + 8 ` goat + 1 ` horse;
  86
  87 subst ([sheep = 3*goat, horse = 10*goat], baz);
  88
  89 baz2 : 1000`gallon/fortnight;
  90
  91 subst (fortnight = 14*day, baz2);
  92
  93 % `` liter/hour;
  94
  95 "... Declared unit conversions ..."$
  96
  97 declare_unit_conversion (MMBtu = 10^6*Btu, kW = 1000*W, kWh = kW*hour, MWh = 1000*kWh);
  98
  99 1 ` MMBtu `` MWh;
 100
 101 1 ` MMBtu/day `` kW;
 102
 103 "... Stefan-Boltzmann constant in customary units ..."$
 104
 105 56697/10000/10^8 ` W/m^2/K^4;
 106
 107 % `` Btu/hour/ft^2/R^4;
 108
 109 %, numer;
 110
 111 "... Velocity of sound in helium at 27 C, He(Cp/Cv) = 5/3, 4 g/mol ..."$
 112
 113 if get ('physical_constants, 'version) = false then load (physical_constants) else 'already_loaded;
 114
 115 v : sqrt (5/3 * %R/(4 ` g/mol) * (27 ` degC `` K));
 116
 117 get (%R, 'description);
 118
 119 constvalue (%R);
 120
 121 constvalue (v);
 122
 123 dimensions (v);
 124
 125 fundamental_units (v);
 126
 127 constvalue (v) `` %;
 128
 129 %, numer;
 130