Solve Exercise 4-5: adding sin(), exp(), pow()

1 files changed, 216 insertions, 0 deletions

diff --git a/ch4/4-05_math-funcs.c b/ch4/4-05_math-funcs.c
new file mode 100644
index 0000000..6c22767
--- /dev/null
+++ b/ch4/4-05_math-funcs.c
@@ -0,0 +1,216 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <ctype.h>
+#include <math.h>
+
+/* The C Programming Language: 2nd Edition
+ *
+ * Exercise 4-5: Add access to library functions like `sin`, `exp`, and
+ * `pow`. See <math.h> in Appendix B, Section 4.
+ *
+ * Answer: Doing this is much like 4-4; we're gonna use more special
+ * characters, and tie them to math.h functions.
+ *
+ * Note that math.h may be a separate library on your system. If you use
+ * GNU/Linux, the library may be called `libm` and can be linked to using
+ * the `-lm` gcc option when you compile. I'll remove this note when I
+ * learn how to write makefiles or something similar.
+ */
+
+#define MAXOP    100
+#define NUMBER   '0'
+#define MAXVAL   100
+#define BUFSIZE  100
+
+int getop(char []);
+void push(double);
+double pop(void);
+int getch(void);
+void ungetch(int);
+void stack_top(void);
+double dupe_top(void);
+void swap_top_two(void);
+void clear_stack(void);
+
+int sp = 0; // Next free stack position
+double val[MAXVAL]; // Value stack
+char buf[BUFSIZE]; // buffer for ungetch
+int bufp = 0; // next free position in buf
+
+/* Reverse Polish calculator:
+ *
+ * Binary operations (+-*\%)
+ *     operand operand operator
+ *
+ * Example: 6 minus 2 in Reverse Polish Notation is "6 2 -"
+ */
+int main() {
+	int type;
+	double op2;
+	char s[MAXOP];
+
+	while ((type = getop(s)) != EOF) {
+		switch (type) {
+			case NUMBER:
+				push(atof(s));
+				break;
+			case '+':
+				push(pop() + pop());
+				break;
+			case '*':
+				push(pop() * pop());
+				break;
+			case '-':
+				op2 = pop();
+				push(pop() - op2);
+				break;
+			case '/':
+				op2 = pop();
+				if (op2 != 0.0) {
+					push(pop() / op2);
+				} else {
+					printf("Error: Cannot divide by zero.\n");
+				}
+				break;
+			/* Yay for modulus! */
+			case '%':
+				op2 = pop();
+				if (op2 != 0.0) {
+					push((int)pop() % (int)op2);
+				} else {
+					printf("Error: Cannot modulo by zero.\n");
+				}
+				break;
+			/* Top of stack */
+			case '?':
+				stack_top();
+				break;
+			/* Dupe the top of the stack */
+			case '#':
+				dupe_top();
+				break;
+			/* Swap the top two */
+			case '~':
+				swap_top_two();
+				break;
+			/* Clear the stack */
+			case '!':
+				clear_stack();
+				break;
+			/* sin() support */
+			case '(':
+				op2 = sin(pop());
+				push(op2);
+				break;
+			/* exp() support */
+			case '{':
+				op2 = exp(pop());
+				push(op2);
+				break;
+			/* pow() support */
+			case '^':
+				op2 = pop();
+				push(pow(pop(), op2));
+				break;
+			case '\n':
+				printf("\t%.8g\n", pop());
+				break;
+			default:
+				printf("Error: Unknown command %s\n", s);
+				break;
+		}
+	}
+	return 0;
+}
+
+void push(double f) {
+	if (sp < MAXVAL) {
+		val[sp++] = f;
+	} else {
+		printf("Error: Stack full. Cannot push %g\n", f);
+	}
+}
+
+double pop(void) {
+	if (sp > 0) {
+		return val[--sp];
+	} else {
+		printf("Error: Stack empty.\n");
+		return 0.0;
+	}
+}
+
+int getop(char s[]) {
+	int i = 0;
+	int c, next;
+
+	while ((s[0] = c = getch()) == ' ' || c == '\t') {
+	}
+	s[1] = '\0';
+	/* The final check is for negative numbers. */
+	if (!isdigit(c) && c != '.' && c != '-') {
+		return c;
+	}
+	/* The second half of this if-statement accounts for negatives */
+	if (c == '-') {
+		next = getch();
+		if (!isdigit(next) && next != '.') {
+			return c;
+		} else {
+			c = next;
+		}
+	} else {
+		c = getch();
+	}
+
+	while (isdigit(s[++i] = c)) {
+		c = getch();
+	}
+	if (c == '.') {
+		while (isdigit(s[++i] = c = getch())) {
+		}
+	}
+	s[i] = '\0';
+	if (c != EOF) {
+		ungetch(c);
+	}
+	return NUMBER;
+}
+
+int getch(void) {
+	return (bufp > 0) ? buf[--bufp] : getchar();
+}
+
+void ungetch(int c) {
+	if (bufp >= BUFSIZE) {
+		printf("ungetch: Too many characters.\n");
+	} else {
+		buf[bufp++] = c;
+	}
+}
+
+void stack_top(void) {
+	if (sp > 0) {
+		printf("Top of stack is %8g\n", val[sp - 1]);
+	} else {
+		printf("Error: Stack empty.\n");
+	}
+}
+
+double dupe_top(void) {
+	double temp = pop();
+	push(temp);
+	push(temp);
+}
+
+void swap_top_two(void) {
+	double tmp1, tmp2;
+	tmp1 = pop();
+	tmp2 = pop();
+	push(tmp1);
+	push(tmp2);
+}
+
+void clear_stack(void) {
+	sp = 0;
+}