REM 'calcul du rayon du cercle dont on connait 3 points A, B et C


CLS
CLEAR


'INPUT "      Valeur de xa (puis appuyez sur ENTREE): ", xa
'INPUT "      Valeur de ya (puis appuyez sur ENTREE): ", ya
'INPUT "      Valeur de xb (puis appuyez sur ENTREE): ", xb
'INPUT "      Valeur de yb (puis appuyez sur ENTREE): ", yb
'INPUT "      Valeur de xc (puis appuyez sur ENTREE): ", xc
'INPUT "      Valeur de yc (puis appuyez sur ENTREE): ", yc

xa = 0
ya = 970
xb = 882
yb = 1196
xc = 1729
yc = 1709

'calcul des coeff a et b des 2  mediatrices de AB et BC


a1 = (xb - xa) / (ya - yb)

a2 = (xc - xb) / (yb - yc)


b1 = (ya ^ 2 - yb ^ 2 + xa ^ 2 - xb ^ 2) / (2 * ya - 2 * yb)

b2 = (yb ^ 2 - yc ^ 2 + xb ^ 2 - xc ^ 2) / (2 * yb - 2 * yc)


'calcul des coord du point d intersection des 2 mediatrices


xi = (b2 - b1) / (a1 - a2)

yi = a1 * xi + b1

'calcul du rayon du cercle par 3 formules


r1 = SQR((xi - xa) ^ 2 + (yi - ya) ^ 2)
r12 = INT(r1)

PRINT "le rayon en pixels du cercle image du soleil est r1 = "; r12


r2 = SQR((xi - xb) ^ 2 + (yi - yb) ^ 2)

'PRINT r2

r3 = SQR((xi - xc) ^ 2 + (yi - yc) ^ 2)

'PRINT r3

'coord du centre I du cercle

'PRINT "xi = "; xi

'PRINT "yi = ";  yi

''''''''''''''''''''''''''''''''''''''''''''''

'calcul de la distance entre les 2 points brillants P et Q

xp = 1051
yp = 325
xq = 1132
yq = 56
'espace de temps entre les 2 images en minutes
t = 12


dist = SQR((xp - xq) ^ 2 + (yp - yq) ^ 2)

PRINT " distance en pixels entre les points brillants P et Q : ", dist

'coeff de conversion pixels en km ;

k = 700000 / r1
dist1 = k * dist
dist2 = INT(dist1)

PRINT " distance en km entre les points brillants P et Q : ", dist2


'calcul de la vitesse d'ejection


v1 = k * dist / (t * 60)
v2 = INT(v1)
PRINT " vitesse des gaz en km/s entre P et Q : "; v2



'calcul de la vitesse de liberation au dessus du soleil au point P :


'distance en pixels du point P au centre du soleil :

d = SQR((xi - xp) ^ 2 + (yi - yp) ^ 2)
d1 = INT(k * d)

PRINT " distance en km du point P au centre du soleil = "; d1
d2 = INT(d1 - 700000)
PRINT " altitude  en km du point P au-dessus de la surface solaire = "; d2

vd = 617.5 * SQR(r1 / d)
vd1 = INT(vd)

PRINT " vitesse de liberation au point P en km/s : "; vd1




END