home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / SPECDEMO.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-01-21  |  6.2 KB  |  293 lines
  1. # include <stdio.h>
  2. # include "miscio.h"
  3. # include <string.h>
  4. # include <math.h>
  5. # include "specfunc.h"
  6.  
  7. int i;
  8. int option;
  9. int InputNum;
  10. int errflag;
  11. float x;
  12. float y;
  13. float z;
  14. float w;
  15. float CalcAns;
  16. char done;
  17.  
  18.  
  19.  
  20. void DisplayOptions()
  21. {
  22.  
  23.    ClrScr();
  24.    GotoXY(20,1);
  25.    printf( "Special Functions Demo Programs");
  26.    GotoXY(5,3);
  27.    printf( "1.  Log Gamma");
  28.    GotoXY(5,4);
  29.    printf( "2.  Gamma");
  30.    GotoXY(5,5);
  31.    printf( "3.  Incomplete Gamma");
  32.    GotoXY(5,6);
  33.    printf( "4.  Incomplete Gamma Complement");
  34.    GotoXY(5,7);
  35.    printf( "5.  Bessel");
  36.    GotoXY(5,8);
  37.    printf( "6.  Tangent");
  38.    GotoXY(5,9);
  39.    printf( "7.  Hyperbolic Cosine");
  40.    GotoXY(5,10);
  41.    printf( "8.  Hyperbolic Sine");
  42.    GotoXY(5,11);
  43.    printf( "9.  Hyperbolic Secant");
  44.    GotoXY(5,12);
  45.    printf( "10. Hyperbolic ArcTangent");
  46.    GotoXY(5,13);
  47.    printf( "11. Error Function (Iterative)");
  48.    GotoXY(5,14);
  49.    printf( "12. Error Function");
  50.    GotoXY(37,3);
  51.    printf( "13. Complement Error Function");
  52.    GotoXY(37,4);
  53.    printf( "14. Error Function (Real)");
  54.    GotoXY(37,5);
  55.    printf( "15. Error Function (Imaginary)");
  56.    GotoXY(37,6);
  57.    printf("16. Hermite Polynomial ");
  58.    GotoXY(37,7);
  59.    printf( "17. Legendre Polynomial ");
  60.    GotoXY(37,8);
  61.    printf( "18. LaGuerre Polynomial");
  62.    GotoXY(37,9);
  63.    printf( "19. Modified Bessel of I ");
  64.    GotoXY(37,10);
  65.    printf( "20. Modified Bessel of K ");
  66.    GotoXY(37,11);
  67.    printf( "21. Jacobi Polynomial ");
  68.    GotoXY(37,12);
  69.    printf( "22. Tchebyshev Polynomial ");
  70.    GotoXY(37,13);
  71.    printf( "23. Beta ");
  72.    GotoXY(37,14);
  73.    printf( "24. Incomplete Beta ");
  74.    GotoXY(37,15);
  75.    printf( "25. END");
  76.    GotoXY(1,25);
  77.    printf( "Refer to Chapter 17 of Manual for details concerning function arguments");
  78. }
  79.  
  80. void AskForOption(int *option)
  81. {
  82.    char done;
  83.  
  84.    do {
  85.       GotoXY(1,17);
  86.       ClrEol();
  87.       GotoXY(5,17);
  88.       ClrEol();
  89.       printf("ENTER SELECTION:  ");
  90.       scanf("%d", option);
  91.    } while ( ! (*option >= 1) && (*option <= 25) );
  92. }
  93.  
  94. void NumOfInputs(int option,int *InputNum)
  95. {
  96.  
  97.    switch ( option ) {
  98.       case 1:
  99.       case 2:
  100.       case 6:
  101.       case 7:
  102.       case 8:
  103.       case 9:
  104.       case 10:
  105.       case 11:
  106.       case 12:
  107.       case 13:
  108.          (*InputNum) = 1;
  109.          break;
  110.       case 3:
  111.       case 4:
  112.       case 5:
  113.       case 14:
  114.       case 15:
  115.       case 16:
  116.       case 17:
  117.       case 19:
  118.       case 20:
  119.       case 22:
  120.       case 23:
  121.          (*InputNum) = 2;
  122.          break;
  123.       case 18:
  124.       case 24:
  125.          (*InputNum) = 3;
  126.          break;
  127.       case 21:
  128.          (*InputNum) = 4;
  129.          break;
  130.    }
  131. }
  132.  
  133. void ReadReals(int InputNum,float *x,float *y,float *z,float *w)
  134. {
  135.  
  136.    GotoXY(1,24);
  137.    printf( "Enter arguments in the order they appear in manual. ");
  138.    GotoXY(5,18);
  139.    printf( "ENTER VALUES FOR :");
  140.    GotoXY(1,20);
  141.    ClrEol();
  142.    printf( "Parameter 1: ");
  143.    scanf("%f", x);
  144.    if ( InputNum > 1 ) {
  145.       GotoXY(20,20);
  146.       printf( "Parameter 2: ");
  147.       scanf( "%f",y);
  148.       if ( InputNum > 2 ) {
  149.          GotoXY(40,20);
  150.          printf( "Parameter 3: ");
  151.          scanf("%f", z);
  152.          if ( InputNum == 4 ) {
  153.             GotoXY(60,20);
  154.             printf( "Parameter 4: ");
  155.             scanf("%f", w);
  156.          }
  157.       }
  158.    }
  159. }
  160.  
  161. void CalculateAnswer(int option,float x,float y,float z,float w,int *errflag)
  162. {
  163.    int Int1 = 21;
  164.    int Int2;
  165.  
  166.    CalcAns = 0;
  167.    (*errflag) = 0;
  168.    switch ( option ) {
  169.       case 1:
  170.          CalcAns = LogGamma(x);
  171.          break;
  172.       case 2:
  173.          CalcAns = Gamma(x,errflag);
  174.          break;
  175.       case 3:
  176.          CalcAns = IncGamma(x,y);
  177.          break;
  178.       case 4:
  179.          CalcAns = IncGammaComp(x,y);
  180.          break;
  181.       case 5:
  182.          CalcAns = Bessel(x,y,errflag);
  183.          break;
  184.       case 6:
  185.          CalcAns = Tan(x);
  186.          break;
  187.       case 7:
  188.          CalcAns = Cosh(x);
  189.          break;
  190.       case 8:
  191.          CalcAns = Sinh(x);
  192.          break;
  193.       case 9:
  194.          CalcAns = Sech(x);
  195.          break;
  196.       case 10:
  197.          CalcAns = ArcTanh(x);
  198.          break;
  199.       case 11:
  200.          CalcAns = ErrFuncIter(x);
  201.          break;
  202.       case 12:
  203.          CalcAns = ErrFunc(x,errflag);
  204.          break;
  205.       case 13:
  206.          CalcAns = ErrFuncComp(x,errflag);
  207.          break;
  208.       case 14:
  209.          CalcAns = ErrFuncR(x,y,errflag);
  210.          break;
  211.       case 15:
  212.          CalcAns = ErrFuncI(x,y);
  213.          break;
  214.       case 16:
  215.          {
  216.             Int1 = floor(x);
  217.             CalcAns = Hermite(Int1,y,errflag);
  218.          }
  219.          break;
  220.       case 17:
  221.          {
  222.             Int1 = floor(x);
  223.             CalcAns = Legend(Int1,y,errflag);
  224.          }
  225.          break;
  226.       case 18:
  227.          {
  228.             Int1 = floor( x);
  229.             CalcAns = Laguerre(Int1,y,z,errflag);
  230.          }
  231.          break;
  232.       case 19:
  233.          CalcAns = ModBesselI(x,y,errflag);
  234.          break;
  235.       case 20:
  236.          CalcAns = ModBesselK(x,y,errflag);
  237.          break;
  238.       case 21:
  239.          {
  240.             Int1 = floor( x);
  241.             CalcAns = Jacobi(Int1,y,z,w,errflag);
  242.          }
  243.          break;
  244.       case 22:
  245.          {
  246.             Int1 = floor(x);
  247.             CalcAns = Tcheb(Int1,y,errflag);
  248.          }
  249.          break;
  250.       case 23:
  251.          CalcAns = Beta(x,y,errflag);
  252.          break;
  253.       case 24:
  254.          CalcAns = IncBeta(x,y,z);
  255.          break;
  256.    }
  257.    GotoXY(5,22);
  258.    printf("ANSWER: %f", CalcAns);
  259.    printf("\n");
  260.    getch();
  261. }
  262.  
  263.  
  264. void main()
  265. {
  266.    option = -1;
  267.    InputNum = -1;
  268.    x = -1.0;
  269.    y = -1.0;
  270.    z = -1.0;
  271.    w = -1.0;
  272.    DisplayOptions();
  273.    do {
  274.       AskForOption(&option);
  275.       if ( option != 25 ) {
  276.          NumOfInputs(option,&InputNum);
  277.          ReadReals(InputNum,&x,&y,&z,&w);
  278.          CalculateAnswer(option,x,y,z,w,&errflag);
  279.          if ( errflag != 0 ) {
  280.             DisplayErrorMessage(errflag);
  281.          }
  282.          GotoXY(1,24);
  283.          ClrEol();
  284.          printf( "Press Return to Continue  ");
  285.          scanf("");
  286.          for ( i = 16; i <= 24; ++i ) {
  287.             GotoXY(1,i);
  288.             ClrEol();
  289.          }
  290.       }
  291.    } while ( ! (option == 25) );
  292. }
  293.