Télécharger la plaquette Cast3M


fr en

Télécharger snv.eso

Retour à la liste

Numérotation des lignes : 

snv

  1. C SNV       SOURCE    CB215821  16/04/21    21:18:29     8920
  2.       FUNCTION SNV(AJV,ITYPE,Gamma,Delta,Xlam,Xi,IFAULT)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5. C
  6. C
  7. C   Fonction de convertion v.a. Johnson -> v.a Normale
  8. C
  9. C
  10.       REAL*8 AJV, Gamma, Delta, Xlam, Xi, V, W, cons
  11. C
  12.       DATA cons /-63.D0/
  13. C
  14. C
  15.       SNV = 0.D0
  16.       IFAULT = 1
  17.       IF (ITYPE .LT. 1 .OR. ITYPE .GT. 4) RETURN
  18.       IFAULT = 0
  19.       GOTO (10, 20, 30, 40), ITYPE
  20. C
  21. C   Distribution Sl
  22. C
  23. C  10  W = Xlam * (AJV - Xi)
  24.   10  W = Xlam * AJV - Xi
  25.       IF (W .LE. 0.D0) GOTO 15
  26. C      SNV = Xlam * (LOG(W) * Delta + Gamma)
  27.       SNV = LOG(W) * Delta + Gamma
  28.       RETURN
  29.   15  IFAULT = 2
  30.       RETURN
  31. C
  32. C   Distribution Su
  33. C
  34.   20  W = (AJV - Xi) / Xlam
  35.       IF (W .GT. cons) GOTO 23
  36.       W = -0.5D0 / W
  37.       GOTO 27
  38.   23  W = SQRT(W * W + 1.D0) + W
  39.   27  SNV = LOG(W) * Delta + Gamma
  40.       RETURN
  41. C
  42. C   Distribution Sb
  43. C
  44.   30  W = AJV - Xi
  45.       V = Xlam - W
  46.       IF (W .LE. 0.D0 .OR. V .LE. 0.D0) GOTO 35
  47.       SNV = LOG(W / V) * Delta + Gamma
  48.       RETURN
  49.   35  IF (W .LE. 0.D0) SNV=-10.
  50.       IF (V .LE. 0.D0) SNV=10.
  51.       RETURN
  52. C
  53. C   Distribution Normale
  54. C
  55.   40  SNV = Delta * AJV + Gamma
  56.       RETURN
  57.       END
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales