00001
00020
00021
00022 # include <stdio.h>
00023 # include <stdlib.h>
00024 # include <math.h>
00025
00026 # include "global.h"
00027 # include "sub.h"
00028 # include "rand.h"
00029
00030
00031 long double maximum (long double a, long double b)
00032 {
00033 if (a>b)
00034 {
00035 return(a);
00036 }
00037 return (b);
00038 }
00039
00040
00041 long double minimum (long double a, long double b)
00042 {
00043 if (a<b)
00044 {
00045 return (a);
00046 }
00047 return (b);
00048 }
00049
00050
00051 long double modulus (long double *x, int n)
00052 {
00053 int i;
00054 long double res;
00055 res = 0.0;
00056 for (i=0; i<n; i++)
00057 {
00058 res += x[i]*x[i];
00059 }
00060 return (sqrt(res));
00061 }
00062
00063
00064 long double dot (long double *a, long double *b, int n)
00065 {
00066 int i;
00067 long double res;
00068 res = 0.0;
00069 for (i=0; i<n; i++)
00070 {
00071 res += a[i]*b[i];
00072 }
00073 return (res);
00074 }
00075
00076
00077 long double mean (long double *x, int n)
00078 {
00079 int i;
00080 long double res;
00081 res = 0.0;
00082 for (i=0; i<n; i++)
00083 {
00084 res += x[i];
00085 }
00086 return (res/(long double)n);
00087 }