|
Specification of a function of your own:
Please type the evaluation program of your function here
using FORTRAN rules.
Your final statement must be
fu= some expression you computed before or just here depending on x
You may use the constants pi, e(=exp(1)), sqrt2(=1.414...),
the integer variables i,j,k, the logicals bool1,bool2,bool3
and the double precision variables sum,h1,h2,h3,h4,y(100),z(100),a(100,100)
which are all intialized with zero resp. .false. .
The routine has the parameters x
(double, input) and fu (double out). never change x!.
first is a local integer and set 0 before calling the function the
first time.
You may use this in order to initialize some local data and set it 1
afterwards to avoid multiple such initialization. Your settings of the local
variables are preserved during program execution.
if you want to compute f'(a), f'(b)analytically
your piece of code must follow the restrictions of JAKEF
|