#include <stdio.h>
#include "glob.h"

// function prototype
void prob_01();

int main()
{
  printf("c outout My problem 1\n"); 
  prob_01();
}

void prob_01()
{
  double M; // Mass in which h = c = me = 1.0
  double L; // Length in which h = c = me = 1.0
  double T; // Time unit in which h = c - me = 1.0

  // Derivation (done in algebra by hand.)

  // the rules here me, h, and c are just numbers

  // M / kg = me
  // M L^2 / T * sec / (kg meters^2) = h
  // L/T sec/meters = c

  // So
  // M = me kg
  // M L^2 / T = h kg meters^2/sec
  // (M L ^2) / T * (T / L) meters/sec = h / c kg meters^2 sec <- this was
  // cancelling                                                !! error 
  // M L = h c kg meters                                       ! c was in 
  // and                                                       ! numerator
  // M L / M = h / c kg meters / (me kg)                         
  // cancelling 
  // L = h / c / me meters
  // (L / T) (sec / meters) = c
  // T c = L sec / meters
  // T = (L / c) (sec / meters)
  // T = (h / c / me meters) / c (sec / meters)
  // cancelling
  // T = h / (c ^ 2 / me) seconds
  //

  M = me; 
  L = h / (c * me);
  T = h / (c * c * me);  // multiplication is much more efficient
  
  printf("Units so c = 1.0, h = 1.0 and me = 1.0\n"); 
  printf("Unit of mass in kilograms = %16.10e\n",M);
  printf("Unit of length in meters  = %16.10e\n",L);
  printf("Unit of time in seconds   = %16.10e\n",T);
}