/* kp.c  copyright (c) by James N. Hall 08/02/2000 */

#include <stdlib.h>
#include <math.h>


int main () {
       float vel = 500; /* velocity */
       float mass = 20; /* mass */
       float z = 225218; /* slug to grains constant may be accurate enough*/
       float en   ; /* energy */
       float mv   ; /* momentum */
       float kp   ; /* kinetic pulse */
       float sem   ; /* sem standard EM */
       int x,y; /* utility variable to increment*/
     sem = 0; 
    for ( y = 0; y <= 4000 ; y = y + 50 ) {
       printf("KINETIC PULSE   WEIGHT   VELOCITY       ENERGY      MOMENTUM  Standard E/M \n");
       for ( x = 0; x <= 600 ; x = x + 5 ) {
        en = ( (mass + x ) * ( vel + y ) * ( vel + y ) * .5)  / ( z ) ;
        mv = ( (mass + x )  * ( vel + y ))  / (  z ) ;
        kp = ( (mass + x ) * ( mass + x ) * 
             ( vel + y ) * ( vel + y ) * (vel + y) * .5  )  / ( z * z ) ;
	if ( mv == 0 ) sem  = 0 ; 
	else sem = ( (vel + y) / 2 ) * ( 1 / mv );
       printf("   %6.2f kp %6d gr %6d f/s %6.0f ft-lb %3.2f slug f/s %6.2f SE/M \n",
           kp , (int) mass + x , (int) vel + y , en, mv , sem  );
       }
    }
}