Tracé d'un ellipsoïde

rem ============================================================================
rem                  Ellipsoïde
rem ============================================================================

Init()
Calcul()
caption 0,"Terminé"
end
rem ============================================================================
SUB Init()
    dim a,b,c,x,y,z,z1,z2,xc,yc,u1,v1, alfa,beta,zoom,u(3),v(3),w(3),xp(3)
    width 0,1000 : height 0,800 : color 0,0,255,255
    picture 10 : width 10,width(0)-100 : height 10,height(0)-100
    top 10,20 : left 10,50 : 2d_target_is 10 : 2d_clear : 2d_pen_color 255,0,0
    color 10,0,0,0
    degrees
    caption 0,"Veuillez patienter           <ESC> pour arrêter"
    xc = width(10)/2 : yc = height(10)/2
    Alfa = 22.5 : ' pi/8
    Beta = 45   : ' pi/4
    Zoom = 25   : z1 = zoom/width(10) : z2 = zoom/height(10)
    u(1) = COS(Alfa) * COS(Beta) : u(2) = SIN(Alfa) * COS(Beta) : u(3) = SIN(Beta)
    v(1) = 0-SIN(Alfa) : v(2) = COS(Alfa) : v(3) = 0
    w(1) = 0-COS(Alfa) * SIN(Beta) : w(2) = 0-SIN(Alfa)* SIN(Beta):w(3)= COS(Beta)
END_SUB
rem ============================================================================
SUB Calcul()
    for u1 = 0 to 360/2 step 2      : ' modifier le step selon la vitesse
        if u1 = 0 then 2d_pen_color 0,0,0 :else : 2d_pen_color 255,0,0
        for v1 = 1 to 360  step 2  : ' ou le finesse du tracé
            Ellipsoide(20,20,30)
          '  Ellipsoide(10,20,30)
          '  Ellipsoide(20,20,5)
          '  Ellipsoide(5,5,10)
          '  Ellipsoide(10,10,10)
          '  Fonction(20,15,5)
            Tracer()
            if scancode = 27 then terminate
       next v1
    next u1
END_SUB
rem ============================================================================
SUB Tracer()
    dim_local xe,ye
    xp(1) = u(1) * x + u(2) * y + u(3) * z
    xp(2) = v(1) * x + v(2) * y + v(3) * z
    xp(3) = w(1) * x + w(2) * y + w(3) * z
    xe = xc * xp(2) * z1 + yc : ye = yc *(1-xp(3) * z2)
    2d_poly_to xe,ye : display
END_SUB
rem ============================================================================
' Ellipsoïde
' Si a = b : ellipsoïde de révolution applati (forme de galet ou de soucoupe volonte)
' Si b = c : ellipsoïde de révolution allongé (forms de ballon de rugby ou de cigare)
' Si a = b = c : Sphère
SUB Ellipsoide(a,b,c)
    X = a * cos(v1) * cos(u1)
    Y = b * cos(v1) * sin(u1)
    Z = c * sin(v1)
END_SUB
rem ============================================================================
SUB Fonction(a,b,c)
    x = a*sin(u1)*cos(v1)
    y = b*sin(u1)*sin(v1)
    z = c*cos(u1)*log(u1+1)
END_SUB
rem ============================================================================

rem ============================================================================
rem                  Ellipsoïde
rem    code à utiliser pour une version inférieure à V 09.27i2
rem ============================================================================

Init()
Calcul()
caption 0,"Terminé"
end
rem ============================================================================
SUB Init()
    dim a,b,c,x,y,z,z1,z2,xc,yc,u1,v1, alfa,beta,zoom,u(3),v(3),w(3),xp(3)
    dim pi2,p
' &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
    dim pi   : ' si votre version de Panoramic_Editor est antérieure à V0.9.27i2
'                laissez cette ligne.
'                Sinon supprimez-la ou mettez-la en REM
' &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
    width 0,1000 : height 0,800 : color 0,0,255,255
    picture 10 : width 10,width(0)-100 : height 10,height(0)-100
    top 10,20 : left 10,50 : 2d_target_is 10
    color 10,0,0,0
    pi = acos(-1) : p = pi/90  : ' <-- = 2°
    caption 0,"Veuillez patienter           <ESC> pour arrêter"
    xc = width(10)/2 : yc = height(10)/2
    pi2 = 2*pi
    Alfa = pi/8
    Beta = pi/4
    Zoom = 25   : z1 = zoom/width(10) : z2 = zoom/height(10)
    u(1) = COS(Alfa) * COS(Beta) : u(2) = SIN(Alfa) * COS(Beta) : u(3) = SIN(Beta)
    v(1) = 0-SIN(Alfa) : v(2) = COS(Alfa) : v(3) = 0
    w(1) = 0-COS(Alfa) * SIN(Beta) : w(2) = 0-SIN(Alfa)* SIN(Beta):w(3)= COS(Beta)
END_SUB
rem ============================================================================
SUB Calcul()
    for u1 = 0 to pi2 step p      : ' modifier le step selon la vitesse
        if u1 = 0 then 2d_pen_color 0,0,0 :else : 2d_pen_color 255,0,0
        for v1 = 0 to pi2  step p  : ' ou le finesse du tracé
            Ellipsoide(20,20,30)
          '  Ellipsoide(10,20,30)
          '  Ellipsoide(20,20,5)
          '  Ellipsoide(5,5,10)
          '  Ellipsoide(10,10,10)
          '  Fonction(20,15,5)
            Tracer()
            if scancode = 27 then terminate
       next v1
    next u1
END_SUB
rem ============================================================================
SUB Tracer()
    dim_local xe,ye
    xp(1) = u(1) * x + u(2) * y + u(3) * z
    xp(2) = v(1) * x + v(2) * y + v(3) * z
    xp(3) = w(1) * x + w(2) * y + w(3) * z
    xe = xc * xp(2) * z1 + yc : ye = yc *(1-xp(3) * z2)
    2d_poly_to xe,ye : display
END_SUB
rem ============================================================================
' Ellipsoïde
' Si a = b : ellipsoïde de révolution applati (forme de galet ou de soucoupe volonte)
' Si b = c : ellipsoïde de révolution allongé (forms de ballon de rugby ou de cigare)
' Si a = b = c : Sphère
SUB Ellipsoide(a,b,c)
    X = a * cos(v1) * cos(u1)
    Y = b * cos(v1) * sin(u1)
    Z = c * sin(v1)
END_SUB
rem ============================================================================
SUB Fonction(a,b,c)
    x = a*sin(u1)*cos(v1)
    y = b*sin(u1)*sin(v1)
    z = c*cos(u1)*log(u1+1)
END_SUB
rem ============================================================================

dim pi : pi = acos(-1)
dim p  : p = pi/180