ssi - Sine Integral Si(x)

#include <msca.h> 
 
double ssi (double x); 
float  ssif(float x);

Given a value $x$ this function computes the Sine Integral:

∫ xsin(x) Si(x)= 0 x dx