Consider two system of nonlinear equations. In this case we will get F(x1, x2)=[c1 c2] ^T; Based on my problem, I need to calculate inverse of F. Also, I have try pinv(F), inv(F) but code gives
............................................
??? Error using ==> sym.svd
Too many input arguments.
Error in ==> pinv at 29
[U,S,V] = svd(A,0);
Error in ==> sys_4th_kt at 48
t=f(y)*pinv(f(x0));
............................
The following one:
ux = Df(x0)\f(x0); %System of linear equation Ax = b, x = inv(A)*b, but better way is x=A\b.
y = x0 - (1)*ux;
t=f(y)*pinv(f(x0));
p=Df(x0)*(I-t)^2;
ux1 = f(y)*pinv(p);
x1 = y-ux1; % 4th order Kung and Traub
You can use Inverse[A] command for finding inverse of A in Mathematica.
The function for computing the pseudo-inverse is called pinv in Matlab, which you have already tried. It requires one matrix as input. From the error message I would guess, your function f(x0) does not give a matrix.
The inverse of a matrix is possible only if it is a square matrix and its determinant should be different from zero. You are finding the impossible inverse.
Dr. Hauke, yes your right. I think mathematica can solve this, so I am trying to right code in Mathematica. Also, tell me any other way in matlab
A simple procedure would be to linearize your F(X) using Taylor expansion and then, having a "good" initial vector X0, you would search for solution for your system, very similar to Newton-Raphson approach. I like Vanicek presentation here: http://www2.unb.ca/gge/Pubs/LN35.pdf
- Badges
- Science topic