/* --------------------------------------------------------------------------
    Function for testing root solver.
----------------------------------------------------------------------------- */
#include "tools.h"
#include "ffam.h"

void 
setup( void )      
{   
    printf( "\n No parameters needed for this function.\n" );
}

/* --------function to compute heat transfer discrepency------------------ */
double 
f( double r )      
{   
    double x = pow( fabs( r ), 0.333 );
    return (r < 0) ? -x : x;
}