Non-linear system of equations Julia

I'm trying to solve a large number (50) of non-linear simultaneous equations in Julia. For the moment I'm just trying to make this work with 2 equations to get the syntax right etc. However, I've tried a variety of packages/tools - NLsolve, nsolve in SymPy and NLOpt in JuMP (where I ignore the objective function and just enter equality constraints)- without much luck. I guess I should probably focus on making it work in one. I'd appreciate any advice on choice of packages and if possible code.

Here's how I tried to do it in NLsolve (using it in mcpsolve mode so I can impose constraints on the variables I am solving for - x[1] and x[2] - which are unemployment rates and so bounded between zero and 1) :

using Distributions using Devectorize using Distances using StatsBase using NumericExtensions using NLsolve beta = 0.95 xmin= 0.73 xmax = xmin+1 sigma = 0.023 eta = 0.3 delta = 0.01 gamma=0.5 kappa = 1 psi=0.5 ns=50 prod=linspace(xmin,xmax,ns) l1=0.7 l2=0.3 wbar=1 r=((1/beta)-1-1e-6 +delta) ## Test code function f!(x, fvec) ps1= wbar + (kappa*(1-beta*(1-sigma*((1-x[1])/x[1])))) ps2= wbar + (kappa*(1-beta*(1-sigma*((1-x[2])/x[2])))) prod1=prod[1] prod2=prod[50] y1=(1-x[1])*l1 y2=(1-x[2])*l2 M=(((prod1*y1)^((psi-1)/psi))+((prod2*y2)^((psi-1)/psi))) K=((r/eta)^(1/(eta-1)))*M pd1=(1-eta)*(K^eta)*(M^(-eta))*prod1 pd2=(1-eta)*(K^eta)*(M^(-eta))*prod2 fvec[1]=pd1-ps1 fvec[2]=pd2-ps2 end mcpsolve(f!,[0.0,0.0],[1.0,1.0], [ 0.3, 0.3])

I get this error messageNon-linear system of equations Julia

Any suggestions very welcome! I appreciate the formulas are pretty ugly so let me know if any further simplifications helpful (I have tried!)

Category:python Views:1 Time:2019-03-13

Related post

  • How to solve random non-linear system of equations? 2011-11-11

    I used matlab, the built-in solve function is amazing. It takes a set of strings with x,y,z or anything in it, and solve any nonlinear system of equations with almost all possible solutions. Like: x^2 + y^3 = 2 sin(y) * diff(log (x)) = 3 It can solve

  • non linear system of differential equations 2014-11-02

    i'm trying to solve the following non linear system of differential equations but i don't understand where is the problem in the code that i wrote. The system is non linear with complex coefficients and complex domain. from pylab import * from scipy.

  • Application for solving linear system of equations 2009-12-16

    I needed an application for solving linear systems of equations (N up to 10), so I got different codes, and compile them, and they seem to work, but I get lots of problems with precision. I mean, the solvers are really very sensitive to small changes

  • parallel iterative algorithms for solving Linear System of Equations 2012-04-28

    Does someone know any library or ready source code of parallel implementation of quick iterative methods (bicgstab, CG, etc) for solving Linear System of Equations for example using MPI or OpenMP? --------------Solutions------------- PetSC is a good

  • Solving linear systems of equations 2010-07-24

    I have a system of 6 equations that I need to solve over and over again in a program (with many different inputs of course). I am currently using the Cramer's rule method of solving the system and it works quite well (it seems that my processor reall

  • How to solve a linear system in Linux shell? 2010-09-06

    Does anyone know of a Linux command that reads a linear system of equations from its standard input and writes the solution (if exists) in its standard output? I want to do something like this: generate_system | solve_system --------------Solutions--

  • What good libraries are there for solving a system of non-linear equations in C++? 2010-11-20

    In a C++ application I'm coding, I need to solve a system of non-linear equations (N equations, N unknowns). The systems I'm solving will be rather small (up to 10 equations/unknowns), so performance is not going to be a real issue. I've searched the

  • Mathematica: help solving a system of non-linear equations with inequality constraints 2011-05-07

    I have a system of 4 non-linear equations in 4 unknowns. In addition I have 1 inequality constraint that I need a function of the four unknowns to satisfy (and then I need the four unknowns to be non-negative). My problem has 12 or so parameters and

  • Solving system of linear and non linear equations? 2011-08-24

    My system of equations have both linear and non linear equations but not quadratic or cubic or equations having variable degree more than one. example : x + y = 3 (linear), y + z = 4 (linear), x * z = 6 (non linear), x / y = 3 (non linear), y * z / x

  • Specifying positive solutions for a system of non-linear equations 2012-04-09

    I want to solve a system of non-linear equations in matlab such that the solutions are positive. Can I specify in matlab that I want positive solutions for unknown variables? How? --------------Solutions------------- Yes, you can define the variables

  • Python non linear equation with lagrangian multipliers estimation 2011-10-19

    i 've been struggling for days on this thing....but to no avail. I'm not very good at difficult math let alone this kind of level of difficulty. I was trying to implement the maximum entropy application for the lottery in python for my graduation ass

  • Solution of non-linear differential equation 2012-03-05

    I don't use Mathematica in general and I need it to compare with an other program. I want to solve a system of three differential and non linear equations. For this I use Dsolve. Everything goes wrong when I put the nonlinear term (exponential). Here

  • How to solve a symbolic non-linear vector equation? (Matlab or other) 2011-07-13

    I'm trying to find a solution to this symbolic non-linear vector equation: P = a*(V0*t+P0) + b*(V1*t+P1) + (1-a-b)*(V2*t+P2) for a, b and t where P, V0, V1, V2, P0, P1, P2 are known 3d vectors. I attempted to do that in Matlab like this: P = sym('P',

  • Fast linear system solver for D? 2011-02-08

    Where can I get a fast linear system solver written in D? It should be able to take a square matrix A and a vector b and solve the equation Ax = b for b and, ideally, also perform explicit inversion on A. I have one I wrote myself, but it's pretty sl

  • How to draw best fit non linear curve in C# 2011-03-18

    after lots of timing searching on net I come back to Stack overflow.can anyone tell me where i will get code for drawing best fit nonlinear curves in c# or wpf & combine curves of linear curve & non linear curves.Thank you --------------Solut

  • Non-linear regression in C# 2011-10-17

    I'm looking for a way to produce a non-linear (preferably quadratic) curve, based on a 2D data set, for predictive purposes. Right now I'm using my own implementation of ordinary least squares (OLS) to produce a linear trend, but my trends are much m

  • Difference between a linear problem and a non-linear problem? Essence of Dot-Product and Kernel trick 2009-07-18

    The kernel trick maps a non-linear problem into a linear problem. My questions are: 1. What is the main difference between a linear and a non-linear problem? What is the intuition behind the difference of these two classes of problem? And How does ke

  • php - generate non linear urls 2009-12-06

    Update: Ok.. I have read a lot of solutions... Thank you so much everybody. I think I will keep it simple and avoid the encryption and just rely on two fields matching in the database. I can keep the id as it is (sequential) but add for example a tim

  • Java : linear algorithm but non-linear performance drop, where does it come from? 2010-02-16

    I am currently having heavy performance issues with an application I'm developping in natural language processing. Basically, given texts, it gathers various data and does a bit of number crunching. And for every sentence, it does EXACTLY the same. T

  • Solving linear systems in Maxima: how to pass arbitrarily sized arrays into linsolve( ) 2010-02-23

    I have a linear system I'm trying to solve in Maxima using linsolve(eqlist, varlist). For a specific problem, I can use linsolve( [ eq[0],eq[1],eq[2] ], [ a[0],a[1],a[2] ]), and this works. But this is rather ugly since it means I have to manually ad

  • reduce system of equations which is under-determined in mathematica 2011-02-14

    I obtained the folloiwng equation (as an example): {2 w11 + 3 w21 == 2 w12, w11 == 4 w12 + 3 w22, w11 + 2 w21 + w22 == 0, 2 w12 + w21 + 2 w22 == 0} And I want to determine w11, w12, w21, w22. However, simply do the following: Solve[{3 w11 + 2 w21 ==

  • List/matrix of coefficients equation (system of equations) 2011-03-02

    I try to extract coefficient from equations (system of equations) into list (matrix). I have tried CoefficientList[poly, {var1, var2, ...}] but without success. This simple example should explain my problem: Eq1 = a D[U[x1, x2], {x1, 2}] + b D[V[x1,

  • Javascript and php-Solving system of equations with one unknown variable 2011-08-04

    I'm looking for a javascript or php-script that is able to solve a system of equations with one unknown variable x. The equation can be of any kind, i.e. linear, quadratic, cubic etc. e.g. 10-2*x = 25*x or 12*sin(x³) = ln(x²) or e^x² = x^5 I hope I'v

  • How to create a slider with a non-linear scale? 2011-08-30

    I have a slider with a minimum value of 0 and maximum of 500. I want to when the slider goes to 100, the thumb be in the middle of the slider. I know it seems wierd, but some programs do it with zoom slider, and I believe it's better. --------------S

  • Quickselect with array of structs has non linear running time 2011-12-12

    Updated, original question below the line: I need to compute a median, and would like to use the O(N) quickselect algorithm. It turns out however that when the array is no longer a flat array of doubles, but rather an array of structs (of which one e

Copyright (C), All Rights Reserved.

processed in 0.157 (s). 11 q(s)