/* --------------------------------- */
/* save as fc.h */
/* --------------------------------- */
#define LOOP 2*300
/* --------------------------------- */
#define A .3
/* --------------------------------- */
double s = .6;
double aa = A;
/* --------------------------------- */
double F(
double t)
{
return( pow(t,2) );
}
char Feq[] = "t**2";
/* --------------------------------- */
double F_mns_a(
double t)
{
return( pow((t-aa),2) );
}
char F_mns_aeq[] = "(t-a)**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( exp(-aa*s)* 2./pow(s,3));
}
char f_seq[] = "exp(-a s) (2/s^3)";
char f2seq[] = "2 exp(-a s)/s^3";
/* ---------------------------------- */
/* ---------------------------------- */
double b = 300.;
double a = A;
/* ---------------------------------- */
/* ---------------------------------- */
char Mathematica_eq_1[] =
"integrate exp(-s *t) * (t- a)**2 dt"
" from t=a to infinity";
/* --------------------------------- */
char Mathematica_eq_2[] =
"integrate exp(-.6*t) * (t-.3)**2 dt"
" from t=.03 to 300";
/* ---------------------------------- */
/* ---------------------------------- */