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