Les N insectes amoureux

' ****************************************************
' *                                                                                       
' *          Les N insectes amoureux                                       
' *                par Papydall                                               
' *                                                 
' ****************************************************
dim titre$ : titre$ = "Les N insectes amoureux"
width 0,700 : height 0,700
caption 0,"Papydall vous présente        " + titre$
dim N,DC,NbPOL
' ******************************************************************************
' Ces constantes peuvent être modifiées
N = 10      : ' Nombre d'insectes
DC = 12      : ' diviseur du coté
NbPol = 200  : ' Nombre de polygones

' ******************************************************************************
dim X(N+1),Y(N+1),c,pi,angle,i,j,x1,x2,y1,y2
dim lx,ly,xmin,xmax,ymin,ymax
pi = 4*atn(1)
lx = 500 : ly = 500                          : ' Dimension de l'écran graphique
xmin = -1 : xmax = 1 : ymin = -1 : ymax = 1 : ' Coordonnées des bords de l'écran
label InitTabCoord, tracer
2d_fill_color 0,0,0 : 2d_rectangle 50,50,650,650
gosub InitTabCoord : gosub tracer
end

' ******************************************************************************
InitTabCoord:
c = 0
for angle = 0 to 2*pi step 2*pi/n
    c = c + 1
    ' Coordonnées des sommets du polygone
    x(c) = cos(angle) : y(c) = sin(angle)
next angle
x(n+1) = x(1) : y(n+1) = y(1) : ' Pour faciliter les calculs et le tracé
return
' ******************************************************************************
tracer:
for i = 1 to NbPol
    x1 = lx*(x(1)-xmin)/(xmax-xmin)
    y1 = ly*(y(1)-ymin)/(ymax-ymin)
    for j = 2 to n+1
        x2 = lx*(x(j)-xmin)/(xmax-xmin) : y2 = ly*(y(j)-ymin)/(ymax-ymin)
        2d_pen_color 0,255,255
        2d_line x1+100,y1+100,x2+100,y2+100
        x1 = x2 : y1 = y2
    next j
    ' Calcul des coordonnées du polygone suivant
    for j = 1 to n
        x(j) = x(j) + (x(j+1) - x(j))/DC
        y(j) = y(j) + (y(j+1) - y(j))/DC
    next j
    x(n+1) = x(1) : y(n+1) = y(1)
next i
return

' ****************************************************
' *
' *          Les N insectes amoureux
' *                par Papydall
' *
' ****************************************************
dim titre$ : titre$ = "Les N insectes amoureux"
width 0,700 : height 0,700 : caption 0,"Papydall vous présente        " + titre$
dim N, DC,NbPOL
' ******************************************************************************
' Ces constantes peuvent être modifiées
N = 4        : ' Nombre d'insectes
DC = 25      : ' diviseur du coté
NbPol = 58  : ' Nombre de polygones

' ******************************************************************************
dim X(N+1),Y(N+1),c,pi,angle,i,j,x1,x2,y1,y2
dim lx,ly,xmin,xmax,ymin,ymax
pi = 4*atn(1)
lx = 500 : ly = 500                          : ' Dimension de l'écran graphique
xmin = -1 : xmax = 1 : ymin = -1 : ymax = 1 : ' Coordonnées des bords de l'écran
label InitTabCoord, tracer
2d_fill_color 0,0,0 : 2d_rectangle 50,50,650,650
gosub InitTabCoord : gosub tracer
end

' ******************************************************************************
InitTabCoord:
c = 0
for angle = 0 to 2*pi step 2*pi/n
    c = c + 1
    ' Coordonnées des sommets du polygone
    x(c) = cos(angle) : y(c) = sin(angle)
next angle
x(n+1) = x(1) : y(n+1) = y(1) : ' Pour faciliter les calculs et le tracé
return
' ******************************************************************************

tracer:
for i = 1 to NbPol
    x1 = lx*(x(1)-xmin)/(xmax-xmin)
    y1 = ly*(y(1)-ymin)/(ymax-ymin)
    for j = 2 to n+1
        x2 = lx*(x(j)-xmin)/(xmax-xmin) : y2 = ly*(y(j)-ymin)/(ymax-ymin)
        2d_pen_color y1/3,255-i,x1/2
        2d_line x1+100,y1+100,x2+100,y2+100
        x1 = x2 : y1 = y2
    next j
    ' Calcul des coordonnées du polygone suivant
    for j = 1 to n
        x(j) = x(j) + (x(j+1) - x(j))/DC
        y(j) = y(j) + (y(j+1) - y(j))/DC
    next j
    x(n+1) = x(1) : y(n+1) = y(1)
next i
return

dim x0,y0,r,Px,Py,Px2,Py2,j%,n%,Pi,a,v
width 0,800
height 0,800
rem coordonnées du centre
x0=400
Y0=400
color 0,0,0,0
rem rayon
r=380

rem nombre de côtés
n%=7
for v=1 to 140
r=r-3
rem angle de début
a=a+50/57.29

Pi = 3.1415926535898

Px = r * Cos(a) + x0
Py = r * Sin(a) + y0

2D_poly_from Px,Py

For j% = 1 To n% :wait 1
    Px2 = r*Cos(a + 2*Pi*j%/n%) + x0
    Py2 = r*Sin(a + 2*Pi*j%/n%) + y0
    2D_poly_to Px2, Py2
    Px = Px2
    Py = Py2
Next j%
2d_pen_color rnd( 255),rnd(125),rnd(240)
next v
end

dim x0,y0,r,Px,Py,Px2,Py2,j%,n%,Pi,a,v
width 0,800
height 0,800

rem coordonnées du centre
x0=400
Y0=400
color 0,0,0,0
rem rayon
r=380

rem nombre de côtés
n%=7
for v=1 to 255
r=r-3
rem angle de début
a=a+50/57.29

Pi = 3.1415926535898

Px = r * Cos(a) + x0
Py = r * Sin(a) + y0

2D_poly_from Px,Py

For j% = 1 To n% :wait 1
    2d_pen_color j% * v/8,rnd(125), v
    Px2 = r*Cos(a + 2*Pi*j%/n%) + x0
    Py2 = r*Sin(a + 2*Pi*j%/n%) + y0
    2D_poly_to Px2, Py2
    Px = Px2
    Py = Py2
Next j%
next v
end

dim beta,gama,f,fx,fy ,R
dim a,ax,ay,b,bx,by
dim c,cx,cy,e,ex,ey
dim x0,y0 ,x
height 0,1000
width  0,1000
R=300
x0=500
y0=500
 ' Angle de depart
beta=8/57
 color 0,0,0,0
'  demi petit cote gama ne depass pas  90

2d_circle x0,y0,r
            for x=1 to 500
  beta=beta+186/57
  gama=gama+ 180/57
  if key_down_code =27 then end
 fx=x0+(sin(beta))
 fy=y0+(cos(beta)*r)
 ' 2d_circle fx,fy,5
 bx=x0-sin(beta+gama)*r
 by=y0+cos(beta+gama)
' 2d_circle bx,by,8
 wait 50
display
 cx= x0-(sin(beta+gama))
 cy= y0-(cos(beta+gama))*r
' 2d_circle cx,cy,10
 2d_line bx,by ,cx,cy
 2d_pen_color rnd( 255),rnd(125),rnd(240)
 
            next x

dim beta,gama,bx,by,cx,cy,R,x0,y0,x

height 0,800 : width  0,800 :  color 0,0,0,0

R=280 : x0=400 : y0=400 : beta=8/57

for x=1 to 1020
  2d_pen_color 255 - x/4,255 - x/4,abs(255 - x/2)
  beta=beta+172/57 :  gama=gama+ 182/57
  if key_down_code =27 then end
  bx=x0-sin(beta+gama)*r : by=y0+cos(beta+gama)
  cx= x0-sin(beta+gama) : cy= y0-cos(beta+gama)*r
  2d_line bx,by,cx,cy
  display : wait 10
next x

' ******************************************************************************
' *
' *          Les N insectes amoureux  (SUITE & FIN)
' *                par Papydall
' *
' ******************************************************************************
dim titre$ : titre$ = "Les N insectes amoureux (suite et fin)"
width 0,700 : height 0,700 : caption 0,"Papydall vous présente        " + titre$
dim N, DC,NbRep ,DegPol
' ******************************************************************************
' Ces constantes peuvent être modifiées
DC = 16      : ' diviseur du coté
NbRep = 20  : ' Nombre de repetitions
DegPol = 12  : ' Degre du polygone

' ******************************************************************************
N = 3        : ' Nombre d'insectes sur les sommets du triangle
dim X(N+1),Y(N+1),ar,pi,angle,i,j,k,x1,x2,y1,y2
dim lx,ly,xmin,xmax,ymin,ymax
pi = 4*atn(1) : ar = 2*pi/DegPol
lx = 500 : ly = 500                          : ' Dimension de l'écran graphique
xmin = -1 : xmax = 1 : ymin = -1 : ymax = 1 : ' Coordonnées des bords de l'écran
' ******************************************************************************


2d_fill_color 0,0,0 : 2d_rectangle 50,50,650,650

for k = 1 to DegPol
    x(1) = 0 : y(1) = 0 : angle = (k-1)* ar
    x(2) = cos(angle) : y(2) = sin(angle) : angle = k * ar
    x(3) = cos(angle) : y(3) = sin(angle)
    x(4) = x(1) : y(4) = y(1)
    for j = 1 to  NbRep
      x1 = lx*(x(1)-xmin)/(xmax-xmin)
      y1 = ly*(y(1)-ymin)/(ymax-ymin)
      2d_poly_from x1+100,y1+100
        for i = 2 to  N+1
            x2 = lx*(x(i)-xmin)/(xmax-xmin)
            y2 = ly*(y(i)-ymin)/(ymax-ymin)
            2d_pen_color 50*i,10*j,250
            2d_poly_to x2+100,y2+100
        next i
        if (k/2) <> int (k/2)        : ' Triangle impair ?
          for i = N+1 to 2 step -1  : ' Rotation sens trigonométrique
              x(i) = x(i)+(x(i-1)-x(i))/DC
              y(i) = y(i)+(y(i-1)-y(i))/DC
          next i
          x(1) = x(n+1) : y(1) = y(n+1)
        else
          for i = 1 to n : ' Triangle pair --> Rotation sens horraire
              x(i) = x(i) + (x(i+1) - x(i))/DC
              y(i) = y(i) + (y(i+1) - y(i))/DC
          next i
          x(n+1) = x(1) : y(n+1) = y(1)
        end_if
    next j
next k

end

' ******************************************************************************
' *
' *          Les N insectes amoureux  (SUITE & FIN)
' *                par Papydall
' *
' ******************************************************************************
dim titre$ : titre$ = "Les N insectes amoureux (suite et fin)"
width 0,700 : height 0,700 : caption 0,"Papydall vous présente        " + titre$
dim N, DC,NbRep ,DegPol
' ******************************************************************************
' Ces constantes peuvent être modifiées
DC = 16      : ' diviseur du coté
NbRep = 50  : ' Nombre de repetitions
DegPol = 20  : ' Degre du polygone

' ******************************************************************************
N = 3        : ' Nombre d'insectes sur les sommets du triangle
dim X(N+1),Y(N+1),ar,pi,angle,i,j,k,x1,x2,y1,y2
dim lx,ly,xmin,xmax,ymin,ymax
pi = 4*atn(1) : ar = 2*pi/DegPol
lx = 500 : ly = 500                          : ' Dimension de l'écran graphique
xmin = -1 : xmax = 1 : ymin = -1 : ymax = 1 : ' Coordonnées des bords de l'écran
' ******************************************************************************


2d_fill_color 0,0,0 : 2d_rectangle 50,50,650,650

for k = 1 to DegPol
    x(1) = 0 : y(1) = 0 : angle = (k-1)* ar
    x(2) = cos(angle) : y(2) = sin(angle) : angle = k * ar
    x(3) = cos(angle) : y(3) = sin(angle)
    x(4) = x(1) : y(4) = y(1)
    for j = 1 to  NbRep
      x1 = lx*(x(1)-xmin)/(xmax-xmin)
      y1 = ly*(y(1)-ymin)/(ymax-ymin)
      2d_poly_from x1+100,y1+100
        for i = 2 to  N+1
            x2 = lx*(x(i)-xmin)/(xmax-xmin)
            y2 = ly*(y(i)-ymin)/(ymax-ymin)
            if x2 < 350
              2d_pen_color abs(160 - 80*i),5*j,abs(200-k*20)
            else
              2d_pen_color 40*i,5*j,abs(200-k*20)
            end_if
            2d_poly_to x2+100,y2+100
        next i
        if (k/2) <> int (k/2)        : ' Triangle impair ?
          for i = N+1 to 2 step -1  : ' Rotation sens trigonométrique
              x(i) = x(i)+(x(i-1)-x(i))/DC
              y(i) = y(i)+(y(i-1)-y(i))/DC
          next i
          x(1) = x(n+1) : y(1) = y(n+1)
        else
          for i = 1 to n : ' Triangle pair --> Rotation sens horraire
              x(i) = x(i) + (x(i+1) - x(i))/DC
              y(i) = y(i) + (y(i+1) - y(i))/DC
          next i
          x(n+1) = x(1) : y(n+1) = y(1)
        end_if
    next j
next k

end

dim i,j,k,m,r,x1,x2,y1,y2,largx,hauty
width 0,600:height 0,600:color 0,0,0,0:display
largx=width(0)/2:hauty=height(0)/2:k=57.29577951308
2d_pen_width 2
for m=4 to 12 step 2
  for i=0.05 to 360 step 0.05
      j=sin(m*i/k)
      r=100*sqr(abs(j))
  if j<0 then r=-1*r
      2d_pen_color 70+i/2,0,70+i/2
      x1=largx+r*cos(i/k):y1=hauty+r*sin(i/k)
      x2=largx+cos(i):y2=hauty+sin(i)
      2d_line x1,y1,x2,y2
      if scancode=27 then terminate
  next i
  2d_fill_color 200,200,0
  2d_circle largx,hauty,15
  wait 800:cls
next m
terminate

dim i,j,k,m,r,x1,x2,y1,y2,largx,hauty
width 0,600:height 0,600:color 0,0,0,0:display
largx=width(0)/2:hauty=height(0)/2:k=57.29577951308
2d_pen_width 2
for m=4 to 12 step 2
  for i=0.05 to 360 step 0.05
      j=sin(m*i/k)
      r=100*sqr(abs(j))
      if j<0 then r=-1*r
      2d_pen_color abs(255-i*1.4),0, 0
      ' 2d_pen_color abs(180-i),0,abs (i-180)
      x1=largx+r*cos(i/k):y1=hauty+r*sin(i/k)
      x2=largx+cos(i):y2=hauty+sin(i)
      2d_line x1,y1,x2,y2
      if scancode=27 then terminate
  next i
  for i = 0 to 200 step 2
    2d_pen_color i,i,0 : 2d_fill_color i,i,0 : 2d_circle largx,hauty, 20 - (i / 10)
    wait 5
  next i
  if m < 12 then wait 1200:cls
next m
end

' ****************************************************
' *                                                                                       
' *          L'esthétique polaire
' *                par Papydall                                                   
' *                                                                                       
' ****************************************************
dim titre$ : titre$ = "L'esthétique polaire"
caption 0,"Papydall vous présente        " + titre$
width 0,700 : height 0,700
dim LX,LY,Xmin,Xmax,Ymin,Ymax,xp,yp,x,y,pi,t,f,periode,pas
pi = 4*atn(1)
LX = 500 : LY = 500  :                        ' dimentions de l'écran graphique
Xmin = -1 : XMax = 1 : Ymin = -1 : YMax = 1 : ' coordonnées des bords de l'écran
label  fonction
2d_pen_color 255,0,0 : 2d_fill_color 0,0,0 : 2d_rectangle 0,0,700,700
2d_poly_from lx+50,ly/2+50

periode = 20*pi : pas = pi/180
for t = 0 to periode +.1*pas step pas
    gosub fonction
    x = f*cos(t) : y = f*sin(t)
    xp = lx*(x-xmin)/(xmax-xmin) : yp = ly*(y-ymin)/(ymax-ymin)
    2d_poly_to  xp+50, yp+50
next t
end

fonction:
f = cos(9*t/10)
' f = cos(2*t/9) * cos(t)
' f = cos(t/10)
' f = 1/(sin(5*t)+.00001)
' f = sin(2*t)/(2+cos(t))
' f = sin(11*t/13)*cos(t)
' f = cos(sin(t/2))*sin(t/10)
return

' ****************************************************
' *
' *          L'esthétique polaires
' *                par Papydall
' *
' ****************************************************
dim titre$ : titre$ = "L'esthétique polaire"
caption 0,"Papydall vous présente        " + titre$
width 0,700 : height 0,700
dim LX,LY,Xmin,Xmax,Ymin,Ymax,xp,yp,x,y,pi,t,f,periode,pas
pi = 4*atn(1)
LX = 500 : LY = 500  :                        ' dimentions de l'écran graphique
Xmin = -1 : XMax = 1 : Ymin = -1 : YMax = 1 : ' coordonnées des bords de l'écran
label  fonction
2d_pen_color 0,0,0 : 2d_fill_color 0,0,0 : 2d_rectangle 0,0,700,700
2d_poly_from lx+50,ly/2+50

periode = 20*pi : pas = pi/180
for t = 0 to periode +.1*pas step pas
    gosub fonction
    x = f*cos(t) : y = f*sin(t)
    xp = lx*(x-xmin)/(xmax-xmin) : yp = ly*(y-ymin)/(ymax-ymin)
    2d_poly_to  xp+50, yp+50
    if yp + 250 > LY
      2d_pen_color abs(255*cos(t)/3),0,abs(200-(LY-YP))
    else
      2d_pen_color abs(350-(LY-YP)),0,abs(255*cos(t)/3)
    end_if

next t
end

fonction:
 f = cos(9*t/10)
' f = cos(2*t/9) * cos(t)
' f = cos(4*t/9) * cos(2*t)
' f = sin(2*t/9) * cos(4*t)
' f = cos(10*t) * cos(5*t)
' f = cos(2*t/13) * cos(13*t)
' f = cos(t/10)
' f = sin(2*t)/(2+cos(t))
' f = sin(11*t/13)*cos(t)
' f = cos(sin(t/2))*sin(t/10)
return

dim a,b,x,y
b=20
width 0,500:height 0,500:color 0,0,0,0:display
2d_pen_color 255,0,0
for a=0 to 200 step .01
  b=b+0.01:x=(b*sin(a/0.01))+225:y=(b*(tan(sin(a)/2)))+225
  2d_point x,y:2d_point y,x
next a

dim a,b,x,y : b=30
width 0,400:height 0,420:color 0,0,0,0
for a=0 to 400 step .025
  b=b+0.01:x=(b*sin(a/0.01))+190:y=(b*(tan(sin(a)/2)))+190
' 2d_pen_color abs(250 - x * 0.8), abs(255 - a), abs( 180 - y)
 2d_pen_color abs(b*sin(a/0.01)),255-a/2,a/2
 2d_point x,y:2d_point y,x:display
next a