Option DEFAULT INTEGER
Option EXPLICIT
MODE 3
CLS

Sub ShowCurve(x(), y(), n, colour, title$)
   CLS
   Text 10, 10, title$, , , 1
   Bezier x(), y(), n, colour
   Pause 3000
End Sub
' One petal, closed
Dim fx%(3) = (320, 380, 450, 320)
Dim fy%(3) = (240, 100, 100, 240)
Math SCALE fx%(),MM.HRES/640,fx%()
Math SCALE fy%(),MM.VRES/480,fy%()

'-------------------------------------------------------
' 2. SPIRAL (multi-point Bzier path)
'    (approximation using many control points)
'-------------------------------------------------------
Dim sx(10) = (320, 380, 420, 450, 430, 360, 260, 200, 180, 230, 320)
Dim sy(10) = (240, 230, 200, 140, 90, 70, 100, 180, 260, 300, 320)
Math SCALE sx(),MM.HRES/640,sx()
Math SCALE sy(),MM.VRES/480,sy()

'-------------------------------------------------------
' 3. WAVY RIBBON (smooth horizontal wave)
'-------------------------------------------------------
Dim wx(6) = (20, 120, 220, 320, 420, 520, 620)
Dim wy(6) = (240, 120, 360, 120, 360, 120, 240)
Math SCALE wx(),MM.HRES/640,wx()
Math SCALE wy(),MM.VRES/480,wy()

'-------------------------------------------------------
' 4. SIGNATURE-STYLE CURVE (fancy stroke)
'-------------------------------------------------------
Dim sigx(7) = (40, 120, 200, 320, 420, 500, 600, 580)
Dim sigy(7) = (300, 260, 320, 200, 260, 180, 260, 300)
Math SCALE sigx(),MM.HRES/640,sigx()
Math SCALE sigy(),MM.VRES/480,sigy()

'-------------------------------------------------------
' 5. Previous demo shapes (S-curve, loops, etc.)
'-------------------------------------------------------
Dim x1(3) = (50, 200, 440, 590)
Dim y1(3) = (400, 50, 430, 80)
Math SCALE x1(),MM.HRES/640,x1()
Math SCALE y1(),MM.VRES/480,y1()

Dim x2(5) = (100, 300, 300, 200, 440, 540)
Dim y2(5) = (240, 480, 0,   400, 80,  240)
Math SCALE x2(),MM.HRES/640,x2()
Math SCALE y2(),MM.VRES/480,y2()

Dim x3(12) = (320,390,441,460,441,390,320,250,199,180,199,250,320)
Dim y3(12) = (240,309,309,240,171,171,240,309,309,240,171,171,240)
Math SCALE x3(),MM.HRES/640,x3()
Math SCALE y3(),MM.VRES/480,y3()

Dim x4(5) = (320, 600, 40, 560, 100, 320)
Dim y4(5) = (450, 380, 140, 200, 60, 30)
Math SCALE x4(),MM.HRES/640,x4()
Math SCALE y4(),MM.VRES/480,y4()

Dim hx1(3) = (320, 190, 190, 320)
Dim hy1(3) = (350, 200, 100, 160)
Math SCALE hx1(),MM.HRES/640,hx1()
Math SCALE hy1(),MM.VRES/480,hy1()

Dim hx2(3) = (320, 450, 450, 320)
Dim hy2(3) = (160, 100, 200, 350)
Math SCALE hx2(),MM.HRES/640,hx2()
Math SCALE hy2(),MM.VRES/480,hy2()

Dim x6(4) = (50, 300, 200, 580, 400)
Dim y6(4) = (420,350, 50, 120, 300)
Math SCALE x6(),MM.HRES/640,x6()
Math SCALE y6(),MM.VRES/480,y6()

'-------------------------------------------------------
'                     MAIN LOOP
'-------------------------------------------------------
Do
   DrawFlower()

   ShowCurve sx(), sy(), 11, RGB(YELLOW),  "Spiral Sweep (11 points)"
   ShowCurve wx(), wy(), 7,  RGB(CYAN),    "Wavy Ribbon Pattern"
   ShowCurve sigx(), sigy(), 8, RGB(WHITE),"Signature-Style Stroke"

   ShowCurve x1(), y1(), 4, RGB(YELLOW),  "Smooth S-Curve Ribbon"
   ShowCurve x2(), y2(), 6, RGB(CYAN),    "Loop-de-Loop"
   ShowCurve x3(), y3(), 13, RGB(RED),     "Figure-Eight"
   ShowCurve x4(), y4(), 6, RGB(GREEN),   "Spiral-Sweep Path"

   CLS
   Text 10,10,"Heart Shape",,,1
   Bezier hx1(), hy1(), 4, RGB(RED)
   Bezier hx2(), hy2(), 4, RGB(RED)
   Fill MM.HRES\2,MM.VRES\2,RGB(red)
   Pause 3000

   ShowCurve x6(), y6(), 5, RGB(WHITE), "Abstract Dragon Stroke"

Loop
Sub DrawFlower()
   Local i, p
   Local float angle, r
   Local dx, dy, x, y
   Local rx!(3), ry!(3)   ' FLOATING ARRAYS
   Local xi%(3), yi%(3)

   CLS
   Text 10,10,"Flower (12 petals)",,,1

   For i = 0 To 11
       angle = i * 30
       r = angle * Pi / 180

       For p = 0 To 3
           dx = fx%(p) - MM.HRES\2
           dy = fy%(p) - MM.VRES\2

           rx!(p) = MM.HRES\2 + dx * Cos(r) - dy * Sin(r)
           ry!(p) = MM.VRES\2 + dx * Sin(r) + dy * Cos(r)
       Next p

       ' Convert float
       For p = 0 To 3
           xi%(p) = Cint(rx!(p))
           yi%(p) = Cint(ry!(p))
       Next p

       Bezier xi%(), yi%(), 4, RGB(MAGENTA)
       x=0
       y=0
       For p= 0 To 3
         Inc x,xi%(p)
         Inc y,yi%(p)
       Next

       Fill x\4,y\4,RGB(yellow)
   Next i

   Pause 3000
End Sub

' Demo: unified SINC smoothing vs interpolation
OPTION EXPLICIT
mode 2
CLS

CONST N = 40          ' input points
CONST M = 120         ' resampled output points (m != n triggers interpolation)
CONST WIN = 17        ' window size (odd)
CONST FREQ = 0.12     ' normalized cutoff (0<freq<=0.5)

DIM x_in(N), y_in(N)
DIM x_s(N), y_s(N)       ' smoothed outputs (m omitted => smoothing)
DIM x_out(M), y_out(M)   ' resampled outputs (m provided => interpolation)
DIM i%, xi%, yi%

' Build a test waveform that fits the screen
FOR i% = 0 TO N - 1
 x_in(i%+1) = (MM.HRES - 1) * i% / (N - 1)
 y_in(i%+1) = (MM.VRES\2) + (MM.VRES\3) * SIN(2 * PI * i% / 12) + 0.15 * MM.VRES * SIN(2 * PI * i% / 5)
 IF y_in(i%+1) < 0 THEN y_in(i%+1) = 0
 IF y_in(i%+1) > MM.VRES - 1 THEN y_in(i%+1) = MM.VRES - 1
NEXT

' Smoothing (m omitted => m=n internally)
MATH SINC x_in(), y_in(), N, WIN, FREQ, x_s(), y_s()

' Interpolation/resampling (m provided)
MATH SINC x_in(), y_in(), N, M, WIN, FREQ, x_out(), y_out()

' Clear and draw axes
CLS
LINE 0, MM.VRES\2, MM.HRES-1, MM.VRES\2, 1, RGB(GREY)
LINE 0, 0, 0, MM.VRES-1, 1, RGB(GREY)

' Plot original points in red
FOR i% = 1 TO N
 xi% = INT(x_in(i%))
 yi% = INT(y_in(i%))
 IF xi% >= 0 AND xi% < MM.HRES AND yi% >= 0 AND yi% < MM.VRES THEN PIXEL xi%, yi%, RGB(RED)
NEXT

' Plot smoothed points in yellow
FOR i% = 1 TO N
 xi% = INT(x_s(i%))
 yi% = INT(y_s(i%))
 IF xi% >= 0 AND xi% < MM.HRES AND yi% >= 0 AND yi% < MM.VRES THEN PIXEL xi%, yi%, RGB(YELLOW)
NEXT

' Plot resampled points in green, slightly offset upward for visibility
FOR i% = 1 TO M
 xi% = INT(x_out(i%))
 yi% = INT(y_out(i%)) - 2  ' small offset for visibility
 IF yi% < 0 THEN yi% = 0
 IF yi% >= MM.VRES THEN yi% = MM.VRES - 1
 IF xi% >= 0 AND xi% < MM.HRES THEN PIXEL xi%, yi%, RGB(GREEN)
NEXT

>location "18/12/2025 11:03:00",44.8488,-74.2937
>astro "Vega",a,az,ra,dec
Error : Argument count
>

location "18/12/2025 11:03:00",44.8488,-74.2937
astro "Vega",a,az,ra,dec
Error : Argument count