example.c

#include <stdio.h>
#include <stdlib.h>
#include <pds/pdsmath.h>
#include <float.h>
    
int main(int argc, char** argv)
{
        double x,y,z;
    
        for(x=0.001;x<0.5;x=x+0.05)
        {
                y=pds_hb(x);
                z=pds_hbinv(y);         
    
                printf("x=%f\ty=hb(x)=%f\t",x,y);
                printf("z=hbinv(y)=%f\t",z);
                printf("ERROR=%6.4f ppm\n",(x-z)*1000000/x);
        }
    
        return EXIT_SUCCESS;
}