Critére de Dang Van en fatigue multiaxiale

Bonjour,

bonjour,
Je vous écrit pour solliciter votre aide sur la programmation du critére de Dang Van en fatigue multiaxiale.
J'ai réussi a programmer le critére en lui meme, j'ai pas encore réussi a programmer le plus petit cercle circonscrit, en voila l'algorithem :

c etape1 : on cherche les 2 pts p1 et p2 du trajet de charge les
c plus eloign's l''un de l''autre
c soit C le cercle de diam'tre [p1,p2]

c etape2 : on teste si tous les points du trajet de charge sont
c ... l''int'rieur du C(I,p1p2/2)
c si c''est le cas, on retient C
c sinon on cherche le point p3 le plus 'loign' de C
c et on determine le cercle passant par (p1,p2,p3)

c etape3 : on teste si tous les pts du trajet de charge st ...
c l''int'rieur de C
c si oui c''est bon
c sinon on cherche le pt p4 le plus 'loign' de C
c parmi les cercles construits sur p1p2p3 , p1p3p4 , p2p3p4
c on retient le + petit cercle contenant les 4 pts
c p1,p2,p2 et p4 => c le nouveau cercle
c le point disparu (ex : p3 si c''est le cercle bas' sur p1,p2,p4)
c est remplac' par p4
c on reboucle (retourne) sur l''etape3

Je suis débutant sur FORTRAN 77.
Merci d'avance