home *** CD-ROM | disk | FTP | other *** search
/ World of Shareware - Software Farm 2 / wosw_2.zip / wosw_2 / CPROG / ANUMR5.ZIP / EDIC1TST.C < prev    next >
C/C++ Source or Header  |  1991-08-12  |  1KB  |  74 lines
  1. /* OK mémoire */
  2.  
  3. #include <math.h>
  4. #include <conio.h>
  5. #include <stdio.h>
  6. #include <stdlib.h>
  7. #include <alloc.h>
  8.  
  9. #include "anum.h"
  10. #include "sysio.h"
  11.  
  12.  
  13.  
  14. double bf(double t, double x)
  15.  
  16. {    double res;
  17.  
  18.     res = x/t + t - 1;
  19.     return(res);
  20. }
  21.  
  22.  
  23.  
  24. void results(int nret, double *xval, double *tval)
  25.  
  26. {    int i;
  27.  
  28.     puts("       T          X");
  29.     puts("---------------------------");
  30.     for (i=0; i<nret; i++)
  31.         printf("% 10.5le  % 10.5le\n", *(tval+i), *(xval+i));
  32. }
  33.  
  34.  
  35.  
  36. void main(void)
  37.  
  38. {    double lb,ub,x0,*xval,*tval;
  39.     int nint,nret,errcode;
  40.  
  41.     clrscr();
  42.     puts("Test program for function initcond1storder");
  43.     puts("------------------------------------------\n");
  44.  
  45.     lb=2.0;
  46.     ub=3.0;
  47.     x0=4.0;
  48.     nret=5;
  49.     nint=180 ;
  50.  
  51.     printf("Lower bound            : %+10.5lf\n",lb);
  52.     printf("Upper bound            : %+10.5lf\n",ub);
  53.     printf("x at t = %+lf     : %+10.5lf\n\n",lb,x0);
  54.     printf("Number of asked values : %3d\n",nret);
  55.     printf("Number of intervals    : %3d\n\n",nint);
  56.  
  57.     xval=farcalloc(nret,sizeof(double));
  58.     tval=farcalloc(nret,sizeof(double));
  59.     if (xval==NULL || tval==NULL)
  60.     {    fputs("Not enough core", stderr);
  61.         exit(-1);
  62.     }
  63.  
  64.     initcond1storder(lb, ub, x0, nret, nint, tval, xval, &errcode, bf);
  65.     SYSMSG(errcode, stderr);
  66.  
  67.     results(nret, xval, tval);
  68.  
  69.     farfree(tval);
  70.     farfree(xval);
  71.  
  72. }
  73.  
  74.