Quantcast
Channel: MaplePrimes - Questions and Posts tagged with solve
Viewing all articles
Browse latest Browse all 1098
↧

Solving numerically a system of non-linear equations

August 30, 2015, 9:14 pm
$
0
0

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

↧
Search
Viewing all articles
Browse latest Browse all 1098

Trending Articles

Bath man appears in court charged with attempted murder of a man...

March 16, 2015, 7:37 am

MACLEAN, Allan

July 30, 2019, 6:00 am

Black Angus Grilled Artichokes

July 16, 2016, 4:37 pm

Practice Sheet of Right form of verbs for HSC Students

September 22, 2019, 11:40 pm

Police blotter for Jan. 12

January 12, 2018, 3:30 am

99 God Status for Whatsapp, Facebook

June 5, 2016, 11:46 pm

Rajasthan Board 12th Science Result 2018 name wise- RBSE 12th commerce result...

May 26, 2018, 9:35 pm

Notorious Naushad of Ippa gang nabbed

July 19, 2019, 6:37 am

Child Kidnapping: Amy McNeil was kidnapped on her way to school by 5 adults;...

February 5, 2017, 10:40 am

Sonible Smartlimit v1.1.5-R2R

April 16, 2024, 7:10 am

NCERT Solutions for Class 9th Sanskrit Chapter 3 पाथेयम्

December 22, 2016, 3:50 am

मतलबी दोस्त स्टेट्स | Matlabi Dost Status in Hindi – Selfish Friends Status

February 13, 2020, 3:12 am

Arrow Flash 2 – Sinhala Dubbed – Episode 23 – 20th March 2016

March 20, 2016, 9:39 am

[GET] AI Traffic Goldmine

July 6, 2025, 4:23 am

[E² Plugin] HDF-Radio

January 26, 2025, 9:02 am

Universal Multi-Patch v1.3 By RADIXX11

January 29, 2018, 2:45 pm

IWAN – Thanks and Praise ( Throw Back Thursday )

March 9, 2016, 11:43 pm

RONALD P SONDERGAARD Arrested by Miami-Dade County Corrections on Mar 03, 2017

March 3, 2017, 6:25 am

मुख मैथुन से उठाएं सेक्स का भरपूर मज़ा, जानें क्या है इसका सही तरीकामुख मैथुन...

May 17, 2020, 2:04 pm

HSSC Excise & Taxation Inspector Result 2017 Scorecard/ Category Wise Merit List

July 29, 2017, 2:44 am
Search