1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
#include <u.h>
#include <libc.h>
#include "hoc.h"
double errcheck(double, char*);
double
Log(double x)
{
return errcheck(log(x), "log");
}
double
Log10(double x)
{
return errcheck(log10(x), "log10");
}
double
Sqrt(double x)
{
return errcheck(sqrt(x), "sqrt");
}
double
Exp(double x)
{
return errcheck(exp(x), "exp");
}
double
Asin(double x)
{
return errcheck(asin(x), "asin");
}
double
Acos(double x)
{
return errcheck(acos(x), "acos");
}
double
Sinh(double x)
{
return errcheck(sinh(x), "sinh");
}
double
Cosh(double x)
{
return errcheck(cosh(x), "cosh");
}
double
Pow(double x, double y)
{
return errcheck(pow(x,y), "exponentiation");
}
double
integer(double x)
{
if(x<-2147483648.0 || x>2147483647.0)
execerror("argument out of domain", 0);
return (double)(long)x;
}
double
errcheck(double d, char* s) /* check result of library call */
{
if(isNaN(d))
execerror(s, "argument out of domain");
if(isInf(d, 0))
execerror(s, "result out of range");
return d;
}
|