Images of 黄経
tŞEÄEHŞE~Eń\lßC@@@ dim solam(50,50),haru(9),natsu(9),aki(9),fuyu(9), yy,dddim xxx(9)w=weekday(now)document.write "{úÍ "& year(now) &" N"& month(now) &""& day(now) &"ú "& weekdayname(w)& " "document.write "@@MJD= "& mjd(year(now),month(now),day(now))&"
"Function mjd(yy,mm,dd) if mm<3 then yy=yy-1 mm=mm+12 end if mjd=int(365.25*yy)+int(yy/400)-int(yy/100)+int(30.59*(mm-2))+dd-678912end function NtŞ´žJ§ÄŹäíÄŹĺ§HIHŞŚI~§~Źáĺá~ŹŚĺŚ§tJ
[ĺžzŠo ß@@@@@@@tŞÍÄÍHŞÍ~Í Dim a(50) Dim ss(50)dim h(30)sub window_onload() fm.t0.value=year(now) com()end subsub com() for k=1 to 9 yy =fm.t0.value y0=fm.t0.value ' xx=fm.t9.value if k<5 then d=17 m=k*3 elseif k<9 then d=3 m=3*K-13 else d=3 m=3+xx/30 end if for i=0 to 9 dd=i+d c=(mjd(yy,m,dd-9/24)-mjd(1899,13,1.5))/36525 SML=280.6824+36000.769325*c+7.22222*0.0001*c*c SPL=281.2206+1.717697*c+4.83333*0.0001*c*c+2.77777*0.000001*c*c*c g=SML-SPL gr=g*3.141592/180 lamda=SML+1.91946*sin(gr)+0.0200939*sin(2*gr)-4.78889*0.001*sin(gr)*c-1.44444*0.00001*sin(gr)*c*c solam(k,i)=lamda-int(lamda/360)*360 if solam(k,i)>350 then solam(k,i)=solam(k,i)-360 next next for i=0 to 8 haru(i)=(solam(1,i))*(solam(1,i+1)) natsu(i)=(solam(2,i)-90)*(solam(2,i+1)-90) aki(i)=(solam(3,i)-180)*(solam(3,i+1)-180) fuyu(i)=(solam(4,i)-270)*(solam(4,i+1)-270) if haru(i)<0 then fm.t1.value ="3 "&i+17&"ú " & weekdayname(weekday(dateserial(y0,3,i+17))) if natsu(i)<0 then fm.t2.value ="6 "&i+17&"ú "&weekdayname(weekday(dateserial(y0,6,i+17))) if aki(i)<0 then fm.t3.value ="9 "&i+17&"ú "&weekdayname(weekday(dateserial(y0,9,i+17))) if fuyu(i)<0 then fm.t4.value="12 "&i+17&"ú "&weekdayname(weekday(dateserial(y0,12,i+17))) nextend subfunction clock()' fm.tt.value=time st=settimeout("clock()",1000)end function'===================================Sub sel_onClick s=fm.sel.Value if s mod 2 =0 then chuset="" else chuset="ß" end if fm.ttt.value=15*s & "" & chuset a(s)=15*s If s >= 19 And s <= 23 Then aq = 15*s - 360 Else aq = 15*s End If yy =fm.t0.value y0=fm.t0.value mjd1 = mjd(yy, 3, 21) + aq mjd2 = mjd1 - mjd(1899, 12, 30) For i = 0 To 20 dd = Day(mjd2) + i - 5 c = (mjd(Year(mjd2), Month(mjd2), dd - 9 / 24) - mjd(1900, 1, 1.5)) / 36525 SML = 280.6824 + 36000.769325 * c + 7.22222 * 0.0001 * c * c SPL = 281.2206 + 1.717697 * c + 4.83333 * 0.0001 * c^2 + 2.77777 * 0.000001*c^3 g = SML - SPL gr = g * 3.141592 / 180 lamda = SML + 1.91946 * Sin(gr) + 0.0200939 * Sin(2 * gr) - 4.78889 * 0.001 * Sin(gr) * c - 1.44444 * 0.00001 * Sin(gr) * c * c solam(s, i) = lamda - Int(lamda / 360) * 360 If solam(s, i) > 359 Then solam(s, i) = solam(s, i) - 360 End If Next For i = 0 To 15 ss(i) = (solam(s, i) - a(s)) * (solam(s, i + 1) - a(s)) If ss(i) < 0 Then fm.hizuke.value = Month(mjd2) & " " & Day(mjd2) + i - 5 & "ú " & WeekdayName(Weekday(DateSerial(Year(mjd2), Month(mjd2), Day(mjd2) + i - 5))) h(s) = Month(mjd2) & " " & dd & "ú " mjdd = mjd(Year(mjd2), Month(mjd2), Day(mjd2) + i - 5) End If Next th=3.141592/180*lamda ' g2.style.postop=406+250*sin(th) ' g2.style.posleft=648-250*cos(th)End sub'=======================================
