Quasi-Newton least squares method: user defined fit

You have the possibility to fit data of your own or let the program generate the data from your specifications. In the first case your input x serves as your initial guess for the optimal parameter and otherwise it is used for the data generation as the unperturbed solution.

Please type here an identification text used for the plot:

Please specify the number of coefficients n here:
Number of fit coefficients n = Important: 1 <= n <= 20 !

Please write a piece of code for the computation of f(x;t) in the following field.
The coefficients appear as x(1) ... x(n) and the ''free'' variable is t:
The input parameters are n, x, t and the result must be named as yt! You may use here prefdefined real variables x1h,x2h,x3h,x4h,x5h,x6h,x7h,x8h,x9h,y(100), the integers i,j,k and logical variables bool1,bool2,bool3. All these variables are initialized with 0 resp. .false. There are also available the constants pi, sqrt2(=sqrt(2.d0)),e1(=exp(1.0d0)). The restrictions of JAKEF apply here!

Initial (for own data) or true vector (for data generation): x(1),...,x(n)
x=

Do you want to fit data of your own or should the program generate data artificially?
internally generated

Number m of data points:
m = Important: n <= m <= 200 !

The interval for t:
a = b =

The perturbation level for the data y:
level = Important: 0 <= level <= 1 !

Own data

Number m of data points:
m = Important: n<= m <= 200 !

Please type your data in the textarea, m pairs of numbers, not necessarily each on a new line
for example:
0,0, 1,1
2,1.4 and so on.

use the standard parameters given below
My own settings:
xtol = terminate if solution x changes less than this (rel)
ftol = terminate if sum of squares changes less than this (rel)
maxstep = max number of steps allowed
maxfev = max number of function evaluations allowed
iprint = -1 = none , 0 = final , 1 = final and short iteration protocol
index1 = index for a component of x to be printed over stepnumber
index2 = index for a second component of x to be printed over stepnumber

Warning!!! - This may take some time.

Click on "evaluate", in order to submit your input.

Back to the theory page

 Back to the top!

28.09.2015