What would be the Julia equivalent for python scipy.optimize.fsolve, to resolve function F (x) f1 = 1- x [1] - x [2] f2 = 8 - x [1] - 3 * x [2] return f1, f2 end x0 = [1,1] Solves a linear system of equations set in not matrix notation

G(x) = [3x[1] + 5x[1] + 2x[3] - 12, 2x[1] + 5x[2] - 4x[3] - 9, 4x[1] - x[2] + 5x[3] + 3] J = ForwardDiff.jacobian(G, [0,0,0]) x = -J \ G([0,0,0]) 7:36:28->>3-element Array{Float64,1}: 2.129032258064516 -1.064516129032258 -2.5161290322580645 The expected results are: x=-33/37 => -0.89189189189…

[image] HerAdri: The expected results are: … Why? G(x) gives you zero and using x = [-33/37, 99/37, 24/37] does not give a zero.

[image] HerAdri: Julia>m = model(GLPK.Optimizer) apologies. It should be Model, not model

In MATLAB, this is done in a single line of code, I also wonder if there is a similar thing in Julia: >> syms x y z >> [x,y,z] = solve(3*x + 5*y + 2*z == 12, 2*x + 5*y - 4*z == 9,4*x - y + 5*z == -3) x = -33/37 y = 99/37 z = 24/37

the result shown in link M \ -b 3-element Array{Float64,1}: -0.8918918918918918 2.675675675675676 0.6486486486486486 The expected results are: x=-33/37 => -0.8918918918918919 y=99/37 => 2.675675675675676 z=24/37 =>0.6486486486486487

You have 5x[1] instead of 5x[2]

sys := { 3x + 5y + 2z = 12, 2x + 5y - 4z = 9, 4x - y + 5z = -3}; [sys := {3x+5y+2z = 12, 2x+5y-4z = 9, 4x-y+5...] solve( sys, {x,y,z}); [{y = 99/37, z = 24/37, x = -33/37}] Maple will automatically use fractions. However, you can force decimal answers using the evalf command. evalf(%); [{…

Julia>m = Model(GLPK.Optimizer) ERROR: MethodError: no method matching Model(::Type{GLPK.Optimizer}) Closest candidates are: Model(::Any, ::Any, ::Any, ::Any, ::Any, ::Any, ::Any, ::Any) at C:\Users\hermesr\.julia\packages\JuMP\iGamg\src\JuMP.jl:126 Model(; caching_mode, solver) at C:\Users\herm…

Yes, in MATLAB, you use double([x y z]) to get decimal answer.

[image] Seif_Shebl: I also wonder if there is a similar thing in Julia There are various packages for symbolic manipulations, but that’s probably overkill for this application. I am curious about the context here. If one is talking about a single linear system, it is trivial to transcribe it…

Julia equivalent for python scipy.optimize.fsolve

General Usage

ChrisRackauckas September 6, 2020, 9:45pm 29

Very very soon.

2 Likes

Topic		Replies	Views
Calling Mathematica into Julia to symbolically solve a system of non-linear equations Specific Domains question , nonlinear , symbolics , mathematica	18	1418	May 11, 2022
How to solve two symbolic nonlinear equations in Julia? Optimization (Mathematical)	27	9535	June 27, 2018
How to solve these equations? General Usage sympy	33	1607	January 12, 2023
Solve system of nonlinear equations New to Julia	12	5930	October 10, 2018
Speed up solution of large system of equations Performance nlsolve , nonlinear , speed-optimization	18	966	August 19, 2023

Julia equivalent for python scipy.optimize.fsolve

Related topics