/* --------------------------------- */
/* save as fc.h */
/* --------------------------------- */
#define LOOP 2*300
/* --------------------------------- */
double s = 1./2.;
double aa = 2.;
/* --------------------------------- */
double F(
double t)
{
return( pow(t,2) );
}
char Feq[] = "t**2";
/* --------------------------------- */
double Fa(
double t)
{
return( pow(aa*t,2) );
}
char Faeq[] = "(a*t)**2";
/* --------------------------------- */
/* ---------------------------------
Laplace transform of F(t)
--------------------------------- */
double f(
double s)
{
return( 2./pow(s,3));
}
char feq[] = "(2/s^3)";
/* ---------------------------------- */
double f_s(
double s)
{
return( (1./aa) * 2./pow(s/aa,3));
}
char f_seq[] = "(1/a) * (2/(s/a)^3)";
char f2seq[] = "2*a^2/(s^3)";
/* ---------------------------------- */
/* ---------------------------------- */
double b = 100.; char beq[] = "100";
double a = 0.; char aeq[] = "0";
/* ---------------------------------- */
/* ---------------------------------- */
char Mathematica_eq[] =
"integrate e**(-s*t)*(a*t)**2 dt"
" from t=0 to infinity";
/* ---------------------------------- */
/* ---------------------------------- */