+
+/* Supply a private version of the power function that returns 1
+ for x**0. This avoids c-1471. Abigail's Japh tests depend
+ on this fix. We leave all the other cases to the VOS C
+ runtime. */
+
+double s_crt_pow(double *x, double *y);
+
+double pow(x,y)
+double x, y;
+{
+ if (y == 0e0) /* c-1471 */
+ {
+ errno = EDOM;
+ return (1e0);
+ }
+
+ return(s_crt_pow(&x,&y));
+}