C_R

3412 Reputation

21 Badges

5 years, 315 days

MaplePrimes Activity


These are replies submitted by C_R

The problem (for choosen numerical parameters) is ill-conditionned because it can be reduced to

which has no solution appart from the a trivial one where a__1=0 and c__1=0 (which solve cannot find).
However, the way to solve, that I have provided below, should give solutions for other numerical parameters. Here is an updatad version using exact numbers
nonlin_sys.mw

@GFY

Yes, that is what I meant by: "For equation 3 you get in the same way 2 more equations sets of solutions".

These solutions are only valid for a subset of equations: equation 3, and equation 2, and equation 4.

I did not check if there is in all subsets of solutions a theta_3 that fullfils equation 1 and 3 at the same time.

We are dealing here with a complex argument which is beyond my comfort zone

FunctionAdvisor(plot,tan)

It seems that the real part of the two different solutions for theta__3 is Pi/2 spaced appart.

@acer 

I would not call it a bug. It comes with the design of the labels. Can't the labels be turned on permanetly in help pages. Its essential information when referred to. I use them allot in techincal documentation.

I cannot reprocude the mixed input of 1D Math and 2D Math.

I did in a rather new Maple session (Maple started in 2D):

?D  -> changed to 1D view -> copied the input -> new worksheet -> pasted the input

@dharr 

I will pay attention to this in the next update.

@Axel Vogt 

I could fix it. Thanks

 @rlopez @Carl Love @Ronan 

I had two Maple sessions running for 2 weeks. About a week ago, I noticed that the submenus had reappeared and I could not trigger the disappearance as described above.

Yesterday I restarted the laptop and the first application I started was Maple 2023 followed with a bit delay Maple 2024 (see startup times below). Submenues disappread in both Maple instances, as expected, after the system went to hibernation. After starting a pdf reader (Foxit), the submenus were back.

I could repeatedly toggle submenu dissaperance by: closing Foxit -> Hibernation -> wakeup (&checking if submenus were gone) -> starting Foxit (&checking if submenus were back).

After system wakeup today: No submenues. Toggling as described above did not work. After about an hours submenus were back in the Maple 2023 session. The only activities I can remember were browsing in Firefox, Outlook and web applications running in Microsoft edge.

After about an hour the Submenus appeard again in the Maple2024 session.
These are the tasks running in explorer.exe and their startup times (Maple 2023 java task highlighted in grey)

 

I have now restarted my other Laptop with one Maple 2024 Session. After the first hibernation submenus disappeard. I could not get the menues back as described above.
After writing this post I checked again and ... the menus are back.

I will try to repeat this. This time without launching any other programs appart from Maple 2024.

 

@Christian Wolinski 

Chrome uses information from previous user entries (on all web forms not only MaplePrimes) to help you filling in, autocomplete or correct new entries in web forms. Chrome tries to guess what you want to enter. Whether this happens with-in Chrome or accros applications is difficult to say. The tendency is clear. With AI assistance tools all user actions (including old mails and documents) are targeted to be analysed to make our live easier for "free" (there is no free lunch).

I use mostly Firefox and have so far not seen login credentials suggested in searches.

What is wrong here? I cannot find this package.

@GFY 


My initial guess that theta_3 and theta_4 are free to chose was incroeect. I striked it through in teh above. What describes it better: Theta__3 does not depend on theta__4. I.e., for a valid theta_3 there are two valid values for theta__4.

When you plug in the above 16 solutions for a_1 and c_1 into the 4 equations and solve for theta_3 and theta_4, you will get all combinations of variables (in total 64) that are solutions for the system of equations

Here is the first solution for theta__3 from equation 1

subs(a__1 = 5.454567591 + 0.005714161093*I, c__1 = 0.002689473223 - 3.665495195*I, secular7[1]);
solve(%, [theta__3]);
      /               12                 9  \            
      \-1.888144781 10   - 1.978005270 10  I/ ((0.0063172

         - 0.1217592075 I) cos(theta__3)

         + (0.121759209 + 0.006316869337 I) sin(theta__3))


          [[theta__3 = -0.8090474571 + 6.755468571 I]]

which combines with two solutions for theta_4 from equations 2 and 4 to

[a__1 = 5.454567591 + 0.005714161093*I, c__1 = 0.002689473223 - 3.665495195*I],[theta__3 = -0.8090474571 + 6.755468571*I],[theta__4 = -0.7836748989 + 8.209192730*I];
[a__1 = 5.454567591 + 0.005714161093*I, c__1 = 0.002689473223 - 3.665495195*I],[theta__3 = -0.8090474571 + 6.755468571*I],[[theta__4 = 0.7871214279 + 8.209192730*I]]

For equation 3 you get in the same way 2 more equations sets of solutions

@NIMA112 

Here is the idea for H[1] which computes much faster. Carefully check if the expressions were copied correctly. C must be part of the sum since it contains n. Check if that is really what you want. The resutling expression is still big.

ABC.mw

No error with 2023.2

@GFY 

Below is a way to solve for the 16 roots. As said above theta_3 and theta_4 are free to choose. If you replace the floats by parameters or rational numbers, this solution will be exact.
 

restart

secular7 := -(3.461584716*10^11)*a__1*((1.024337843*a__1^4+(4.073013529*c__1^2+29.75257330)*a__1^2+c__1^4+25.6247057300000*c__1^2-9.31744223200000*10^(-9))*cos(`θ__3`)+(.1403469826*a__1^2+.1233422711*c__1^2-2.396674825)*sin(`θ__3`)), -(3.752578062*10^10)*((1.000000000*c__1^4+(12.21904058*a__1^2+13.28028686)*c__1^2-110.1848865*a__1^2+9.219040578*a__1^4+1.057632194*10^(-8))*cos(`θ__4`)+(-.1226429622*c__1^2-.5582050515*a__1^2+4.766172985)*sin(`θ__4`))*c__1, -(3.461584716*10^11)*a__1*((1.024337843*a__1^4+(4.073013529*c__1^2+29.75257330)*a__1^2+c__1^4+25.6247057300000*c__1^2-9.31744223200000*10^(-9))*sin(`θ__3`)+(-.1403469826*a__1^2-.1233422711*c__1^2+2.396674825)*cos(`θ__3`)), -(3.752578062*10^10)*((1.000000000*c__1^4+(12.21904058*a__1^2+13.28028686)*c__1^2-110.1848865*a__1^2+9.219040578*a__1^4+1.057632194*10^(-8))*sin(`θ__4`)+(.1226429622*c__1^2+.5582050515*a__1^2-4.766172985)*cos(`θ__4`))*c__1

-0.3461584716e12*a__1*((1.024337843*a__1^4+(4.073013529*c__1^2+29.75257330)*a__1^2+c__1^4+25.6247057300000*c__1^2-0.9317442232e-8)*cos(theta__3)+(.1403469826*a__1^2+.1233422711*c__1^2-2.396674825)*sin(theta__3)), -0.3752578062e11*((1.000000000*c__1^4+(12.21904058*a__1^2+13.28028686)*c__1^2-110.1848865*a__1^2+9.219040578*a__1^4+0.1057632194e-7)*cos(theta__4)+(-.1226429622*c__1^2-.5582050515*a__1^2+4.766172985)*sin(theta__4))*c__1, -0.3461584716e12*a__1*((1.024337843*a__1^4+(4.073013529*c__1^2+29.75257330)*a__1^2+c__1^4+25.6247057300000*c__1^2-0.9317442232e-8)*sin(theta__3)+(-.1403469826*a__1^2-.1233422711*c__1^2+2.396674825)*cos(theta__3)), -0.3752578062e11*((1.000000000*c__1^4+(12.21904058*a__1^2+13.28028686)*c__1^2-110.1848865*a__1^2+9.219040578*a__1^4+0.1057632194e-7)*sin(theta__4)+(.1226429622*c__1^2+.5582050515*a__1^2-4.766172985)*cos(theta__4))*c__1

(1)

indets([secular7])

{a__1, c__1, theta__3, theta__4, cos(theta__3), cos(theta__4), sin(theta__3), sin(theta__4)}

(2)

Removing theta__3

secular7[1]/cos(theta__3):
convert(expand(%),tan):
solve(%,{tan(theta__3)})[];

tan(theta__3) = -0.2000000000e-6*(0.1772916110e18*a__1^4+0.7049540690e18*a__1^2*c__1^2+0.1730792358e18*c__1^4+0.5149552650e19*a__1^2+0.4435104486e19*c__1^2-1612655781.)/(4858229699.*a__1^2+4269597205.*c__1^2-0.8296292943e11)

(3)

secular7[3]/cos(theta__3):
convert(expand(%),tan):
solve(%,{tan(theta__3)})[];

tan(theta__3) = 5000000.*(4858229699.*a__1^2+4269597205.*c__1^2-0.8296292943e11)/(0.1772916110e18*a__1^4+0.7049540690e18*a__1^2*c__1^2+0.1730792358e18*c__1^4+0.5149552650e19*a__1^2+0.4435104486e19*c__1^2-1612655781.)

(4)

(3)/(4)

1 = -0.4000000000e-13*(0.1772916110e18*a__1^4+0.7049540690e18*a__1^2*c__1^2+0.1730792358e18*c__1^4+0.5149552650e19*a__1^2+0.4435104486e19*c__1^2-1612655781.)^2/(4858229699.*a__1^2+4269597205.*c__1^2-0.8296292943e11)^2

(5)

(5)*(denom@rhs)((5))

(4858229699.*a__1^2+4269597205.*c__1^2-0.8296292943e11)^2 = -0.4000000000e-13*(0.1772916110e18*a__1^4+0.7049540690e18*a__1^2*c__1^2+0.1730792358e18*c__1^4+0.5149552650e19*a__1^2+0.4435104486e19*c__1^2-1612655781.)^2

(6)

((x -> x)-rhs)((5)*(denom@rhs)((5)));# removing denominators and bringing everything to the lhs

(4858229699.*a__1^2+4269597205.*c__1^2-0.8296292943e11)^2+0.4000000000e-13*(0.1772916110e18*a__1^4+0.7049540690e18*a__1^2*c__1^2+0.1730792358e18*c__1^4+0.5149552650e19*a__1^2+0.4435104486e19*c__1^2-1612655781.)^2 = 0.

(7)

 

Removing theta__4

secular7[2]/cos(theta__4):
convert(expand(%),tan):
solve(%,{tan(theta__4)})[];

tan(theta__4) = 0.5000000000e-7*(0.3459516943e19*a__1^4+0.4585290362e19*a__1^2*c__1^2+0.3752578062e18*c__1^4-0.4134773878e20*a__1^2+0.4983531313e19*c__1^2+3968847369.)/(0.1047354015e11*a__1^2+2301136447.*c__1^2-0.8942718090e11)

(8)

secular7[4]/cos(theta__4):
convert(expand(%),tan):
solve(%,{tan(theta__4)})[];

tan(theta__4) = -20000000.*(0.1047354015e11*a__1^2+2301136447.*c__1^2-0.8942718090e11)/(0.3459516943e19*a__1^4+0.4585290362e19*a__1^2*c__1^2+0.3752578062e18*c__1^4-0.4134773878e20*a__1^2+0.4983531313e19*c__1^2+3968847369.)

(9)

(8)/(9)

1 = -0.2500000000e-14*(0.3459516943e19*a__1^4+0.4585290362e19*a__1^2*c__1^2+0.3752578062e18*c__1^4-0.4134773878e20*a__1^2+0.4983531313e19*c__1^2+3968847369.)^2/(0.1047354015e11*a__1^2+2301136447.*c__1^2-0.8942718090e11)^2

(10)

((x -> x)-rhs)((10)*(denom@rhs)((10)));# removing denominators and bringing everything to the lhs

(0.1047354015e11*a__1^2+2301136447.*c__1^2-0.8942718090e11)^2+0.2500000000e-14*(0.3459516943e19*a__1^4+0.4585290362e19*a__1^2*c__1^2+0.3752578062e18*c__1^4-0.4134773878e20*a__1^2+0.4983531313e19*c__1^2+3968847369.)^2 = 0.

(11)

Reducing the order of polynominals

a__1=sqrt(A),c__1=sqrt(C)

a__1 = A^(1/2), c__1 = C^(1/2)

(12)

subs((12),[(7),(11)])

 

[(4858229699.*A+4269597205.*C-0.8296292943e11)^2+0.4000000000e-13*(0.1772916110e18*A^2+0.7049540690e18*A*C+0.1730792358e18*C^2+0.5149552650e19*A+0.4435104486e19*C-1612655781.)^2 = 0., (0.1047354015e11*A+2301136447.*C-0.8942718090e11)^2+0.2500000000e-14*(0.3459516943e19*A^2+0.4585290362e19*A*C+0.3752578062e18*C^2-0.4134773878e20*A+0.4983531313e19*C+3968847369.)^2 = 0.]

(13)

Solve

solve(%)

{A = 29.75227495+0.6233655581e-1*I, C = -13.43584779-0.1971650235e-1*I}, {A = .9664284322+0.4214874164e-1*I, C = -28.52004755-.3397994675*I}, {A = -0.1330420049e-4+0.4783720471e-1*I, C = 0.8966865280e-3+0.3797299034e-1*I}, {A = -10.10493915+0.2865831138e-1*I, C = 23.77838669-.1223123774*I}, {A = -10.10493915-0.2865831138e-1*I, C = 23.77838669+.1223123774*I}, {A = -0.1330420049e-4-0.4783720471e-1*I, C = 0.8966865280e-3-0.3797299034e-1*I}, {A = .9664284322-0.4214874164e-1*I, C = -28.52004755+.3397994675*I}, {A = 29.75227495-0.6233655581e-1*I, C = -13.43584779+0.1971650235e-1*I}, {A = 29.75248435+0.5813738538e-1*I, C = -13.43601042-0.1642603678e-1*I}, {A = .9662748803+0.2933062771e-2*I, C = -28.51938427+.2012880136*I}, {A = 0.2187606320e-3+0.2805489963e-1*I, C = 0.2889472001e-3-.1260960552*I}, {A = -10.10522706+0.4580686029e-1*I, C = 23.77849379-.1239405402*I}, {A = -10.10522706-0.4580686029e-1*I, C = 23.77849379+.1239405402*I}, {A = 0.2187606320e-3-0.2805489963e-1*I, C = 0.2889472001e-3+.1260960552*I}, {A = .9662748803-0.2933062771e-2*I, C = -28.51938427-.2012880136*I}, {A = 29.75248435-0.5813738538e-1*I, C = -13.43601042+0.1642603678e-1*I}

(14)

for i from 1 to nops([(14)]) do eval([(12)],(14)[i]) end do;

 

[a__1 = 5.454567591+0.5714161093e-2*I, c__1 = 0.2689473223e-2-3.665495195*I]

 

[a__1 = .9833045159+0.2143219163e-1*I, c__1 = 0.3181338417e-1-5.340511178*I]

 

[a__1 = .1546348964+.1546779085*I, c__1 = .1394278711+.1361743174*I]

 

[a__1 = 0.4507682445e-2+3.178829890*I, c__1 = 4.876324843-0.1254145092e-1*I]

 

[a__1 = 0.4507682445e-2-3.178829890*I, c__1 = 4.876324843+0.1254145092e-1*I]

 

[a__1 = .1546348964-.1546779085*I, c__1 = .1394278711-.1361743174*I]

 

[a__1 = .9833045159-0.2143219163e-1*I, c__1 = 0.3181338417e-1+5.340511178*I]

 

[a__1 = 5.454567591-0.5714161093e-2*I, c__1 = 0.2689473223e-2+3.665495195*I]

 

[a__1 = 5.454586396+0.5329220326e-2*I, c__1 = 0.2240616595e-2-3.665517077*I]

 

[a__1 = .9829939502+0.1491902758e-2*I, c__1 = 0.1884582447e-1+5.340387573*I]

 

[a__1 = .1189001959+.1179766754*I, c__1 = .2513815163-.2508061393*I]

 

[a__1 = 0.7204873763e-2+3.178880144*I, c__1 = 4.876336257-0.1270836686e-1*I]

 

[a__1 = 0.7204873763e-2-3.178880144*I, c__1 = 4.876336257+0.1270836686e-1*I]

 

[a__1 = .1189001959-.1179766754*I, c__1 = .2513815163+.2508061393*I]

 

[a__1 = .9829939502-0.1491902758e-2*I, c__1 = 0.1884582447e-1-5.340387573*I]

 

[a__1 = 5.454586396-0.5329220326e-2*I, c__1 = 0.2240616595e-2+3.665517077*I]

(15)

NULL

NULL


 

Download solve821_reply.mw

Have you checked for backups?
Depending on how Maple is terminated there should be backups if the options are configured like this

I thought when odetest returns zero without assumptions then the tested solution is valid over the complex domain. In your example this is not the case.

The worring consequence is that plotting for a given parameter the solution of one of the "root ode"s (where the residual of the identical solution of the origianl ode is non zero) does not show an empty plot. For this reason, I am now also of the opinion that one should work with the odetest of the original ode.  There is a flaw in my reasoning that no assumptions are made when using PDEtools:-Solve that I would like to understand.

First 18 19 20 21 22 23 24 Last Page 20 of 67