Matlab - 1D inverse modelling

Hola,
he leído un estudio en el que utilizan algo llamado "1D inverse modelling" y me interesa aplicarlo a mis datos.

Si conozco Ai, h y Ci. Pero desconozco z*, p* y q*.
A1=(C1(z*+h)-C1(h))p*+(1-C1(z*+h))q*;
A2=(C2(z*+h)-C2(h))p*+(1-C2(z*+h))q*;
A3=(C3(z*+h)-C3(h))p*+(1-C3(z*+h))q*;
A4=(C4(z*+h)-C4(h))p*+(1-C4(z*+h))q*;

Teniendo en cuenta que Ci es una función de z*, por ejemplo:
C1=1-(4(((z*)^2)/(r^2))+1)^(-0.5)
siendo 'r' conocido.

En el estudio dicen que lo resuelven con 'levenberg marquardt' en matlab ¿Podéis ayudarme? Estoy completamente perdida.
Gracias, LP.

Es decir se trata de un sistema de ecuaciones no-lineales, que necesito resolver en 'matlab', mediante la función 'fsolve' y utilizando el algoritmo de 'Levenberg-Marquard'

Gracias! y please, animaros a responder!
LP

Hola,

he cambiado de nombre xq había alguien que ya se apodaba LP.

LauPra.