The content of my post has bin lost apparently... 

This is my Maple-input:

>restart;
>BondMarket:= B(i, rho[b], y, epsilon[S, B]) - beta(i, rho[b], epsilon[D, beta])*(1-tau)*T(i, y, epsilon[D, M]);
>LoanMarket:= L(rho[l], i, y, epsilon[D, L]) - lambda(rho[l], i, epsilon[S, lambda])*(1-tau)*T(i, y, epsilon[D, M]);
>MoneyMarket:= T(i, y, epsilon[D, M]) - m(i)*R;
>GoodsMarket:= y - Y(i, rho[b], rho[l]);

>LogBondMarket:= log(BondMarket); 
>LogLoanMarket:= log(LoanMarket); 
>LogMoneyMarket:= log(MoneyMarket); 
>LogGoodsMarket:= log(GoodsMarket);

>TaylorLogBondMarket:= mtaylor(LogBondMarket, [i=i[e], rho[b]=rho[b][e], y=y[e], epsilon[S, B]=epsilon[S, B][e], epsilon[D, beta]=epsilon[D, beta][e], epsilon[D, M]=epsilon[D, M][e]], 2); 
>TaylorLogLoanMarket:= mtaylor(LogLoanMarket, [rho[l]=rho[l][e], i=i[e], y=y[e], epsilon[D, L]=epsilon[D, L][e], epsilon[S, lambda]=epsilon[S, lambda][e], epsilon[D, M]=epsilon[D, M][e]], 2); 
>TaylorLogMoneyMarket:= mtaylor(LogMoneyMarket, [i=i[e], y=y[e], epsilon[D, M]=epsilon[D, M][e]], 2); 
>TaylorLogGoodsMarket:= mtaylor(LogGoodsMarket, [y=y[e], i=i[e], rho[b]=rho_[b][e], rho[l]=rho[l][e]], 2);

I want to solve these 4 equations with a gauss matrix for the variables i, rho[b], rho[l], y, epsilon... It's okay that the set does not have one unique solution (less expressions than vars). 

Can somebody help me?

Thanks!

S.