source: doc/theses/jiada_liang_MMath/test2.cfa@ 520fa9e

Last change on this file since 520fa9e was 7f8c232, checked in by Peter A. Buhr <pabuhr@…>, 15 months ago

formatting

  • Property mode set to 100644
File size: 1.6 KB
Line 
1#include <fstream.hfa>
2#include <stdlib.hfa>
3#include <enum.hfa>
4
5struct MR { double mass, radius; };
6
7enum( MR ) Planet { // typed enumeration
8 // mass (kg) radius (km)
9 MERCURY = { 0.330_E24, 2.4397_E6 },
10 VENUS = { 4.869_E24, 6.0518_E6 },
11 EARTH = { 5.976_E24, 6.3781_E6 },
12 MOON = { 7.346_E22, 1.7380_E6 }, // not a planet
13 MARS = { 0.642_E24, 3.3972_E6 },
14 JUPITER = { 1898._E24, 71.492_E6 },
15 SATURN = { 568.8_E24, 60.268_E6 },
16 URANUS = { 86.86_E24, 25.559_E6 },
17 NEPTUNE = { 102.4_E24, 24.746_E6 },
18 PLUTO = { 1.303_E22, 1.1880_E6 }, // not a planet
19};
20
21enum( double ) { G = 6.6743_E-11 }; // universal gravitational constant (m3 kg-1 s-2)
22
23static double surfaceGravity( Planet p ) with( p ) {
24 return G * mass / ( radius \ 2 ); // no qualification, exponentiation
25}
26static double surfaceWeight( Planet p, double otherMass ) {
27 return otherMass * surfaceGravity( p );
28}
29
30int main( int argc, char * argv[] ) {
31 if ( argc != 2 ) exit | "Usage: " | argv[0] | "earth-weight"; // terminate program
32
33 double earthWeight = convert( argv[1] );
34 double earthMass = earthWeight / surfaceGravity( EARTH );
35
36 Planet rp = fromInt( prng( countof( Planet ) ) ); // select random orbiting body
37 choose( rp ) { // implicit breaks
38 case MERCURY, VENUS, EARTH, MARS:
39 sout | rp | "is a rocky planet";
40 case JUPITER, SATURN, URANUS, NEPTUNE:
41 sout | rp | "is a gas-giant planet";
42 default:
43 sout | rp | "is not a planet";
44 }
45
46 for ( p; Planet ) { // enumerate
47 sout | "Your weight on" | ( p == MOON ? "the" : " " ) | p
48 | "is" | wd( 1,1, surfaceWeight( p, earthMass ) ) | "kg";
49 }
50}
Note: See TracBrowser for help on using the repository browser.