Equation de la pomme

rem ============================================================================
rem         Equation de la Pomme
rem ============================================================================
picture 10 : full_space 10 : 2d_target_is 10
Pomme(300,120,1,1.2,0.15,0.08)
end
rem ============================================================================
SUB Pomme(xc,yc,a,b,p,q)
    dim_local x,y,t,dt,r,s,pi2
    dt = pi/180 : s = 200 : pi2 = pi/2
    2d_pen_color 255,0,0
    for t = 0-pi/2 to 3*pi/2 step dt
        r = a*(1-sin(t))
        x = r*exp(0-p*(t-pi2)*(t-pi2))*cos(t)
        y = b*r*exp(0-q*(t-pi2)*(t-pi2))*sin(t)
        if t = 0-pi/2
           2d_poly_from xc+s*x,yc-s*y
        else
           2d_poly_to xc+s*x,yc-s*y
        end_if
    next t
    2d_flood xc,2*yc,255,0,0
END_SUB
rem ============================================================================