g++ выдает nan

#include <stdio.h>
#include <math.h>

float plus(float x, float z)
{
  return x+z;
}
float del(float x, float z)
{
  return x/z;
}
float umn(float x, float z)
{
  return x*z;
}
int ravno(float n, float k)
{
  return (n==k)?1:0;
}
int bolshe(float n, float k)
{
  return (n>k)?1:0;
}
int menshe(float n, float k)
{
  return (n<k)?1:0;
}
float kvadrat(float x)
{
  return umn(x,x);
}
float abs(float x)
{
  if(bolshe(x,0))
    return x;
  else
    return -x;
}

float a0(float x) {return x;}

float b0(float z) {return z;}

float anext(float a, float x, int k)
{
  return umn(a, umn(-1, del(umn(kvadrat(x), plus(umn(2, k), 1)), plus(umn(2,k), 3))));
}

float bnext(float b, float z, int k)
{
  return umn(b, umn(-1, del(umn(kvadrat(z), plus(umn(2, k), 1)), plus(umn(2,k), 3))));
}

float y_telo(float x, float z, float acount, float bcount, float a, float b, float a_next, float b_next, float eps, int k)
{
  if(menshe(abs(plus(a_next,b_next)), eps))
  {
    return plus(plus(acount, bcount), plus(a_next, b_next));
  }
  else
  {
    y_telo(x,z, acount+a_next, bcount+b_next, a_next,b_next,anext(a_next,x,plus(k,1)),bnext(b_next,z,plus(k,1)), eps, plus(k,1));
  }
}
float yfunc(float x, float z, float eps)
{
  return y_telo(x,z,a0(x),b0(z),a0(x),b0(z),anext(a0(x),x,0),bnext(b0(z),z,0), eps, 0);
}
int main()
{
  float x = 0.2; float z = 0.3; float eps = 0.1;
  printf("%f\n", yfunc(x,z,eps));
  printf("%f\n", atan((x+z)/(1-x*z)));
  return 0;
}

g++ -o lab1 lab1.c