home *** CD-ROM | disk | FTP | other *** search
/ POINT Software Programming / PPROG1.ISO / pascal / swag / math.swg / 0058_Gravity.pas < prev    next >
Encoding:
Pascal/Delphi Source File  |  1994-01-27  |  2.5 KB  |  66 lines
  1. (*
  2. >Does anyone have any equations for gravity??
  3.  
  4. It's not as tough as you probably think it is.  The way I work motion in
  5. my programs is, I keep track of the acceleration, velocity, and position
  6. of an object in both the x and y directions.  In other words, I have these
  7. variables:
  8.  
  9.   var ax, ay, vx, vy, px, py: integer;
  10.  
  11. When you have a force -- like gravity, or wind resistance, or whatever --
  12. you need to recalculate the accelerations every game round.  Then you
  13. alter the velocities accordingly, and after that you change the positions.
  14. For example, each round you execute code like this:
  15.  
  16.   ax := {formula for force in the "x" direction};
  17.   ay := {formula for force in the "y" direction};
  18.  
  19.   vx := vx + ax;
  20.   vy := vy + ay;
  21.  
  22.   px := px + vx;
  23.   py := py + vy;
  24.  
  25. Notice how simple it is to keep track of motion: all you need to do is
  26. supply a formula for acceleration, and the program runs "blind" after
  27. that point.
  28.  
  29. So gravity is just a matter of supplying the right "acceleration" formulas.
  30. If you are talking gravity near the surface of the earth, gravity provides
  31. very nearly a constant acceleration.  In which case:
  32.  
  33.   ax := 0;  {no "horizontal" gravity}
  34.   ay := g;  {a constant -- assign whatever value you like}
  35.  
  36. For objects to fall "down" the screen, "g" should be positive.  Motion
  37. towards the top of the screen would mean a negative velocity.  That's
  38. because "y" coordinates increase from top to bottom, and frankly that
  39. confuses me and it confuses the numbers.  You might do well to do this:
  40. have your calculations assume that "y" coordinates increase from bottom
  41. to top, and then draw at position (px, GetMaxY + 1 - py).  With coordinates
  42. increasing from bottom to top, "g" should be negative and upward motion
  43. means positive "vy".
  44.  
  45. If you want gravity as applies to celestial objects in orbit, the formulas
  46. for acceleration would be:
  47.  
  48.   x := px - sx;  { new variables: sx and sy are the locations of the sun or }
  49.   y := py - sy;  { whatever, and x and y are thus the distances from it }
  50.  
  51.   ax := g*x / exp(3*ln(x*x + y*y)/2);
  52.   ay := g*y / exp(3*ln(x*x + y*y)/2);
  53.  
  54. Again, I recommend plotting at (px, GetMaxY + 1 - py); and again, "g"
  55. should be negative.
  56.  
  57. Be advised that there is a singularity at the location of the sun or
  58. whatever: the "ln" calculations will fail.  Another gravity formula I've
  59. seen used is "bowl" gravity, like a marble rolling around in a bowl.  It's
  60. unrealistic, but it "feels" good and doesn't have a singularity.  In which
  61. case:
  62.  
  63.   ax := g*x;  { negative "g" again }
  64.   ay := g*y;
  65.  
  66. *)