F=exp[a*(N-Z+1)/A]*b*(1+A^1/3)^2, How we can determine the coefficients a, b, by using MATLAB or other program, the data of F, S, A it in the attached file?
Be careful when doing transformations, the solution given above gives the coefficients for the line closest to Y= aX+lnb, which does not necessarily give the function with the smallest error in F for the stated form of F. Squared error sum in F for the given points in this solution is >100000, closer solutions can be readily found manually. a=-26, b=16 is pretty close; sumsq = 10950,4.
Sorry, I am no MATLAB expert, but defining the sum of the error squared from the four arrays of data should not be very difficult, it's easy in EXCEL, finding the (a,b) point that minimizes this value is probably just some solver. fminunc looks like a good candidate.
You can use Genetic algorithm or other optimization method such as PSO or simulated annealing, By choosing a and b as your optimization variable and optimizing the difference between new generated path and the desired one (F) to zero.