Tag Info

Hot answers tagged

6

The following remarks are crucial when we are to solve equations symbolically: So far (in general) Mathematica cannot solve transcendental equations when two unknowns are involved, nevervetheless in some exceptional cases it may seem like it could (see e.g. How do I solve this equation?). This is also the case when some symbolic constants are involved ...


5

There are a few problems with your code. Allow me to highlight them and guide you to a solution. First of all, the proper way to define a function in Mathematica is using :=. So your code should read F[x_] := NSum[Exp[-x BesselJZero[0, a]^2]/BesselJZero[0, a]^2, {a, 1, Infinity}] Furthermore, you should note that the zeroes of the Bessel function are ...


5

Define the helper functions and variables like this: F[x_] := {x[[2]] x[[3]], -x[[1]] x[[3]], -0.51 x[[1]] x[[2]]} dt = 0.1; xi = {0, 1, 1}; Do not define xj since you need to use that as a single symbol in FindRoot. Now this equation, xj == xi + dt F[xj] makes sense if we substitute a concrete vector value for xj. If we don't the equation will still ...


5

a = 0.00228797 ; b = 0.00138781/2; c = 0.00209387; d = -0.281261; e = -0.832702; f = -1.43328; ellipse1 = a x^2 + 2 b x y + c y^2 + d x + e y + f; ContourPlot[ellipse1 == 0, {x, -500, 500}, {y, -500, 500}, ImageSize -> 200] uv =First@Solve[ForAll[{x, y}, a (x - u)^2 + 2 b (x - u) (y - v) + c (y - v)^2 + f2 == ellipse1], {u, v, f2}, Reals] ...


4

This is a way to perhaps get approximate solutions and perhaps refine numerically. Using setup: a = 0.25 F = 0.015641 q1 = 5.3415 c[b_] :=Log[b+1] Then f[x_, y_] := (x^2 + q1^2 - y^2)/2*(x + q1 - y) + (1 - 2*a)*x + 2*F - 2*(1 - a)*c[x] g[x_, y_] := (x^2 + q1^2 - y^2)/2*(x + q1 - y) + (1 - 2*a)*y + 2*F - 2*(1 - a)*c[y] p1 = ...


4

Here is some preprocessing that might be useful. We will get rid of radicals "by hand", and eliminate variables using GroebnerBasis. The end result will be a single (complicated) expression in several variables. eqns = {x == (Etot - M - MI)/( Etot + 1/4 (Kd - Sqrt[Kd] Sqrt[8 Etot + Kd])), Etot == 2 Di + M + Inh M, Itot == Inh + MI, Kd == M^2/Di, ...


4

Numerical solution: sol[a_] := {a, α} /. First@NSolve[ Tan[a*Sqrt[α]] + Tan[(1 - a)*Sqrt[α]] == Sqrt[α] && 5 <= α <= 35, α]} // Quiet ListLinePlot[sol /@ Range[0, 0.5, 0.01]] I don't know why ContourPlot doesn't work well,even though I change some options of ContourPlot.Maybe it is because there are some singular points. ...


4

The problem here can arise because of numerical underflow which appears for sufficiently large dimension of the problem. Some numerically very small number multiplies the parameter "a" and therefore "a" does not appear in the "solution". Consider a simple example Define the matrix m (fill it with random numbers, here exponentially distributed) In[263]:= ...


4

The equation is numerically unstable. The second root is not valid at machine precision so is rejected. Higher precision is required to arrive at a valid second root. Use Solve to see exact solutions as root objects. Then use higher precision when converting root objects to numbers. k = .2; eqn = Rationalize[2 r^2 - .2 r^(k + 2) - 1 == 0]; sol = Solve[eqn, ...


3

The Method option in LQOutputRegulatorGains is not documented yet and seems to be going through several iterations. This is what works for the different versions. v8: Method -> {"RiccatiSolveOptions" -> {Method -> "Eigensystem"}} v9: Method -> {"Riccati", {Method-> "Eigensystem"}} v10:Method ->{"Riccati"->{Method-> ...


3

Just a way to visualize by rearranging (noting Chenminqi 's observation). Uses: $tan(x)+tan(y)=sin(x+y)/(cos(x)cos(y))$. Limiting range important to avoid periodicity and singularities: cp = ContourPlot[ Sin[Sqrt[y]]/Sqrt[y] == Cos[(1 - x) Sqrt[y]] Cos[ x Sqrt[y]], {x, 0, 0.5}, {y, 0, 35}, FrameLabel -> {"a", "\[Alpha]"}, BaseStyle -> 20]; ...


3

This gives the solution to the equation: temp = Solve[eqn[x, y, z], D[y[x, z], x]] This extracts the value of the rule: res = D[y[x, z], x] /. temp This evaluates the result at the requested coordinates: res /. {x -> 1, y -> 1, z -> 1} (* {9} *)


3

I post this only for insights into solving this equation (given the discussion) and ambiguity of interpretation...and for fun. I have alread +1 Artes answer. f[x_] := Sqrt[x] + Sqrt[1 - x^2]; g[x_] := Sqrt[2 - 3 x - 4 x^2]; Note the domain that f(x)-g(x) is real valued. dom = Reduce[x > 0 && Abs[x] < 1 && 2 - 3 x - 4 x^2 > 0] ...


3

We can find the minimum-producing argument x as a function of n2 using NDSolve. We need only find an initial minimum and a differential equation for the trace of the minima. In fact, we can integrate the minimum value of eisl and return a single function that yields the minimum for each n2. Given that the parameters b,e0,g,gamma,p,epsilon`, are (constant) ...


3

Log[e,x] is a syntax error. It can be simplified to Log[x]. Experts may provide numerical approaches to solutions. I post this (which is equivalent to my previous post) to provide insights into your equations. a = 0.25; F = 0.015641; q1 = 5.3415; c[b_] := Log[b + 1] f[x_, y_] := (x^2 + q1^2 - y^2)/2*(x + q1 - y) + (1 - 2*a)*x + 2*F - 2*(1 - a)*c[x] g[x_, ...


3

ToRules is useful here. If we manually set a value to the parameter C[1], we can convert the remaining expression to a list of rules similar to what Solve would return: In[5]:= rules = {ToRules[sol /. C[1] -> 0]} Out[5]= {{phi -> -ArcCos[-(Csc[theta]/Sqrt[3])]}, {phi -> ArcCos[-(Csc[theta]/Sqrt[3])]}, {phi -> ...


3

The right-hand side of your equation is asymptotic to zero, so it has no finite solution. It is a good practice to make a plot of expression before trying to solve expression == 0. Plot[-(2/x^3) - 1/x^(5/2) - 1/x^2 + (4 Log[Sqrt[E] x])/x^3, {x, 0, 100}]


3

I post this for completeness (accepting that it was somewhat opaque). In the series of questions it was known that the quadric form was an ellipse (this could be established using determinants as a previous poster has used). Accepting an ellipse quadratic form,e.g. the test one provided ellipse = -1.43328 - 0.281261 x + 0.00228797 x^2 - 0.832702 y + ...


3

Needs["MultivariateStatistics`"] points = {{31.4799, 432.849}, {-195.826, 356.419}, {210.029, 121.779}, {76.1586, -4.47992}, {26.6061, -6.97711}, {160.236, 27.0398}, {203.248, 241.865}, {-225.351, 218.912}, {-159.899, 109.93}, {-106.465, 58.164}, {-1.98952*10^-11, 442.}, {-229.146, 324.489}, {-31.4799, 9.15114}, {195.826, 85.5807}, ...


3

Something little bit different. a = 0.00228797; b = 0.00138781/ 2; c = 0.00209387; d = -0.281261; e = -0.832702; f = -1.43328; ellipse1 = a x^2 + 2 b x y + c y^2 + d x + e y + f; point1 = {x, y} /. Maximize[{y, ellipse1 == 0}, {x, y}][[2]]; point2 = {x, y} /. Maximize[{-y, ellipse1 == 0}, {x, y}][[2]]; center = (point1 + point2)/2; angle = phi /. ...


3

Plot[2 x - Sinh[x], {x, -Pi, Pi}] FindRoot[2 x == Sinh[x], {x, #}] & /@ {-2, 0, 2} {{x -> -2.17732}, {x -> 0.}, {x -> 2.17732}}


3

Please check the result by doing the following eqn = FullSimplify[ 2*a*b*n*c*Q - 3*a*b*n*p*c*Q + a*b*n*p^2*c*Q - 4*a^2*b*n*p*Q^2 - 4*a*b^2*n*p*Q^2 + 6*a^2*b*n*p^2*Q^2 + 6*a*b^2*n*p^2*Q^2 - 2*a^2*b*n*p^3*Q^2 - 2*a*b^2*n*p^3*Q^2 + a^2*b*c*Q^2 + a*b^2*c*Q^2 - 2*a^2*b*p*c*Q^2 - 2*a*b^2*p*c*Q^2 == 0] // TraditionalForm This returns ...


3

Reduce has its own syntax for specifying constrains and domains. Please try this: Reduce[2*a*b*n*c*Q - 3*a*b*n*p*c*Q + a*b*n*p^2*c*Q - 4*a^2*b*n*p*Q^2 - 4*a*b^2*n*p*Q^2 + 6*a^2*b*n*p^2*Q^2 + 6*a*b^2*n*p^2*Q^2 - 2*a^2*b*n*p^3*Q^2 - 2*a*b^2*n*p^3*Q^2 + a^2*b*c*Q^2 + a*b^2*c*Q^2 - 2*a^2*b*p*c*Q^2 - 2*a*b^2*p*c*Q^2 == 0 && 0 < a ...


2

One alternate approach is FindRoot: FindRoot[ -(2/x^3) - 1/x^(5/2) - 1/x^2 + (4 Log[Sqrt[E] x])/x^3, {x, 1}, MaxIterations -> 10000 ] During evaluation of In[16]:= FindRoot::cvmit: Failed to converge to the requested accuracy or precision within 10000 iterations. >> {x -> 6.626338311108174*10^1760} From there I would be encouraged to ...


2

Clear["Global`*"]; f[rx_, ry_] := Module[{d1, d2, d3, a, b, c, sol1, sol2, data}, {d1, d2, d3} = # - {rx, ry} & /@ {{0, -1}, {Sqrt[3]/2, 1/2}, {-Sqrt[3]/2, 1/2}}; {a, b, c} = Norm[#]^-3 & /@ {d1, d2, d3}; sol1 = NSolve[a b (d2 - d1) Sin[{x, y}.(d1 - d2)] + b c (d3 - d2) Sin[{x, y}.(d2 - d3)] + a c (d3 - d1) Sin[{x, y}.(d1 ...


2

iter = Reap[FindRoot[Sin@x == Cos[x], {x, 0}, EvaluationMonitor :> Sow[x]]] {{x -> 0.785398}, {{0., 1., 0.782042, 0.785398, 0.785398}}} Length@Last@Last@iter - 1 4


2

Options[FindRoot] (* {AccuracyGoal -> Automatic, Compiled -> Automatic, DampingFactor -> 1, Evaluated -> True, EvaluationMonitor -> None, Jacobian -> Automatic, MaxIterations -> 100, Method -> Automatic, PrecisionGoal -> Automatic, StepMonitor -> None, WorkingPrecision -> MachinePrecision} *) ?? EvaluationMonitor (* ...


2

This does not solve the equation (in the following I just let $w=1$), of which there are infinite, but I hope provides insights. The surface of solutions is sectioned by desired x value. sol[g_, x0_] := Module[{p, l, sn}, p = g[[1, 1]]; l = Cases[g, Line[w___] :> w, Infinity]; sn = Select[p[[#]] & /@ l, #[[1, 1]] == x0 &]; ...


2

Symbolically solve this equation: Solve[Tan[a Sqrt[α]]+Tan[(1-a) Sqrt[α]]==Sqrt[α]//TrigToExp,a] Solve[Tan[a Sqrt[α]]+Tan[(1-a) Sqrt[α]]==Sqrt[α]/.a->ArcTan[x]/Sqrt[α] //FullSimplify,x] Updated: expr=Tan[a*Sqrt[α]]+Tan[(1-a)*Sqrt[α]]-Sqrt[α]; expr=FullSimplify@ReleaseHold@Numerator@Together[expr/.Tan[x_]:>Sin[x]/Hold@Cos[x]] ...


2

Funny problem. I like it. But there is must be something else in the condition of the problem? Or I misunderstand you. You say that I am guaranteed that w1, w2 exist Ok. Let's set v1 = {0, 0, 0, 1}; v2 = {1, 0, 0, 0}; it's independent vectors. Then find w1 and w2 in general case for this example: w1 = v1*a1 + v2*b1 w2 = v1*a2 + v2*b2 ...



Only top voted, non community-wiki answers of a minimum length are eligible