Ballistic Theory, Progress report...Suitable for 5yo Kids
in [Relativity]
|
Prev: OWLS is not equal to c
Next: Mathematical Inconsistencies in Einstein's Derivation of the Lorentz Transformation
From: George Dishman on 4 Oct 2005 18:06 "Henri Wilson" <H@..> wrote in message news:0tt5k1tb7afbn1f7d7avltq14eof0kujqe(a)4ax.com... > On Tue, 4 Oct 2005 11:27:40 +0100, "George Dishman" > <george(a)briar.demon.co.uk> > wrote: > >> >><jgreen(a)seol.net.au> wrote in message >>news:1128419173.735037.311150(a)g49g2000cwa.googlegroups.com... >>> > >>> >>> George, my old computer died on me, losing the email I remember quite >>> vividly (from you) saying how sagnac machine works because the TIME OF >>> TRAVEL of the signal alters. If you cannot see/understand that this >>> refers to VELOCITY, as the DISTANCE does NOT ALTER, that is too bad! >> >>Sadly it seems you also lost my response to your comment, >>the distance does alter Jim, the detector MOVES while the >>light is in transit. In the experiment, the length is >>known to change because we know the speed of rotation of >>the table, the time is measured to change and when you >>calculate the speed as distance/time you always get c. > > George, that's the aether explanation. No, it isn't an explanation at all Henri, it is a statement of fact. The detector moves in the lab frame, period. >>In the commercial products, knowing the speed is c, the >>device measures the time difference and calculates the >>rotation and they work perfectly. >> >>As for your stuff on the galaxy and a merry-go-round, the >>correct analogy is that you make the measurement not against >>the horses head but using a gyroscope (or you could say a >>distant mountain if you were on a non-rotating planet). >>Astronomers are well aware of the local proper motion of >>stars. That's why you need to learn how astronomy is done >>before criticising. > > They eventually learnt to allow for a CONSTANT light speed. No, they just got better instruments that were able to measure the very low rotation rates involved. Look up the ICRF. > The still haven't woken up to the fact that all starlight isn't > miraculously > adjusted (by the fairies) to travel to little planet Earth at exactly 'c'. Flaunting your inability to comprehend something as simple as SR isn't helping your case. George
From: Henri Wilson on 4 Oct 2005 18:02 On Sun, 02 Oct 2005 21:53:08 +0200, "Paul B. Andersen" <paul.b.andersen(a)deletethishia.no> wrote: >Henri Wilson wrote: >> >> A, I have supported you on this. In Algol's case, the WCH happens to be the >> large planet 'Androcles'. Do you have any objections to that? >> It is likely that flares occuring on the main star are reflected from >> 'Androcles' and mistaken for flares on IT. > >Very likely, indeed. :-) > >There IS no limit to the stupidities you can utter, is there? :-) >You have no idea of what a flare is, do you? :-) Of course I do. You obviously don't. > >> However I might be inclined to agree with the tusselad that a small third >> object does orbit the main Algol star. > >The "third object" is not a particularly small star. >It is an A5 V. That means it is bigger and brighter than the Sun. >But it is smaller than Algol A and B. > >> On the published curve, there is a small departure from our predicted curves >> that might be associated with an object orbiting with the same period but >> lagging in phase behind the main star. > >The orbital period of Algol C is 681 days, that is 235 times >the period of the inner binary. Well it would not affect the brightness curve significantly. > >The "small departure" from your predicted curves, >is the second minimum when Algol A eclipses Algol B. there is probably NO Algol B. > >Your "drawing program" is unable to mimic that because >the degrees of freedom is too small. So you are unable >to invent any fantasy data that will produce the correct >light curve. Here is some of the code. Please tell the world which part constitutes the 'drawing program'. Dim stoppoint, temppos, temppos2, startpoint, lineYs, eccindex, pointindex, molecV1, sinechange Dim maxspeed, minspeed, p, j, K, L, s, t, falsep, lag, Ylag, Xlag, cosvan As Single, ecc As Variant Dim lightspeed(70000), position(70000), circlerun, Pref, prefmod As Integer, frac, velroot As Single Dim velocity(70000) As Double, Vangle(70000) As Double, LS(70000), pos(70000), lastj, lastk Dim imaghite, outerperiod, outerp, velratio, twovelrat, scaling, points, yaw, piyaw, piyawb, dlim, stepint, lineY Dim W, vone, voneA, rone, rtwo, Vonec, Vtwo, Vtwoc, b, c, d, e, f, u, m, n, z, v, v2, vb, q, r, velint As Single Dim timefraction, dfraction, jpi, vell0, vell90, vell180, vell270, velfactor, pointmin As Single Dim orbs As Integer, A As Integer, x, vellxovery, ninet, twenine, lumrat, pimod, Gmod Dim pointstwo As Long, roll As Double, yawadd As Integer, sec, stepsize As Integer Dim Y As Double, YB As Double, cosy As Double, cosroll As Double, cosrollV As Double, cosrollVA As Double Dim cut As Boolean, incmolvec As Boolean, tempfac Dim pi As Double, twopi, elipsize, pointnumber As Integer, expercent, xfactor Dim Tmax As Double, linearcurve As Boolean, logcurve As Boolean Dim Tmin As Double, Xshort, Yshort, ForRev As Boolean, Maginner, magouter, magboth Dim BmaxInner, BmaxOuter, BmaxBoth, BminInner, BminOuter, BminBoth, PP As Double Dim Tat(70000) As Double, Extspeed(70000) As Double, extS(70000) As Double Dim Tat2(70000) As Double, extcircle(70000) As Double, circleS(70000) As Double Dim Tbt(70000) As Double, hite(400) Dim Ta As Double, xco(150), yco(150) Dim Tb As Double, LR, UD, SIZE, WIDE, TAmin As Double, TAmax As Double Dim Stardist As Double Dim dtemp As Single Dim intensity(600) Dim INTB(600) Dim scanintA(4, 12, 800), scanintB(4, 12, 800), added(4, 12, 800) As Integer Dim orbtime(100) Dim vell(4000) Dim vellangle(4000) Dim VelB(4000) Dim velangleB(4000) Dim vellx As Double Dim velly As Double Dim Xtemp As Double Dim Ytemp As Double Dim xstart, Ystart As Integer, FocusX As Integer, FocusY As Integer Dim G As Double, pmax Dim Xacc As Double Dim Yacc As Double Dim Force As Double, Radvector As Double Private Sub Command10_Click() End End Sub Private Sub Command9_Click() If Label31.Visible = False Then Label31.Visible = True Label17.Visible = False Command9.Caption = "Close Instructions" Else Label31.Visible = False Command9.Caption = "Instructions" End If End Sub Private Sub Cpulses_Click() If Pulsemove.Visible = False Then Pulsemove.Visible = True Timer1.Enabled = False Cpulses.Caption = "Normal" Call Pulsepos Else Pulsemove.Visible = False Cpulses.Caption = "Pulses" Timer2.Enabled = False End If End Sub Public Sub Form_Load() Randomize Combo1.AddItem 0.2 'this magnifies the time scale. Combo1.AddItem 0.5 Combo1.AddItem 1 Combo1.AddItem 2 Combo1.AddItem 5 Combo1.AddItem 30 CMaxvel.AddItem 0.000001 'peripheral velocity, c=1 CMaxvel.AddItem 0.00001 CMaxvel.AddItem 0.00003 CMaxvel.AddItem 0.0001 CMaxvel.AddItem 0.0003 CMaxvel.AddItem 0.001 CMaxvel.AddItem 0.01 CMaxvel.AddItem 0.1 Combo3.AddItem 0.003 'period. years Combo3.AddItem 0.01 Combo3.AddItem 0.03 Combo3.AddItem 0.1 Combo3.AddItem 0.3 Combo3.AddItem 0.7 Combo3.AddItem 1 Combo3.AddItem 3 Combo3.AddItem 10 CTempfac.AddItem 100 CTempfac.AddItem 200 CTempfac.AddItem 500 CTempfac.AddItem 1000 CTempfac.AddItem 1500 CTempfac.AddItem 3000 Combo4.AddItem 0 ' Startpoint, x1000 LYs Combo4.AddItem 0.03 Combo4.AddItem 0.06 Combo4.AddItem 0.1 Combo4.AddItem 0.3 Combo4.AddItem 0.6 Combo4.AddItem 1 Combo4.AddItem 2 Combo4.AddItem 5 Combo4.AddItem 20 Combo4.AddItem 100 Combo4.AddItem 1000 Combo5.AddItem 0 'orbit roll, rotation around major axis Combo5.AddItem 15 Combo5.AddItem 30 Combo5.AddItem 60 Combo6.AddItem 20000 Combo6.AddItem 33000 Combo6.AddItem 60000 Combo7.AddItem 1# 'Orbs, number of Orbits for calc. Combo7.AddItem 2# Combo7.AddItem 3# Combo7.AddItem 8# Combo7.AddItem 15# Combo7.AddItem 25# Combo8.AddItem 0.25 Combo8.AddItem 0.5 Combo8.AddItem 1 Combo8.AddItem 2 Combo8.AddItem 5 'a scaling factor to bring output Combo8.AddItem 20# 'into screen Combo9.AddItem 180 'Yaw angle Combo9.AddItem 150 Combo9.AddItem 120 Combo9.AddItem 90 Combo9.AddItem 60 Combo9.AddItem 30 Combo9.AddItem 0 Combo9.AddItem -30 Combo9.AddItem -60 Combo9.AddItem -90 Combo9.AddItem -120 Combo9.AddItem -150 Combo10.AddItem 4 'relative luminosity inner/outer star Combo10.AddItem 3 Combo10.AddItem 2 Combo10.AddItem 1.7 Combo10.AddItem 1.5 Combo10.AddItem 1.3 Combo10.AddItem 1# Combo10.AddItem 0.8 Combo10.AddItem 0.6 Combo10.AddItem 0.4 Combo10.AddItem 0.2 eccentricity.AddItem 0 'eccentricity eccentricity.AddItem 0.06 eccentricity.AddItem 0.1 eccentricity.AddItem 0.15 eccentricity.AddItem 0.2 eccentricity.AddItem 0.25 eccentricity.AddItem 0.3 eccentricity.AddItem 0.4 eccentricity.AddItem 0.5 eccentricity.AddItem 0.6 eccentricity.AddItem 0.7 eccentricity.AddItem 0.75 eccentricity.AddItem 0.8 Combo11(0).AddItem 0 ''velratio' R2/R1 Combo11(0).AddItem 0.1 Combo11(0).AddItem 0.2 Combo11(0).AddItem 0.4 Combo11(0).AddItem 0.6 Combo11(0).AddItem 0.8 Combo11(0).AddItem 1# Clag.AddItem 0 'phase lag Clag.AddItem 30 Clag.AddItem 60 Clag.AddItem 75 Clag.AddItem 90 Clag.AddItem 120 Clag.AddItem 150 Clag.AddItem 180 Clag.AddItem 210 Clag.AddItem 240 Clag.AddItem 270 Clag.AddItem 300 Clag.AddItem 330 Cextinction.AddItem 0.001 'unification factor, fractional unification per LY Cextinction.AddItem 0.002 Cextinction.AddItem 0.005 Cextinction.AddItem 0.01 Cextinction.AddItem 0.02 Cextinction.AddItem 0.05 pi = 3.141592653575 twopi = 2 * pi pimod = pi / 180 cut = False vellx = 0 velly = 0.001 FocusX = 0: FocusY = 0 Ystart = 0 Xtemp = xstart: Ytemp = 0 Radvector = (((FocusX - Xtemp) ^ 2) + ((FocusY - Ytemp) ^ 2)) ^ 0.5 'length of radius vector Label17.Visible = True p = 0 sec = 0 pointindex = 2.6 Erase velocity Erase Vangle star1 = True star2 = True End Sub Public Function combovalues() ecc = eccentricity n = Combo1 velratio = Combo11(0) falsep = Combo3 tempfac = CTempfac / 90 / 10000000 roll = Combo5 * pimod orbs = Combo7 imaghite = Combo8 voneA = CMaxvel vone = voneA * 1000 yaw = Combo9 lag = Clag velroot = Combo10 pointindex = Combo6 expercent = (Cextinction) xfactor = 1 - expercent startpoint = 1000 * Combo4 piyaw = pimod * yaw rootc = 2.512 logrootc = 0.921 c = 2.997 * (10 ^ 5) 'kilometres per second End Function Private Sub Command5_Click() 'read text If Label17.Visible = False Then Timer1.Enabled = False Label17.Visible = True Command5.Caption = "Close Text" Else: Label17.Visible = False Command5.Caption = "Read Text" End If End Sub Private Sub Command7_Click() 'back to startpage Timer1.Enabled = False Unload truevel Startpage.Show End Sub Public Sub Cstart_Click() 'start program. red button Randomize Timer1.Enabled = False Call combovalues Erase lightspeed Erase position Erase orbtime Erase LS Erase pos Erase vell Erase VelB Erase xco Erase yco If Option5.Value = True Then Combo4.Enabled = False Else Combo4.Enabled = True Label17.Visible = False Label44.Visible = False Label31.Visible = False Label19.Left = -1200 Label29.Left = -1200 Command9.Caption = "Instructions" Command2.Caption = "PAUSE" cut = True Stardist = 0 If eccentricity.Text = "" Then eccindex = 20 'this is to generate an error message if eccentricity is changed Label35.Visible = True Else Label35.Visible = False p = 1 eccindex = ecc truevel.Cls details.Cls DrawWidth = 2 Erase velocity Erase Vangle pmax = 0 vellx = 0 velly = 0.001 If ecc = 0 Then details.Show details.Circle (150, 120), 75 details.Line (50, 100)-(150, 100), RGB(0, 0, 255) details.Line (50, 50)-(150, 150), RGB(255, 0, 0) details.PSet (150, 120) points = 20000 GoTo 10: End If '**************************** new 'create and draw ellipse 'Timer1.Enabled = False p = 1 details.Cls DrawWidth = 2 If Combo6.Text = 20000 Then pointnumber = 2.6 stepsize = 20 End If If Combo6.Text = 300 Then pointnumber = 4.54 stepsize = 20 End If If Combo6.Text = 60000 Then pointnumber = 8.7 stepsize = 100 End If Erase velocity Erase Vangle pmax = 0 ecc = eccentricity xstart = pointnumber * (1 - ecc) ^ (5 / 3) 'empirically determined to give the about the right number of points G = xstart / (1 + ecc) / (10 ^ 6) 'ditto elipsize = 750 / xstart * ((1 - ecc) ^ (2 / 3)) 'empirically adjusts print size Ystart = 0 Xtemp = -xstart: Ytemp = 0 Radvector = ((Xtemp ^ 2) + (Ytemp ^ 2)) ^ 0.5 'length of radius vector vellx = 0 velly = 0.001 Force = G / ((Xtemp ^ 2) + (Ytemp ^ 2)) '(Radvector ^ 2) DrawWidth = 1 While Ytemp > -0.00000001 'And Fix(1000 * Xtemp) < 0 If p Mod 4 = 0 Then PSet (2000 + (elipsize * Xtemp), 1500 + (Ytemp * elipsize)), RGB(255, 255, 255) 'draw ellipse Xacc = Force * (-Xtemp) / Radvector Yacc = Force * (Ytemp) / Radvector Xtemp = Xtemp + vellx '+ (0.5 * Xacc) Ytemp = Ytemp + velly '- (0.5 * Yacc) Radvector = ((Xtemp ^ 2) + (Ytemp ^ 2)) ^ 0.5 velocity(p) = ((vellx ^ 2) + (velly ^ 2)) ^ 0.5 If 100000000 * vellx <> 0 Then x = (velly / vellx) Vangle(p) = (pi / 2) - Atn(x) 'Call XYCoords K = Vangle(p - 1) End If vellx = vellx + Xacc velly = velly - Yacc If Abs(velly) < 0.000001 Then 'stop point half way round lineytemp = p 'this is for setting the 90 and 270 degree lines at the right place End If Force = G / (Radvector ^ 2) p = p + 1 Wend pmax = p - 1 points = (2 * p) '+ 1 lineY = lineytemp * 4000 / points Line3.Y1 = 315 + lineY: Line3.Y2 = 315 + lineY Line5.Y1 = 4315 - lineY: Line5.Y2 = 4315 - lineY Label10.Top = Line3.Y1 + 10: Label12.Top = Line5.Y1 + 10 LA.Left = 1980 - (elipsize * xstart) LC.Left = 1980 + ((Xtemp - xstart) * elipsize / 2) LB.Left = 1980 + (Xtemp * elipsize) LF.Visible = True LA.Visible = True LC.Visible = True LB.Visible = True details.Show For m = 1 To 21 details.Print Next details.Print "AC = "; Xtemp + xstart; " FC ="; Xtemp - xstart: details.Print "Minimum velocity= "; CMaxvel * 1000 * velocity(p - 1): details.Print "points="; points: 'pmax; details.Print "eccentricity="; (Xtemp - xstart) / (Xtemp + xstart): If ecc = 0 Then r = CMaxvel * falsep / 2 / pi details.Print "Radius ="; r; "LYs": End If For j = 1 To pmax - 1 'other half of ellipse velocity(p + j - 1) = velocity(p - j) Vangle(p + j - 1) = (2 * pi) - Vangle(p - j) Next Call diagram velocity(0) = 0.001 velocity(p) = 0.001 * (1 - ecc) / (1 + ecc) Vangle(0) = 0: 'XVel(0) = 0 Vangle(p) = pi details.Show piyaw = pimod * yaw Gmod = G * ecc / 0.25 For j = 0 To points Step 20 ''velocity v/s time graph details.PSet (50 + (j * 100 / p), 300 - (0.1 / Gmod * velocity(j))), RGB(0, 0, 255) details.PSet (50 + (j * 100 / p), 200 - (0.05 / Gmod * velocity(j) * Cos(Vangle(j) - piyaw))), RGB(0, 0, 0) details.PSet (50 + (j * 100 / p), 150 - (20 * Vangle(j))), RGB(255, 0, 0) Next 10: End If 'call coords 'prints out x and y coordinates of ellipse End Sub Public Function diagram() If ecc <> 0 Then 'Draw ellipse bottom corner. ecc = eccentricity Shape2.Visible = True: Shape1.Visible = True Shape1.Shape = 2 Shape1.Width = 600 * ((1 - (ecc ^ 2)) ^ 0.5) Shape1.Left = 1260 - (Shape1.Width / 2) Shape2.Top = 7220 - (300 * ecc) Shape2.Left = Shape1.Left - 15 + (Shape1.Width / 2) yaw = Combo9 Line7.Visible = True Line7.X1 = 15 + Shape2.Left: Line7.Y1 = 15 + Shape2.Top Line7.X2 = Line7.X1 + (1000 * Cos(yaw * pimod)) Line7.Y2 = Line7.Y1 + (1000 * Sin(yaw * pimod)) Label34.Top = Line7.Y2: Label34.Left = Line7.X2 End If End Function '*************************************************************************** '*************************************************************************** Public Sub lightfront_click() 'start light front demo. Green button Label17.Visible = False Call lightfronts 'If Option8.Value = True Then Call lightfronts Else Call extinction End Sub Public Function lightfronts() If eccentricity.Text = "" Then 'Stops a crash Label35.Visible = True Else If cut = False Then Label44.Visible = True GoTo out End If ecc = eccentricity.Text If ecc <> eccindex Then 'Stops a crash Label44.Visible = True cut = False GoTo out End If If Option5.Value = True Then Combo4.Enabled = False Else Combo4.Enabled = True 'DISABLE STRTING DISTANCE Call allvalues If ecc = 0 Then GoTo circlerun 'for circle Line1.Visible = False Line9.Visible = False Label43.Visible = True: Label50.Visible = True j = 0 cosyaw = Cos(yaw * pimod) minspeed = 1 maxspeed = 1 For K = 0 To points Step stepsize 'time unit=1 orbit If Option1.Value = True Then molecV1 = tempfac * (Int((830 / ((Rnd * 100) + 1))) - 8) 'introducing molecular velocities If Rnd < 0.5 Then sinechange = -molecV1 Else sinechange = molecV1 End If Cospiyaw = Cos(Vangle(K) - piyaw) * cosrollV lightspeed(K) = 1 + sinechange + (velocity(K) * Cospiyaw) 'emitted light speed. x c LS(K) = 1 + sinechange - (velocity(K) * velratio * Cospiyaw) 'second star If lightspeed(K) > maxspeed Then maxspeed = lightspeed(K) lastk = K End If If lightspeed(K) < minspeed Then minspeed = lightspeed(K) lastj = K End If position(K) = (lightspeed(K) * (1 + (startpoint / falsep))) - (K / points) ' pos(K) = (LS(K) * (1 + (startpoint / falsep))) - (K / points) Next GoTo avoidcircle circlerun: molecV1 = 0 'for circle minspeed = 1 - (cosrollVA) maxspeed = 1 + (cosrollVA) For K = 0 To 360 'time unit=1 orbit If Option1.Value = True Then molecV1 = tempfac * (Int((830 / ((Rnd * 100) + 1))) - 8) '.00001 / Int(20 * Rnd + 1) 'introducing molecular velocities If Rnd < 0.5 Then sinechange = -molecV1 Else sinechange = molecV1 End If lightspeed(K) = 1 + sinechange + (Cos(pimod * K) * cosrollVA) 'emitted light speed. x c LS(K) = 1 + sinechange - (velratio * Cos(pimod * K) * cosrollVA) position(K) = (lightspeed(K) * (1 + (startpoint / falsep))) - (K / 360) ' pos(K) = (LS(K) * (1 + (startpoint / falsep))) - (K / 360) Next alleccopts Shape1.Visible = False Shape2.Top = 7220 Line7.X1 = 1545: Line7.X2 = 2550: Line7.Y1 = 7230: Line7.Y2 = 7230 Label34.Left = 2520: Label34.Top = 7170 avoidcircle: allcirclabs Shape1.Shape = 3: Line7.Visible = True Shape1.Top = 6930: Shape1.Width = 600: Shape1.Height = 600: Shape1.Left = 1000 Call diagram DrawWidth = 1 Timer1.Enabled = True End If sec = 0 'sec=number of orbits. this is an index for counting time events out: End Function Private Sub Pulsemove_Click() If Timer2.Enabled = False Then Timer2.Enabled = True Else Timer2.Enabled = False End Sub Public Sub Timer1_Timer() 'drawing light fronts Label44.Visible = False '1 time event (p)=1 orbit If ecc = 0 Then 'for circle If position(270) >= stoppoint / n Then 'restarts next part of output back at zero truevel.Cls temppos = position(270) For j = 0 To 360 position(j) = position(j) - temppos pos(j) = pos(j) - temppos Next End If Else 'ellipse If position(1000) >= stoppoint / n Then truevel.Cls temppos = position(1000) For j = 0 To points Step stepsize position(j) = position(j) - temppos pos(j) = pos(j) - temppos Next End If End If If ecc <> 0 Then 'for ellipse For j = 0 To points Step stepsize PSet ((100 * n * position(j)), 300 + (j * 4000 / points)), RGB(255, 0, 0) If Option5.Value = True Then 'with extinction Extspeed(j) = lightspeed(j) - 1 position(j) = position(j) + 1 + (Extspeed(j) * xfactor) Else: position(j) = position(j) + lightspeed(j) End If Next If velratio > 0 Then For j = 0 To points Step stepsize If j <= points / 2 Then PSet ((100 * n * pos(j)), 2300 - (j * 4000 / points)), RGB(0, 255, 255) 'inner star If j > points / 2 Then PSet ((100 * n * pos(j)), 4300 - ((j - (points / 2)) * 4000 / points)), RGB(0, 255, 255) 'inner star If Option5.Value = True Then extS(j) = LS(j) - 1 pos(j) = pos(j) + 1 + (extS(j) * xfactor) Else: pos(j) = pos(j) + LS(j) End If Next End If GoTo gobelow Else 'for circle For j = 0 To 359 Step 1 PSet ((100 * n * position(j)), 300 + (11.11 * j)), RGB(255, 0, 0) If Option5.Value = True Then extcircle(j) = lightspeed(j) - 1 position(j) = position(j) + 1 + (extcircle(j) * xfactor) Else: position(j) = position(j) + lightspeed(j) End If If velratio > 0 Then If j <= 179 Then PSet ((100 * n * pos(j)), 2300 - (11.11 * j)), RGB(0, 255, 255) If j > 179 Then PSet ((100 * n * pos(j)), 4300 - (11.11 * (j - 179))), RGB(0, 255, 255) If Option5.Value = True Then circleS(j) = LS(j) - 1 pos(j) = pos(j) + 1 + (circleS(j) * xfactor) Else: pos(j) = pos(j) + LS(j) End If End If Next End If If ecc = 0 Then 'orbit diagram , bottom left PSet (1260 + (320 * Sin(sec / 5)), 7230 - (320 * Cos(sec / 5))), RGB(0, 0, 0) PSet (1260 + (320 * Sin(sec / 5)), 7230 - (320 * Cos(sec / 5))), RGB(255, 250, 250) End If gobelow: sec = sec + 1 xfactor = xfactor * (xfactor ^ falsep) End Sub Public Sub Command2_Click() 'pause/restart combovalues Label17.Visible = False If Timer1.Enabled = True Then Timer1.Enabled = False Command2.Caption = "CONTINUE" Label19.Visible = True Label29.Visible = True 'show distance Label22.Visible = True LWait.Visible = True If Option5.Value = True Then Combo4.Enabled = False Else Combo4.Enabled = True LWait.Caption = sec + (startpoint / falsep) If ecc = 0 Then If Option5.Value = False Then ninet = startpoint + (minspeed * (falsep * sec)) 'no extinction (circle twenine = startpoint + (maxspeed * (falsep * sec)) Else ninet = (falsep * (sec - 0.75)) 'extinction (circle) twenine = (falsep * (sec - 0.25)) End If Else 'ellipse If Option5.Value = False Then 'no extinction ellipse ninet = startpoint + (falsep * minspeed * (sec - (lastj / points))) twenine = startpoint + (falsep * maxspeed * (sec - (lastk / points))) 'startpoint + (falsep * maxspeed * sec) Else: ninet = (falsep * (sec - 0.5)) 'extinction ELLIPSE twenine = (falsep * sec) End If End If ninet = Format(ninet, "#.00") twenine = Format(twenine, "#.00") Label19.Caption = ninet 'number of L.Y.s Label29.Caption = twenine If ecc = 0 Then 'circle distance labels If Option5.Value = False Then Label29.Top = 330 'no extinction Label19.Top = 2400 Label29.Left = 50 + (100 * n * (position(0) - lightspeed(0))) 'position of label Label19.Left = 50 + (100 * n * (position(180) - lightspeed(180))) Else Label29.Top = 1330 'extinction Label19.Top = 3350 Label29.Left = 50 + (100 * n * (position(90) - 1)) 'position of label Label19.Left = 50 + (100 * n * (position(270) - 1)) End If Else 'ellipse. distance labels If Option5.Value = False Then Label19.Top = 300 + (4000 * lastj / points) Label29.Top = 300 + (4000 * lastk / points) Label19.Left = 50 + (100 * n * (position(lastj) - minspeed)) 'position of distance labels Label29.Left = 50 + (100 * n * (position(lastk) - maxspeed)) 'no extinction Else Label19.Top = 300 + (4000 * ((3 / 4) + (yaw / 360))) Label29.Top = 300 + (4000 * ((1 / 4) + (yaw / 360))) Label19.Left = 50 + (100 * n * position(359 * ((3 / 4) + (yaw / 360)))) 'position of distance labels Label29.Left = 50 + (100 * n * position(359 * ((1 / 4) + (yaw / 360)))) 'with extinction End If End If outerperiod = falsep 'setting import values for brightness curve. outerp = outerperiod * 3.15576 * (10 ^ 7) dtemp = startpoint + (falsep * sec) 'lys dtemp = Format(dtemp, "0000.00") Else 'restart Timer1.Enabled = True Randomize If eccentricity.Text = "" Then 'Stops a crash Label17.Visible = False Label35.Visible = True Else If cut = False Then Label44.Visible = True GoTo outagain End If End If Label35.Visible = False Label22.Visible = False Label19.Left = -1200 Label29.Left = -1200 LWait.Visible = True Label29.Visible = False Label19.Visible = False LWait.Caption = "wait" ecc = eccentricity.Text If ecc <> eccindex Then 'Stops a crash Label44.Visible = True cut = False GoTo outagain End If Label31.Visible = False Command9.Caption = "Instructions" Command2.Caption = "PAUSE" details.Cls eccindex = ecc velratio = Combo11(0) Line1.Visible = False Line9.Visible = False Line12.Visible = False Timer1.Enabled = True Call diagram End If outagain: End Sub '*********************************************************************** '*********************************************************************** '*********************************************************************** Public Sub Frame6_Click() 'start brightness curve Command2.Caption = "CONTINUE" brightness End Sub 'OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 'OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO Public Function brightness() Call combovalues Call diagram Label31.Visible = False Scan.Label16.Visible = False Scan.Label17.Visible = False Scan.Caption = "SCAN OF BRIGHTNESS CURVES VERSUS YAW ANGLE" Command9.Caption = "Instructions" Timer1.Enabled = False If cut = False Then 'Stops a crash Label44.Visible = True GoTo outyetagain End If If ecc <> eccindex Then 'Stops a crash Label44.Visible = True cut = False GoTo outyetagain End If outerperiod = falsep 'setting import values for brightness curve. outerp = outerperiod * 3.15576 * (10 ^ 7) 'seconds PP = outerp / points dtemp = startpoint + (falsep * sec) 'lys dtemp = Format(dtemp, "0000.00") eccindex = ecc Gmod = G * ecc / 0.25 cosroll = Cos(roll) cosrollV = cosroll * vone cosrollVA = cosroll * voneA If velratio <> 0 Then velroot = Combo10 Else velroot = 1 'the slower star must be bigger and therefore brighter If eccentricity.Text = "" Then 'Stops a crash Label35.Visible = True Else Label35.Visible = False 'Output.Cls 'importing values for brightness curve. dtemp = startpoint + (falsep * sec) 'LYs Stardist = dtemp * 9.46 * (10 ^ 12) 'kilometres Erase intensity Erase INTB Erase Tat Erase Tbt BmaxOuter = 0: BminOuter = 100000000 BmaxInner = 0: BminInner = 100000000 BmaxBoth = 0: BminBoth = 100000000 Maginner = 0: magouter = 0: magboth = 0 '************************************************** If Option6.Value = True Then GoTo scanyaw '************************************************** scaling = 12000 * imaghite / points '/ orbs If ecc = 0 Then GoTo circlebright 'for circle v = velocity(0) * c * cosrollV 'main ellipse 'km/sec. 3.15576*10^7 secs per year outerperiod = Format(outerperiod, "0.00000") vone = Format(vone, "0.0000000"): Vtwo = Format(Vtwo, "0.0000000") rone = Format(rone, "00.0000000") 'DETERMINE MAX AND MIN TIMES Tmin = 2 * Stardist / (c + v) 'in seconds. This sets starting angle =0, for min travel time. Inner star is approaching Tmax = 1 'this is the max time (sec) TAmin = 2 * Stardist / (c + v) TAmax = 1 Erase orbtime For A = 0 To orbs 'time lag between subsequent orbits orbtime(A) = A * outerp 'secs Next '************************** '************************** 'main calculation Xlag = Int(points * lag / 360) ' For r = 0 To points - 1 Step 1 'determine T minimum. Y = Cos(Vangle(r) - (pimod * yaw)) If lag <> 0 Then Ylag = Cos(Vangle((r + Xlag) Mod (points - 1)) - (pimod * yaw)) If Option1.Value = True Then molecV1 = tempfac * (Int((830 / ((Rnd * 100) + 1))) - 8) 'molecular velocities - mean velocity= 30*tempfac x c If Rnd < 0.5 Then sinechange = -molecV1 Else sinechange = molecV1 End If Tat(r) = (Stardist / (c * (1 + sinechange + (velocity(r) * Y * cosrollV)))) + (PP * r) 'travel time plus starting lag, one orbiTraveltime(n) = Tat If Tat(r) < Tmin Then Tmin = Tat(r) pointmin = r End If If Tat(r) > Tmax Then Tmax = Tat(r) If lag = 0 Then 'lag = Clag Tbt(r) = (Stardist / (c * (1 + sinechange - (velocity(r) * velratio * Y * cosrollV)))) + (PP * r) 'travel time plus starting lag, outer star. one orbiTraveltime(n) = Tat Else 'Ylag = Cos(Vangle(r + xlag) - (pimod * yaw)) Tbt(r) = (Stardist / (c * (1 + sinechange + (velocity((r + Xlag) Mod (points - 1)) * velratio * Ylag * cosrollV)))) + (PP * r) End If If Tbt(r) < Tmin Then Tmin = Tbt(r) pointmin = r End If If Tbt(r) > Tmax Then Tmax = Tbt(r) 'Maximum light travel time over 1 orbit If orbs = 2 Then If Tat(r) < TAmin Then TAmin = Tat(r) If Tat(r) > TAmax Then TAmax = Tat(r) End If Next If orbs = 2 Then Call Phasecheck GoTo outyetagain Else dfraction = (Tmax + (outerp * (orbs - 1)) - Tmin) / 500 'more than the number of screen pixels. Call Brightcalc If compare.Visible = False Then Call UseOutput Output.Cls 'print combined brightness graph, ellipse, no scan Call Magdiff Call both GoTo outyetagain 'CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC circlebright: Vtwo = voneA * velratio rone = outerperiod * voneA / 2 / pi 'LYs rtwo = rone * velratio cosroll = Cos(roll) cosrollV = cosroll * vone cosrollVA = cosroll * voneA molecV1 = 0 Erase Tat, intensity, Tbt, INTB Tmin = Stardist / (c * (1 + (cosrollVA + 0.0000001))) Tmax = (Stardist / (c - (c * cosrollVA))) + ((orbs) * outerp) dfraction = (Tmax - Tmin) / 500 Erase orbtime For A = 0 To orbs 'time lag between subsequent orbits orbtime(A) = A * outerp 'secs Next For r = 0 To points - 1 If Option1.Value = True Then molecV1 = tempfac * (Int((830 / ((Rnd * 100) + 1))) - 8) 'introducing molecular velocities If Rnd < 0.5 Then sinechange = -molecV1 Else sinechange = molecV1 End If Tat(r) = (Stardist / (c * (1 + sinechange + (Cos(twopi * r / points) * cosrollVA)))) + (PP * r) If Clag = 0 Then Tbt(r) = (Stardist / (c * (1 + sinechange - (Cos(twopi * (r / points)) * velratio * cosrollVA)))) + (PP * r) Else Tbt(r) = (Stardist / (c * (1 + sinechange + (Cos(twopi * ((r / points) + (lag / 360))) * velratio * cosrollVA)))) + (PP * r) End If Next Call Brightcalc If compare.Visible = False Then Call UseOutput Output.Print "Time(max)="; t; "years": Output.Print "Time(min)="; u; "years": Output.Print "Outer Orbit velocity="; voneA; "x c" Output.Print "Inner Orbit velocity="; Vtwo; "x c" End If Output.Cls Call Magdiff Call both End If GoTo outyetagain '********************************************************************************** '***************************************************************************** '******************************************************************** scanyaw: 'repeat but with scan Output.Hide details.Hide Scan.Cls Scan.Command1.Caption = "Wait...." Label31.Visible = False If ecc = 0 Then GoTo circlescan 'for circle v = velocity(0) * c * cosrollV 'km/secs. 3.15576*10^7 secs per year Ylag = Int(points * lag / 360) 'DETERMINE MAX AND MIN TIMES Tmin = Stardist / (c - v) 'in seconds. This sets starting angle =0, for min travel time. Inner star is approaching Tmax = (Stardist / (c - v)) + ((orbs + 1) * outerp) 'this is the max time (sec) 'an orbit has been added to make certain. 'DIV = (outerp) / points 'sets size of time division, in seconds 'divdays = DIV / 86400 'divdays = Format(divdays, "0.000000") dfraction = (Tmax - Tmin) / 500 Erase orbtime For A = 0 To orbs 'time lag between subsequent orbits orbtime(A) = A * outerp 'secs Next '************************** '************************** 'main calculation Erase scanintA, scanintB, Tat, Tbt twovelrat = 0 cosroll = Cos(roll) cosrollV = cosroll * vone cosrollVA = cosroll * voneA For yaw = 0 To 11 'yaw loop yawadd = yaw * 36 piyawb = (pi * yaw / 6) For r = 0 To points 'determine T minimum. Y = Vangle(r) - piyawb Tat(r) = (Stardist / (c * (1 + (velocity(r) * Cos(Y) * cosrollV)))) + (PP * r) 'travel time plus starting lag, one orbiTraveltime(n) = Tat If Tat(r) < Tmin Then Tmin = Tat(r) Next For twovelrat = 0 To 3 'starts velratio loop velratio = twovelrat / 3 Erase Tbt For r = 0 To points Y = Vangle(r) - piyawb If Clag = 0 Then Tbt(r) = (Stardist / (c * (1 - (velocity(r) * velratio * Cos(Y) * cosrollV)))) + (PP * ((r + Ylag) Mod points)) 'travel time plus starting lag, one orbiTraveltime(n) = Tat Else Tbt(r) = (Stardist / (c * (1 + (velocity(r) * velratio * Cos(Y) * cosrollV)))) + (PP * ((r + Ylag) Mod points)) End If If Tbt(r) < Tmin Then Tmin = Tbt(r) For x = 0 To (orbs - 1) 'Y = angle/point in radians If twovelrat <> 0 Then Tb = Tbt(r) + orbtime(x) 'second star of binary pair timefraction = Fix((Tb - Tmin) / dfraction) 'timefraction has maximumm value Dfraction scanintB(twovelrat, yaw, timefraction) = scanintB(twovelrat, yaw, timefraction) + velroot 'sums the number of times each division End If If twovelrat = 0 Then Ta = Tat(r) + orbtime(x) 'adds times for subsequent orbits timefraction = Fix((Ta - Tmin) / dfraction) 'timefraction has approx. maximumm value 3500' scanintA(twovelrat, yaw, timefraction) = scanintA(twovelrat, yaw, timefraction) + 1 'sums the number of times each division End If Next Next Next Next Scan.Show Scan.Label2.Visible = True Scan.Label3.Visible = True Scan.Label13.Visible = False Scan.Label14.Visible = False Scan.Label15.Visible = False Scan.Label16.Visible = False Scan.Label17.Visible = False DrawWidth = 1 Call combined Scan.Command1.Caption = "Back To main Screen" GoTo outyetagain circlescan: 'scan for circle, V2/V1=0.4,0.7,1....luminosity ratio=0.5, 1, 1.5 Vtwo = voneA * velratio rone = outerperiod * voneA / 2 / pi 'LYs rtwo = rone * velratio Output.Cls Output.Print "Outer Orbit velocity="; voneA; "x c" Output.Print "Inner Orbit velocity="; Vtwo; "x c" rone = Format(rone, "0.0E+00") rtwo = Format(rtwo, "0.0E+00") Scan.Label16.Visible = True Scan.Label17.Visible = True Scan.Label13.Visible = True Scan.Label14.Visible = True Scan.Label15.Visible = True Scan.Label16.Visible = False Scan.Label17.Visible = False Scan.Caption = "FOR CIRCULAR ORBITS" Erase scanintA, scanintB, Tat, Tbt, orbtime cosroll = Cos(roll) cosrollV = cosroll * vone cosrollVA = cosroll * voneA For A = 0 To orbs 'time lag between subsequent orbits orbtime(A) = A * outerp 'secs Next Tmin = Stardist / (c * (1 + cosrollVA)) Tmax = (Stardist / (c * (1 - (cosrollVA)))) + ((orbs + 1) * outerp) 'DIV = outerp / points 'sets size of time division, in seconds dfraction = (Tmax - Tmin) / 500 For r = 0 To points Tat(r) = (Stardist / (c * (1 + (Cos(2 * pi * r / points) * cosrollVA)))) + (PP * r) For z = 0 To orbs - 1 Ta = Tat(r) + orbtime(z) 'adds times for subsequent orbits timefraction = Fix((Ta - Tmin) / dfraction) 'timefraction has approx. maximumm value 3500' For f = 0 To 5 scanintA(0, f, timefraction) = scanintA(0, f, timefraction) + 1 'sums the number of times each division Next Next Next '1st for loop starts For velrat = 0 To 5 '2nd for loops start velratio = velrat / 5 Vtwo = vone * velratio Erase Tbt For r = 0 To points Tbt(r) = (Stardist / (c * (1 - (Cos(2 * pi * ((r / points) + (lag / 360))) * velratio * cosrollVA)))) + (PP * r) 'For z = 0 To points For b = 0 To (orbs - 1) 'Y = angle/point in radians Tb = Tbt(r) + orbtime(b) 'second star of binary pair timefraction = Fix((Tb - Tmin) / dfraction) 'timefraction has maximumm value Dfraction For lumrat = 1 To 3 velroot = lumrat / 2 scanintB(lumrat, velrat, timefraction) = scanintB(lumrat, velrat, timefraction) + velroot 'sums the number of times each division Next Next Next Next Scan.Show Scan.Label16.Visible = True Scan.Label17.Visible = True DrawWidth = 1 Call combined Scan.Command1.Caption = "Back To main Screen" outyetagain: End If velratio = 0 End Function Public Function combined() 'linear output, scan Scan.Cls scaling = 20000 * imaghite / points / orbs If ecc <> 0 Then Scan.Label2.Visible = True Scan.Label3.Visible = True Scan.Label1.Visible = True Scan.Label13.Visible = False Scan.Label14.Visible = False Scan.Label15.Visible = False Scan.Label16.Visible = False Scan.Label17.Visible = False Scan.Caption = "COMBINED BRIGHTNESS VERSUS TIME (each cycle repesents one period)" If Option3.Value = True Then For f = 0 To 11 For m = 0 To 405 ''print brightness graph, ellipse, scan Scan.PSet (30 + (m / 2.2), 100 + (39 * f) - (scanintA(0, f, m) * scaling)), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 100 + (39 * f) - ((scanintA(0, f, m) + scanintB(e, f, m)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next Else For f = 0 To 11 For m = 0 To 405 ''print brightness graph, ellipse, scan Scan.PSet (30 + (m / 2.2), 180 + (39 * f) - (50 * Log((scanintA(0, f, m) * scaling) + 1))), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 180 + (39 * f) - (50 * Log(((scanintA(0, f, m) + scanintB(e, f, m))) + 1) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next End If End If If ecc = 0 Then Call labels Scan.Caption = "COMBINED BRIGHTNESS VERSUS TIME (each cycle repesents one period)" If Option3.Value = True Then 'print brightness graph, circle, scan For f = 0 To 5 For m = 0 To 405 Scan.PSet (30 + (m / 2.2), 250 + (50 * f) - (scanintA(0, f, m) * scaling)), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 250 + (50 * f) - ((scanintA(0, f, m) + scanintB(e, f, m)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next Else For f = 0 To 5 'scan, circle, log brightness graphs For m = 0 To 405 Scan.PSet (30 + (m / 2.2), 250 + (50 * f) - (50 * Log((scanintA(0, f, m) * scaling) + 1))), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 250 + (50 * f) - (50 * Log(((scanintA(0, f, m) + scanintB(e, f, m)) + 1)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next End If End If End Function Public Function individual() 'individual star output, scan Scan.Cls Scan.Caption = "INDIVIDUAL STAR CONTRIBUTION" scaling = 20000 * imaghite / points / orbs If ecc <> 0 Then Scan.Label1.Visible = True Scan.Label2.Visible = True Scan.Label3.Visible = True Scan.Label13.Visible = False Scan.Label14.Visible = False Scan.Label15.Visible = False Scan.Label16.Visible = False Scan.Label17.Visible = False If Option3.Value = True Then For f = 0 To 11 For m = 0 To 405 'scan, ellipse, linear brightness graph Scan.PSet (30 + (m / 2.2), 100 + (39 * f) - (scanintA(0, f, m) * scaling)), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 100 + (39 * f) - (scanintB(e, f, m) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next Else For f = 0 To 11 For m = 0 To 405 'scan, ellipse, log brightness graphs Scan.PSet (30 + (m / 2.2), 180 + (39 * f) - (50 * Log((scanintA(0, f, m) + 1)) * scaling)), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 180 + (39 * f) - (50 * Log((scanintB(e, f, m) + 1)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next End If End If If ecc = 0 Then Call labels If Option3.Value = True Then For f = 0 To 5 For m = 0 To 405 'scan, circle, linear brightness graph Scan.PSet (30 + (m / 2.2), 250 + (50 * f) - ((scanintA(0, f, m)) * scaling)), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 250 + (50 * f) - ((scanintB(e, f, m)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next Else For f = 0 To 5 For m = 0 To 405 'scan, circle, log brightness graphs Scan.PSet (30 + (m / 2.2), 250 + (50 * f) - (50 * Log((scanintA(0, f, m) * scaling) + 1))), RGB(84, 255, 0) For e = 1 To 3 Scan.PSet (30 + (190 * e) + (m / 2.2), 250 + (50 * f) - (50 * Log((scanintB(e, f, m) + 1)) * scaling)), RGB(e * 84, 255, 255 - (e * 84)) Next Next Next End If End If End Function Public Function outer() Output.Cls 'no scan, print brightness graph For m = 0 To 500 If Option3.Value = True Then Output.PSet (151 + m, 500 - (intensity(m) * scaling)), RGB(0, 255, 0) Else Output.PSet (151 + m, 500 - (50 * Log((intensity(m) + 1)) * scaling)), RGB(0, 255, 0) End If Next Output.Label19.Caption = magouter End Function Public Function inner() 'no scan, print brightness graph Output.Cls Bmax = 0: Bmin = 10000000 For m = 0 To 500 If Option3.Value = True Then Output.PSet (151 + m, 500 - (INTB(m) * scaling)), RGB(0, 255, 0) Else Output.PSet (151 + m, 500 - (50 * Log((INTB(m) + 1)) * scaling)), RGB(0, 255, 0) End If Next Output.Label19.Caption = Maginner End Function Public Function both() 'no scan, print brightness graph If compare.Visible = False Then Output.Cls Output.Show If ecc <> 0 Then For m = 0 To (points - 1) Step 20 Output.PSet (151 + (m * 62.5 / p), 219 - (0.05 / Gmod * velocity(m) * Cos(Vangle(m) - piyaw))), RGB(0, 0, 255) 'velocity diagram. Output.PSet (151 + ((m + points - 1) * 62.5 / p), 219 - (0.05 / Gmod * velocity(m) * Cos(Vangle(m) - piyaw))), RGB(0, 0, 255) Output.Line (151, 220)-(400, 220) Next End If If Option3.Value = True Then For m = 0 To 500 Output.PSet (151 + m, 500 - ((intensity(m) + INTB(m)) * scaling)), RGB(0, 255, 0) Next Else: For m = 0 To 500 'print combined brightness graph,ellipse Output.PSet (151 + m, 500 - (50 * Log(((intensity(m) + INTB(m)) + 1)) * scaling)), RGB(0, 255, 0) Next End If Output.Label19.Caption = magboth 'Output Else 'Output.Hide truevel.Show compare.Cls For m = 0 To 500 LR = CLR: UD = CUD: SIZE = CSize: WIDE = CWide 'comparing curves compare.PSet (LR + (m / WIDE), UD - (SIZE * Log(((intensity(m) + INTB(m)) + 1)) * scaling)), RGB(0, 155, 0) '+ 2.8 Next Lmags.Caption = magboth End If End Function Public Function alleccopts() vell0 = lightspeed(360) 'velocity labels vell90 = lightspeed(90) vell180 = lightspeed(180) vell270 = lightspeed(270) vell0 = Format(vell0, "0.00000") vell180 = Format(vell180, "0.00000") vell90 = Format(vell90, "0.00000") vell270 = Format(vell270, "0.0000") Label3.Caption = vell0 Label4.Caption = vell0 Label5.Caption = vell180 Label6.Caption = vell90 Label7.Caption = vell270 Label3.Visible = True Label4.Visible = True Label43.Visible = False: Label50.Visible = False rone = falsep * voneA / 2 / pi * 365 'Ldays rone = Format(rone, "00.000000") Label37.Caption = rone Line3.Y1 = 1330: Line3.Y2 = 1330 Line5.Y1 = 3330: Line5.Y2 = 3330 End Function Private Sub Form_MouseDown(Button As Integer, Shift As Integer, x As Single, Y As Single) Line12.Visible = False Line12.Y1 = 0 'yellow line disappears from screen Line12.Y2 = 0 End Sub Private Sub Form_MouseUp(Button As Integer, Shift As Integer, x As Single, Y As Single) Line12.Y1 = 70 Line12.Y2 = 4440 'yellow line appears on screen Line12.X1 = x Line12.X2 = x If Timer1.Enabled = False Then Line12.Visible = True End Sub Public Sub Option3_Click() linearcurve = True End Sub Public Sub Option4_Click() logcurve = True End Sub Public Function allcirclabs() If ecc = 0 Then Label3.Visible = True Label4.Visible = True Label5.Visible = True Label6.Visible = True Label7.Visible = True Else Label3.Visible = False Label4.Visible = False Label5.Visible = False Label6.Visible = False Label7.Visible = False End If End Function Public Function allvalues() Randomize DrawWidth = 1 Label17.Visible = False Label31.Visible = False LF.Visible = False LA.Visible = False LC.Visible = False LB.Visible = False Command9.Caption = "Instructions" Command2.Caption = "PAUSE" Label35.Visible = False Erase lightspeed Erase position Erase orbtime Erase LS Erase pos Erase vell Erase VelB Label17.Visible = False truevel.Cls details.Cls Label19.Left = -1200 LWait.Visible = False Label29.Visible = False Line12.Visible = False Call combovalues eccindex = ecc pi = 3.1415926536 piyaw = pimod * yaw c = 2.997 * (10 ^ 5) 'kilometres per second cosroll = Cos(roll) cosrollV = cosroll * vone cosrollVA = cosroll * voneA stoppoint = 108 Timer1.Interval = 50 LR = -200 UD = 680 SIZE = 200 End Function Public Function coords() 'write some ellipse cordinates to file. Open "c:\ellipseXY" For Output As #1 Print #1, "ecc = "; ecc For n = 1 To 60 Print #1, "X= "; xco(n); " Y= "; yco(n) Next Close #1 End Function Public Function labels() Scan.Label1.Visible = False Scan.Label2.Visible = False Scan.Label3.Visible = False Scan.Label13.Visible = False Scan.Label14.Visible = False Scan.Label15.Visible = False Scan.Label16.Visible = True Scan.Label17.Visible = True End Function Public Function Magdiff() BmaxOuter = 0: BminOuter = 100000000 BmaxInner = 0: BminInner = 100000000 BmaxBoth = 0: BminBoth = 100000000 Maginner = 0: magouter = 0: magboth = 0 For m = 100 To 400 If intensity(m) <> 0 Then If intensity(m) > BmaxOuter Then BmaxOuter = intensity(m) If intensity(m) < BminOuter Then BminOuter = intensity(m) 'establishing brightness ratio End If If INTB(m) <> 0 Then If INTB(m) > BmaxInner Then BmaxInner = INTB(m) If INTB(m) < BminInner Then BminInner = INTB(m) 'establishing brightness ratio End If If intensity(m) <> 0 And INTB(m) <> 0 Then If intensity(m) + INTB(m) > BmaxBoth Then BmaxBoth = intensity(m) + INTB(m) If intensity(m) + INTB(m) < BminBoth Then BminBoth = intensity(m) + INTB(m) End If Next If velratio <> 0 Then Maginner = Log(BmaxInner / BminInner) / 0.921 magouter = Log(BmaxOuter / BminOuter) / 0.921 If velratio = 0 Then magboth = magouter Else magboth = Log(BmaxBoth / BminBoth) / 0.921 magouter = Format(magouter, "0.0000") Maginner = Format(Maginner, "0.0000") magboth = Format(magboth, "0.0000") End Function Private Sub CHide_Click() If compare.Visible = False Then compare.Visible = True CHide.Caption = "HIDE CURVE" If compare.Picture = Empty Then compare.Picture = LoadPicture("c:\testcurves\curve1.jpg") Else compare.Visible = False CHide.Caption = "SHOW CURVE" End If End Sub Public Sub Cgetvalues_Click() Call Getvalues CLR.Text = LR: CUD.Text = UD: CWide.Text = WIDE: CSize.Text = SIZE End Sub Public Sub compare_Click() compare.Cls CHide.Caption = "HIDE CURVE" LR = CLR: UD = CUD: SIZE = CSize: WIDE = CWide If CUD <> "" And CLR <> "" And CSize > "" And CWide > "" Then If UD <> 0 And CLR <> 0 And CSize > 0 And CWide > 0 Then Call keepvalues Call both Else Call Getvalues End If End If End Sub Public Function keepvalues() Open "c:\oldvalues.txt" For Output As #1 Write #1, LR, UD, WIDE, SIZE Close #1 End Function Public Function Getvalues() Open "c:\oldvalues.txt" For Input As #1 Input #1, LR, UD, WIDE, SIZE Close #1 End Function Public Function UseOutput() Output.Cls 'Output.Print "Time Interval for Monitoring Intensity ="; divdays; "days" t = Tmax / 3.15576 / (10 ^ 7) u = Tmin / 3.15576 / (10 ^ 7) t = Format(t, "0.00E+00"): u = Format(u, "0.00E+00") Output.Print "Time(max)="; t; "years": Output.Print "Time(min)="; u; "years": Output.Print "pointmin=", pointmin: dtemp = Format(dtemp, "0.00E+00") Output.Text1.Text = dtemp: Output.Text2.Text = rone: Output.Text3.Text = outerperiod Output.Text4.Text = orbs: Output.Text5.Text = points: Output.Text8.Text = rtwo: Output.Text9 = voneA: Output.Text10 = velratio Output.Text11.Text = yaw: Output.Text12.Text = ecc: Output.Text13 = Combo5.Text Output.Show End Function Public Function Phasecheck() Output.Cls Output.Show dfraction = (TAmax + (outerp) - TAmin) / 500 'more than the number of screen pixels. For r = 0 To points - 1 For x = 0 To 1 Ta = Tat(r) + orbtime(x) 'adds times for subsequent orbits timefraction = Fix((Ta - TAmin) / dfraction) intensity(timefraction) = intensity(timefraction) + 1 'sums the number of times each division Next Next For m = 0 To 500 If Option3.Value = True Then Output.PSet (151 + m, 500 - (intensity(m) * scaling)), RGB(0, 255, 0) Else Output.PSet (151 + m, 500 - (50 * Log((intensity(m) + 1)) * scaling)), RGB(0, 255, 0) End If Next For m = 0 To (points - 1) Step 20 Output.PSet (151 + (m * 125 / p), 219 - (0.05 / Gmod * velocity(m) * Cos(Vangle(m) - piyaw))), RGB(0, 0, 255) 'velocity diagram. Output.PSet (151 + ((m + points - 1) * 125 / p), 219 - (0.05 / Gmod * velocity(m) * Cos(Vangle(m) - piyaw))), RGB(0, 0, 255) Output.Line (150, 220)-(650, 220) Next Output.Label19.Caption = magouter End Function Public Function XYCoords() If p Mod (400) = 0 Then 'print out sample coordinates. Xshort = 100 * Xtemp Yshort = 100 * Ytemp Xshort = Format(Xshort, "000.0000") Yshort = Format(Yshort, "000.0000") xco(p / 400) = Xshort yco(p / 400) = Yshort End If End Function Public Function Brightcalc() For r = 0 To points - 1 'Step 50 For x = 0 To (orbs - 1) Ta = Tat(r) + orbtime(x) 'adds times for subsequent orbits timefraction = Fix((Ta - Tmin) / dfraction) If timefraction >= 0 And timefraction < 501 Then intensity(timefraction) = intensity(timefraction) + 1 'sums the number of times each division If velratio <> 0 Then 'second star of binary pair Tb = Tbt(r) + orbtime(x) timefraction = Fix((Tb - Tmin) / dfraction) If timefraction >= 0 And timefraction < 501 Then INTB(timefraction) = INTB(timefraction) + velroot 'sums the number of times each division End If Next Next End Function Public Function Pulsepos() Erase lightspeed, position Pulsemove.DrawWidth = 1 'orbs = 3 For K = 0 To 12000 lightspeed(K) = (Sin(twopi * K / 12000) * voneA) 'emitted light speed. x c position(K) = (lightspeed(K) * (1 + (startpoint / falsep))) - (K / 12000) L = Int(K / 30) hite(L) = hite(L) + 1 Next End Function Public Sub Timer2_Timer() For K = 0 To 12000 position(K) = position(K) + lightspeed(K) L = Int(K / 30) hite(L) = (hite(L) + 1) / hite(L) Next For m = 0 To 400 Pulsemove.Line (410 - m, 300 - hite(m))-(410 - m, 300), RGB(Int(127 * (1 + Sin(pi * m / 66))), Int(127 * (1 + Cos(pi * m / 33))), Int(127 * (1 - Sin(pi * m / 66)))) Next End Sub > >> This star, if it exists, plays little part in the basic observed brightness >> curve. > >Right. >C is not eclipsing A and B, and it is too far >away to have a great effect. >It does have a small observable effect on A and B, though. > > > > Algol is obviously not an eclipsing binary. > >Quite. >The light curve is exactly as it should be if it were >an eclipsing binary - with the second minimum and all - >So according to Wilsonian logic, it obviously isn't. How anyone can find evidence of a secondary eclipse from the rough data available is beyond me. > >Paul HW. www.users.bigpond.com/hewn/index.htm see: www.users.bigpond.com/hewn/variablestars.exe "Sometimes I feel like a complete failure. The most useful thing I have ever done is prove Einstein wrong".
From: George Dishman on 4 Oct 2005 18:11 "Henri Wilson" <H@..> wrote in message news:64u5k1p9ehdpu5t9bppems72m33cke7hsg(a)4ax.com... > On Tue, 4 Oct 2005 10:45:06 +0100, "George Dishman" > <george(a)briar.demon.co.uk> > wrote: > >> >>"Henri Wilson" <H@..> wrote in message >>news:pth4k1tckn79et44jthevig0h9lcohvih5(a)4ax.com... > >>>>If your theory isn't capable of making a >>>>prediction for this test then it isn't a >>>>theory, and if its prediction is wrong, as >>>>Ritz's was, then the theory is falsified. >>> >>> George, I told you. Each component is moving noramally in the frame of >>> the >>> next >>> component. >> >>And I told you, and you agreed, that the light >>is not moving normally to the mirrors, but that >>is beside the point. > > George, the 'c' component is not....but the 'v' component IS. And as I told you the result is the vector sum for Ritzian theory which produces a null prediction. If you disagree, show your prediction and how you get it. >>Henri, what is needed is >>for you to take the experimental setup of Sagnac >>and derive a prediction using your theory. Show >>your working so that other people can understand >>how to apply your equations for their own >>experiments. That is what a theory is, something >>that anyone can use as a predictive tool. > > I have already told you why I don't believe the sagnac effect is related > to > light speed. I don't care what you believe, if you want to claim you have a theory, you have to make it available for people to use. That's what a theory is! >>My expectation is that applying your BaT to the >>Sagnac will give the same null prediction that >>using Ritzian theory gives, but you have never >>said how what you call "BaT" differs from Ritz. > > It is the same. > It is easier to write 'BaT'. Sure, but knowing they are the same means now I don't have to keep qualifying my answers to say Ritz is falsified but BaT might not be if it differs. They are both falsified by Sagnac. George
From: Henri Wilson on 4 Oct 2005 18:52 On 4 Oct 2005 04:17:04 -0700, "Paul B. Andersen" <paul.b.andersen(a)hia.no> wrote: >Henri Wilson wrote: >> On 3 Oct 2005 02:06:27 -0700, "Paul B. Andersen" <paul.b.andersen(a)hia.no> >> wrote: >> >> >Henri Wilson wrote: >> > >> >> >> >> A, I have supported you on this. In Algol's case, the WCH happens >> >> to be the >> >> large planet 'Androcles'. Do you have any objections to that? >> >> It is likely that flares occuring on the main star are reflected from >> >> 'Androcles' and mistaken for flares on IT. >> > >> >Very likely, indeed. :-) >> > >> >There IS no limit to the stupidities you can utter, is there? :-) >> >You have no idea of what a flare is, do you? :-) >> >> I should do. I used to take lovely photos of them. > >Sorry, Henri. >I don't believe you. >Not very smart to bluff about something you don't >understand the consequence of. Paul, whether or not you believe me is of little concern. While studying, I once worked as a tech assistant in the solar physics section of CSIRO here in OZ. I helped make, assemble and use a very neat solar telescope. It took very nice pictures of the sun's surface including soem very spectacular flares from the edge. I'll bet you have never seen moving shots of a flare. >> Flare stars are usually pretty dark red but every now and then emit a very >> large flare. >> I would categorize them as intrinsically variable. >> >> So a large flare is emitted from the star and is reflected off >> the large planet >> 'Androcles', which has a much higher radial velocity than the star. > >Of course, Henri. :-) >Seen from Sirius, a Solar flare will obviously >easily be mistaken for a flare at Jupiter. Paul, flare stars are usually pretty dull red. They have a habit of spontaneously becoming very bright presumably as a huge flare erupts. If a planet as large as jupiter was orbiting at close radius, it would also become very bright at the time. The flare could easily be on the far side of the star and so the reflected light would be more prominant to us than the red star. Doppler shift of the flare's light from the star could easily and wrongly be taken as radial velocity of the star itself. .....but all this is probably too hard for you.... >> >> However I might be inclined to agree with the tusselad that a small third >> >> object does orbit the main Algol star. >> > >> >The "third object" is not a particularly small star. >> >It is an A5 V. That means it is bigger and brighter than the Sun. >> >But it is smaller than Algol A and B. >> >> It IS small and there is NO Algol B. > >Of course there is no K2 star, Henri. >The observed K2 spectrum is predicted by the BaT, isn't it? >That's what the BaT predicts an Anrocles planets radiates, isn't it? >And the K2 spectrum is observed by Einstainania anyway, >and can not be trusted. Let me phrase that correctly. The observed K2 spectrum is assumed from information that arrives here as a consequence of the BaT. The K2 spectrum is INTERPRETED according to Einsteiniana and certainly cannot be trusted.. >> >> On the published curve, there is a small departure from our predicted curves >> >> that might be associated with an object orbiting with the same period but >> >> lagging in phase behind the main star. >> > >> > >> >The orbital period of Algol C is 681 days, that is 235 times >> >the period of the inner binary. >> >> I wont argue. >> >> > >> >The "small departure" from your predicted curves, >> >is the second minimum when Algol A eclipses Algol B. >> >> Nah! There is NO eclipse. > >Of course there isn't, Henri. :-) >No K2 star - no eclipse. Correct. Oh, there could be. The brightness curves would be just about the same. The main difference is that the one involving an eclipse is dead flat between the troughs. ..unlike Algol's , which has a distinct concavity....just as the BaT predicts. >> I accept that the BaT curves for single stars are very hard >> to distinguish from >> the occasional GENUINE eclipsing binary. > >But Algol isn't a GENUIN eclipsing binary, is it? :-) >So the second minimum is obviously not really there. >It is observed by Einsteiania only. See the curve for yourself. www.users.bigpond.com/hewn/group1.jpg > >> >> > >> >> This star, if it exists, plays little part in the basic observed brightness >> >> curve. >> > >> > >> >Right. >> >C is not eclipsing A and B, and it is too far >> >away to have a great effect. >> >It does have a small observable effect on A and B, though. >> >> Only if you are addicted to Einsteiniana. > >Of course, Henri. >The BaT can obviously explain the observed effects >which the stupid astronomers assign to Algol C. Algol C could easily exist. I never said it didn't. It is more than likely that Algol B doesn't exist. > >> >> Algol is obviously not an eclipsing binary. >> > >> >Quite. >> >The light curve is exactly as it should be if it were >> >an eclipsing binary - with the second minimum and all - >> >So according to Wilsonian logic, it obviously isn't. >> >> see: www.users.bigpond.com/hewn/group1.jpg > >No Algol there. Try U Cep. > >> The Algol curve is identical to the BaT prediction. > >Except for the second minimum, of course. > >But the second minimum is obviously not there, >it is only observed. >And as we know, observations can never be trusted. A long dip is observed. It is wrongly interpreted as a secondary eclipse in an attempt to make observations match Einsteiniana.. > >It's not like it is invented by Henri Wilson. > >> >> Androcles has produced the same curve. > >Indeed he have. > >Seriously: >The fact is of course that you and Androcles >have falsified the ballistic theory by demonstrating >that you have to use fantasy parameters which are wildly >wrong to mimic the light curve. >And even then, you miss the second minimum in the light curve. We do not miss it. There is a continuous dip in the curve between the main dips. Because it isn't dead flat, as required by Einsteiniana, astronomers concocted the idea of a secondary eclipse. > >To claim that all the observations made by Algol >are wrong is just too stupid. >Algol is observed at just about all wavelengths, >And every observation fits the fact that Algol is >a B8 star and a K2 star eclipsing each other. >And as always most of the information is found >in the spectra. Illusions can be very strange Paul. > >For example, the K2 star is very bright in the X-rays, >while the B8 is dark. The X-ray "light curve" has a deep >minimum when the visual light curve has its secondary >minimum. The reason is obvious. >http://www.edpsciences.org/articles/aa/abs/2003/49/aa0057/aa0057.html "No eclipse of the quiescent X-ray emission is observed leading us to the conclusion that the overall coronal filling factor of Algol B is small." This beautifully epitomizes how astronomy has been completely stuffed by Einteiniana. >You can retrieve data from the XMM-Newton satellite archive >and see for yourself. >http://xmm.esac.esa.int/external/xmm_data_acc/xsa/index.shtml > > >The BaT can't even get the light curve right, >not to mention the spectra. Different parts of the Spectra comes from different radii, which have different radial velocities and therefore different brightness characteristics. That is why IR brightness usually varies a lot less than visible. > >Paul HW. www.users.bigpond.com/hewn/index.htm see: www.users.bigpond.com/hewn/variablestars.exe "Sometimes I feel like a complete failure. The most useful thing I have ever done is prove Einstein wrong".
From: Henri Wilson on 4 Oct 2005 18:58
On Tue, 4 Oct 2005 23:06:51 +0100, "George Dishman" <george(a)briar.demon.co.uk> wrote: > >"Henri Wilson" <H@..> wrote in message >news:0tt5k1tb7afbn1f7d7avltq14eof0kujqe(a)4ax.com... >> On Tue, 4 Oct 2005 11:27:40 +0100, "George Dishman" >> They eventually learnt to allow for a CONSTANT light speed. > >No, they just got better instruments that were able >to measure the very low rotation rates involved. >Look up the ICRF. > >> The still haven't woken up to the fact that all starlight isn't >> miraculously >> adjusted (by the fairies) to travel to little planet Earth at exactly 'c'. > >Flaunting your inability to comprehend something >as simple as SR isn't helping your case. George, the SRian postulate that all starlight travels to little planet Earth at the same speed 'c' is a direct consequence of combining LET with christianity. > >George > HW. www.users.bigpond.com/hewn/index.htm see: www.users.bigpond.com/hewn/variablestars.exe "Sometimes I feel like a complete failure. The most useful thing I have ever done is prove Einstein wrong". |