Télécharger gamma.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gamma.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5.  
  6. ************************************************************************
  7. * NOM : GAMMA
  8. * DESCRIPTION : Teste l'operateur 'FONC' 'GAMMA' sur des LISTREEL
  9. *
  10. *
  11. * LANGAGE : GIBIANE-CAST3M
  12. * AUTEUR : Serge PASCAL (CEA/DEN/DM2S/SEMT/LM2S)
  13. * REVISION : Clément BERTHINIER (CEA/DEN/DM2S/SEMT/LM2S)
  14. **********************************************************************
  15. * VERSION : v1, 28/10/2009, version initiale
  16. * HISTORIQUE : v1, ??/??/????, création
  17. * HISTORIQUE : v2, 11/02/2015, Des valeurs Inf apparaissaient...
  18. * Correction de la plage de validite
  19. * HISTORIQUE :
  20. ************************************************************************
  21. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  22. * en cas de modification de ce sous-programme afin de faciliter
  23. * la maintenance !
  24. ************************************************************************
  25.  
  26.  
  27.  
  28. * Test de quelques valeurs remarquables :
  29. L1 = 'PROG' -1.5 -0.5 0.5 1. 1.5 2. 2.5 3. 3.5 4. ;
  30. RCPI1 = PI ** 0.5 ;
  31. L2 = 'PROG' (4. * RCPI1 / 3.) (-2. * RCPI1) RCPI1 1.
  32. (0.5 * RCPI1) 1. (0.75 * RCPI1) 2.
  33. (0.125 * 15. * RCPI1) 6. ;
  34. L3 = 'FONC' 'GAMMA' L1 ;
  35. ERR1 = 'MAXI' 'ABS' ((L3 - L2) / L2) ;
  36.  
  37. * Test egalite : G(x)G(1-x)=PI/sin(PI*x) :
  38. E1 = EXP 1. ;
  39. L1 = 'PROG' (-1.D0 * (E1** 5.05)) PAS (E1 / 140.) (E1 ** 5.1) ;
  40. L2 = SIN (180. * L1) ;
  41. L3 = 'FONC' 'GAMMA' (1 - L1) ;
  42. L4 = 'FONC' 'GAMMA' (L1) ;
  43. L5 = L2 * L3 * L4 - ('PROG' ('DIME' L1) * PI);
  44.  
  45. ERR2 = 'MAXI' 'ABS' L5 ;
  46. 'LIST' (ABS L5) ;
  47. 'LIST' ERR2 ;
  48. *'DESS' ('EVOL' 'MANU' L1 L23) ;
  49.  
  50. ERR0 = 'MAXI' ('PROG' ERR1 ERR2) ;
  51. 'SI' (ERR0 '>' 2.E-10) ;
  52. 'ERRE' 5 ;
  53. 'FINS' ;
  54.  
  55. 'FIN' ;
  56.  
  57.  
  58.  
  59.  
  60.  

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