Fitting a conic to a given set of points using Trust Region method

A general conic can be uniquely describe by the following equation up to a scale factor: Ax^2+Bxy+Cy^2+Dx+Ey+F=0Then (A,B,C,D,E,F) is often called algebraic parameter vector of the conic.Usage: [ParA,RSS,iters,Jg] = TR_conic(XY,ParAini,DeltaIni)Child functions:Residuals_ellipse, Residuals_hyperbola, AtoG(can be found from previous submissions) , JmatrixLMA (included in the main function)Input:XY:given points i=1 to nParAini = [A,B,C,D,E,F]'- the initial parameter vectorDeltaIni: the initial size of the trust region.(this is optional; if it is missing, TR sets it to 1)Output:ParA: vector of algebraic parameters of the conic foundRSS: the Residual Sum of Squares (the sum of squares of the distances)iters: the number of (main) iterationsJg: the norm of the first derivative