Natural Science Forum / Physics / Relativity / September 2005

It's not irrelevant. It is most pertinent to what you are doing.
If you were to calculate the areas of the triangles between your point
pairs and the focus you'd see a trend.
I know, because I went that way myself years ago.

Don't be silly, I developed it years ago. I know exactly what you are
doing.
When you first began you were programming circles.

They'll do that with a circle, but you won't get a sensible curve.

That's a silly statement to make. If you don't have equal areas in equal
times,
AND YOU DONT, you'll change the shape of the light curve.

Yours will precisely precess.

Bullshit.

Divisions are:
v=delta_x/interval, that can't be zero
Phase =  (EinsteinTime-GalileoTime + time)/Period; That can't be zero.
Doppler =(c_dash-v)/c_dash, that can't be zero;
EinsteinTime = D/c;  That can't be zero.
GalileoTime = D/(c+v); That can't be zero.

You must think I was born yesterday. I've has a career as a software
engineer,
Show my divide by zero.
You can't because you don't know what you are babbling about.

void Orbit()
{
double time;
double Angle,A = 0;
double x,y,z,old_x,delta_x,v,D,EinsteinTime, GalileoTime,Doppler,line,
c_dash, count;
long i,j,k, step, arrive;
char s[80];
double Frequency;

RECT Rect;
HDC hDC;
BOOL flag;

ClearArrays(); // empty out previous data.
CosPitch = cos(pitch);  //precompute to save time
SinPitch = sin(pitch);
CosYaw = cos(yaw);
SinYaw = sin(yaw);
CosRoll = cos(roll);
SinRoll = sin(roll);
SemiMinorAxis = SemiMajorAxis * sqrt(1.0 - eccentricity*eccentricity);
MajorAxis = 2.0*SemiMajorAxis;
interval = Period/iterations;
GetClientRect(hWnd,&Rect);
hDC= GetDC(hWnd);
FillRect(hDC, &Rect, GetStockObject(BLACK_BRUSH));
ReleaseDC(hWnd, hDC);
out_of_range = FALSE;

Angle = 0;
x = SemiMajorAxis;
y = 0;
z = 0;
Rotate(&x,&y,&z);
old_x = x;  // initial value
for(k=0;k<2;k++)
{
step = iterations/ARRAYSIZE;                  // the "distance" between
points saved
for(i = 0,count=0.0;i<iterations; i++, count=i*ARRAYSIZE/iterations) //
set up one period in increments
  {
 time = Period *(double)i/(double)iterations;
 Angle = Kepler((double)i/(double)iterations*2.0*PI); // get one
position of the ellipse
 x = SemiMajorAxis*cos(Angle);       // compute as x and y.
 y = SemiMinorAxis*sin(Angle);
 z = 0;
 Rotate(&x,&y,&z);                // pitch and yaw the ellipse
    delta_x = old_x-x;         // distance source move toward observer
 old_x = x;               // save for next computation
 v=delta_x/interval;         // speed source moves toward observer
       if(!out_of_range)
       if ((v > 3000) || (v< -3000))
        {
        sprintf(s,"Parameters give an improbable velocity of%6.0lf
Km/sec. This may result in a long computation time",v);
        flag = MessageBox(hWnd,s,"   Velocity out of range",
MB_OKCANCEL);
        if(flag == 1)  out_of_range = TRUE;
        else return;
        }
 D = Distance+x;               // Distance to observer
 EinsteinTime = D/c;
 GalileoTime = D/(c+v);               // Time it takes light to reach
observer
    if(Androcles)
     {
     Phase =  (EinsteinTime-GalileoTime + time)/Period;
     }
    else Phase = time/Period;
    arrive = (int)((Phase)*(double)ARRAYSIZE);
 ArrivalTime[(long)count] = arrive+ARRAYSIZE/3;
    while (Phase<0.0) Phase+=1.0;
    while (Phase>=1.0) Phase-=1.0;
    //don't arrive before we leave
    arrive = (int)((Phase)*(double)ARRAYSIZE);
    LightCurveArray[arrive]++;// Drop a "photon" into the indexed
array, not in any order
 if(i%step == 0) // index1 goes from 0 to iterations, count from 0 to
ARRAYSIZE
  {
  EllipseX[(long)count] = x;
  EllipseY[(long)count] = y;
  EllipseZ[(long)count] = z;
  Velocity[(long)count] = v;
  if(MoveSystem(x,y,z)) return; // request to halt processing recieved
  }
       if(k==0)      // only on the first pass
        {
        if(magnify > 1.0)      // artificially make c smaller to
exaggerate Doppler
   {
   c_dash = c/(magnify);   // arbitrary
   if(Androcles) Doppler =(c_dash-v)/c_dash;
         else if(Einstein) Doppler =
sqrt((1.0+v/c_dash)/(1.0-v/c_dash));
         else Doppler = c_dash/(c_dash-v);  //Lorentz aether;
         }
        else
         {
   if(Androcles) Doppler =(c-v)/c;
         else if(Einstein) Doppler = sqrt((1.0+v/c)/(1.0-v/c));
         else Doppler = (c)/(c-v);  //Lorentz aether;
         }

        //now populate the visible Spectrum array from 0-255 with bytes
proportional to Doppler shift
    for(j=0;j<10;j++)
      {
      if(arrive/4>ARRAYSIZE/4-10)
       {

       //line = (246.0 - 10.0)*line/(750.0-400.0)-246.0; //scaled to
pixels,256 to 0
       line = (246.0 - 10.0)*line/(750.0-400.0)-246.0; //scaled to
pixels,256 to 0
       line = 255.0-line;
      }
      else
       {
      Frequency = Doppler * c/series[j]; // series is wavelength
                   line = c/Frequency; // back to wavelength
       line = (246.0 - 10.0)*line/(750.0-400.0)-246.0; //scaled to
pixels,256 to 0
       line = 255.0-line;
       }

         if(line > 0.0 && line <256.0)
          Spectrum[arrive/4][(int)line]= (BYTE)TRUE;
         }
  }// end first pass, k = 0
 } //end for(i... loop
} // end for(k.... loop
Magnitude();
FilterCurve();
}

That will get you a padded cell, but go ahead, set yourself up.

I've always said it was a ternary system, H, but you've got 48 cycles in
70 years
or 0.686 cycles a year, 532 days. Plenty of time and ROOM for a
Sun-Earth-Moon
triple, which you have not modelled. Doesn't the rest of the curves past
the three you modeled mean something without you jumping to conclusions?

Using Kepler's equation, move the barycentre of the Earth-Moon system
around the Sun, then move the Moon around its own local ellipse. Then
see what kind of curve you'll get from the Sun.

They can only come from the spectra.
I'd like to help you, but my ankle won't let me get to a decent library.
Doc is going to operate in 5 weeks, perhaps I'll be less of a grouch
by Christmas. Don't you have anything in Sydney? There must be a
major university there.

Hmmm.... the core of a star is cooler... I think I'll reserve judgment
on that,
mainly because I've never modelled IR anyway. Where did you get IR
spectral velocities?
Androcles