Advanced Search

Question: Solving numerically a system of non-linear equations

Question:Solving numerically a system of non-linear equations

Delali Accolley 20
solve fsolve
August 31 2015
1

Hi, I'm trying to solve without success numerically the following system of 15 nonlinear equations. Could anyone help, please? Thanks
 

> 

restart
> 

n := 0.27231149e-1:
> 

x := 0.5116034663e-1:
> 

F := .1561816797:
> 

eq1 := sigma*C0 = pgamma*W*H1*(1-E0-L0)/(1+n):
> 

eq2 := sigma*C1 = W*H1*(1-L1):
> 

eq3 := (1+R)*C0 = (1+rho)*exp(x)*C1:
> 

eq4 := (1+R)*C1 = (1+rho)*exp(x)*C2:
> 

eq5 := C1 = (1+phi)*C0:
> 

eq6 := pgamma*L0+pgamma*(1+(1+n)*F/(pgamma*W*H1))*E0+L1 = (1+R)*(1+(1+n)*F/(pgamma*W*H1))/(ppsi*exp(x))-pgamma*(1+(1+n)*F/(pgamma*W*H1))/ppsi:
> 

eq7 := 1 = pgamma*(1+ppsi*E0)/(1+n):
> 

eq8 := exp(x)*A1 = pgamma*W*L0*H1/(1+n)+Epsilon1-C0-F*E0:
> 

eq9 := exp(x)*A2 = W*L1*H1+(1+R)*A1-C1-(1+n)*Epsilon1:
> 

eq10 := (1+R)*A2 = C2:
> 

eq11 := Y = H^alpha*K^(1-alpha):
> 

eq12 := alpha*Y = W*H:
> 

eq13 := (1-alpha)*Y = (1+R)*K:
> 

eq14 := K = A1/(1+n)+A2/(1+n)^2:
> 

eq15 := H = (pgamma*L0+L1)*H1/(1+n):
> 

eq := {eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9, eq10, eq11, eq12, eq13, eq14, eq15}:
> 

vars := {A1, A2, C0, C1, C2, E0, H, H1, K, L0, L1, R, W, Y, Epsilon1}:
> 

NULL
> 

fsolve(eq, vars); 1; assign(%)

fsolve({1 = .6865382886+.1072247031*E0, C1 = 1.475639047*C0, H = .9734907289*(.7052335150*L0+L1)*H1, K = .9734907289*A1+.9476841993*A2, Y = H^.6874443*K^.3125557, (1+R)*A2 = C2, (1+R)*C0 = 1.121850394*C1, (1+R)*C1 = 1.121850394*C2, 1.052491643*A1 = .6865382886*W*L0*H1+Epsilon1-C0-.1561816797*E0, 1.052491643*A2 = W*L1*H1+(1+R)*A1-C1-1.027231149*Epsilon1, 5.171201776*C0 = .6865382886*W*H1*(1-E0-L0), 5.171201776*C1 = W*H1*(1-L1), .3125557*Y = (1+R)*K, .6874443*Y = W*H, .7052335150*L0+.7052335150*(1+.2274915796/(W*H1))*E0+L1 = 6.083468374*(1+R)*(1+.2274915796/(W*H1))-4.515468884-1.027231149/(W*H1)}, {A1, A2, C0, C1, C2, E0, H, H1, K, L0, L1, R, W, Y, Epsilon1})

 (1)
> 

``

 

Download DDGE.mw
Please Wait...