serf - Error function erf(x)

#include <msca.h> 
 
double serf (double x); 
float  serff(float x);

Given a value $x$ this function computes the error function:

∫ ⎷2- ze-t2 dt π 0