Maple solve or solve with matrix

Maple solve or solve with matrix - Messages

#1 Posted: 4/13/2016 3:13:53 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Doing something wrong ...
Trying to solve empirical formula with matrix input using either solve or maple(solve)Solve Block.sm (12 KiB) downloaded 68 time(s).
#2 Posted: 4/13/2016 7:10:44 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Before applying some suite, you must have something to process.
With only two temp. BTU means nothing. Two values are missing.
With the two you provide: who knows to what BTU rate they apply ?

Jean

Forum BTU.gif
#3 Posted: 4/18/2016 3:05:15 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Wrote

Before applying some suite, you must have something to process.

I do. Based on empirical formula, I want to solve for the "flow" based on the total energy. Delta T is constant, so I don't understand your comment.

Here is a revised file, hopefully clearer as to what I am looking for...
Solve Block.sm (7 KiB) downloaded 56 time(s).
#4 Posted: 4/18/2016 5:22:30 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Maybe like this ?

Forum BTU Solved.gif
#5 Posted: 4/18/2016 5:44:16 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Jean,
The first descriptor shows ME doing the algebraic manipulation and getting the correct answer. I am trying to use "solve" to solve the equation with no success. As always, thank you for looking and responding
#6 Posted: 4/18/2016 6:19:49 PM
Mike Kaganski

Mike Kaganski

184 likes in 434 posts.

Group: User

I'm not sure why would you need the solve() when you may directly calculate the result as you did,
but here it goes Solve Block.sm (7 KiB) downloaded 52 time(s).
solve() cannot find units, so you need to "help" it allowing it to search for unitless values and adding required units by hands. Also, it cannot search matrices - so iterate it (I used vectorize() ).
Third - you need to pay attention when you use assignments, and when - comparisons. You happened to use := in solve() instead of ≡.
Fourth - solve(#,#) searches inside the range that is set in program's settings - use solve(#,#,#,#) to define custom range
С уважением, Михаил Каганский
1 users liked this post
ElSid 4/18/2016 6:31:00 PM
#7 Posted: 4/18/2016 6:37:29 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Thank you!
See attached.Solve.PNG
I like documenting and showing the empirical equations for others to "see", whenever possible. Yes, easy for me to re-arrange, put I wanted to keep the look of the empirical equations when possible. This is one for Jean and his units argument.
Why do we need to put the range? (solve(4)) I thought Solve two went only on the positive side so range was not needed. I guess I should de-fault to more instead of less. I had something similar without the units corrected as you did it...at one point
#8 Posted: 4/18/2016 7:20:15 PM
Davide Carpi

Davide Carpi

1416 likes in 2873 posts.

Group: Moderator

Wrote

Why do we need to put the range? (solve(4))I thought Solve two went only on the positive side so range was not needed.I guess I should de-fault to more instead of less.I had something similar without the units corrected as you did it...at one point


solve(2) is perfect for calculations "on the fly"; if you need consistent and reproducible calculations I strongly reccomend solve(4); this because solve(2) it depends on settings that may change from computer to computer or in the same computer they may changes on the time (since the range you need in some application may not be the same of another, and someone - even you - may have changed these settings since last time you have executed the worksheet). However if you still want to use solve(2), I suggest you to document somewhere near the function the range to use.
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
1 users liked this post
Mike Kaganski 4/18/2016 8:09:00 PM
#9 Posted: 4/18/2016 8:34:44 PM
Mike Kaganski

Mike Kaganski

184 likes in 434 posts.

Group: User

I want to add to Davide's excellent explanation.

It's not always good to set solve range (in settings, or in solve(4) ) to very large size, because it may lead to inability to find roots, or poor precision.
solve() first splits the range to 100 even pieces and evaluates the expression for each point (excluding high edge, btw). You may see this if you try solve(x≡x;x) for default range, or similar with solve(4). After that, slove analyzes the data, using limited number of iterations. So, for very large range, no points may help it to find roots, or too many iterations may be required, so no results may be returned.
С уважением, Михаил Каганский
1 users liked this post
Davide Carpi 4/18/2016 8:52:00 PM
#10 Posted: 4/18/2016 11:02:28 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Elsid,

Revisit "Samples" the solve stuff I posted within last 2 weeks.
"solve", "Solve" are brute force of unpredictable inaccuracy.
In one of the examples, the acuracy is e-9, that is enough to
shoot Jupiter instead of Mars. OK, HVAC is empirical. Further
temperature controlled, not to avail so much "pain & prestige".

My Smath 5346 does not have the vectorise operator, survive w/o.
BTW, such huge numbers with so many decimals, make it integer,
and maybe set the last integer value to 0 , to meet empirical.

Cheers, Jean
1 users liked this post
ElSid 4/19/2016 12:24:00 AM
#11 Posted: 4/19/2016 12:24:48 AM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Agree about truncating and rounding ... another reason I am not high precision. Looking up that "brute force" discussion
#12 Posted: 8/3/2017 3:39:08 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Follow up on this post. On 0.98.6398, I am now seeing a problem where assigning the solve block to a variable fails. solve block 2.smz (2 KiB) downloaded 34 time(s).
#13 Posted: 8/3/2017 4:19:33 PM
Martin Kraska

Martin Kraska

1222 likes in 2150 posts.

Group: Moderator

Solve() from the Maxima plugin
2017-08-03 21_14_16-SMath Studio 0.98.6398 - [solve%20block%202.smz_].png
solve%20block%202.sm (4 KiB) downloaded 37 time(s).
Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
#14 Posted: 8/3/2017 5:04:27 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

mkraska - Thank you, but not what I was looking for. I have the power already calculated/known ... I am looking to solve for a rate of time. As you can see in my example, the program will solve as long as I do NOT define the solve block with a variable. It also looks like you are using the Maple plugin.

Update, I thought that Vectorize eliminated the nesting issue. I think that is what is going on as my solution may be a nested item and not a matrix. Definitely a nesting issue. How do I "un-nest" a Vectorized equation? Solve Block 2.smz (2 KiB) downloaded 34 time(s).without jMG Matrix manipulation. I was hoping Vectorize would simplify the issue.
#15 Posted: 8/4/2017 10:55:57 AM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

The solve block is not compatible for assigning objects that are vectorized.
This was first noticed solving Gamma(x), that solves un-assigned.
erf(x), Gamma(x) are two reported, that solve from vectorising, not assignable.
Nothing to do => copy/pate ... Does not spool to file either

Jean
1 users liked this post
ElSid 8/4/2017 2:50:00 PM
#16 Posted: 8/4/2017 2:52:58 PM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

Wrote

The solve block is not compatible for assigning objects that are vectorized....


Any recommendations? For now, I will do copy and paste, but I would like to "automate" the solution since these numbers are used downstream. I'm not seeing a programmed solve program, unless there is an example I missed in the forum.

#17 Posted: 8/4/2017 3:00:40 PM
Martin Kraska

Martin Kraska

1222 likes in 2150 posts.

Group: Moderator

Wrote

Any recommendations? For now, I will do copy and paste, but I would like to "automate" the solution since these numbers are used downstream. I'm not seeing a programmed solve program, unless there is an example I missed in the forum.



Maybe, I don't understand the question. To me, it looks like you try to repeatedly solve an equation for various values of the constants in the equation. And you have vectors of these constants.

So what I did was to generate a symbolic solution and use vectorize to apply the various input values. But you didn't seem to like it.

Also, I am a bit helpless faced with your use of units. I'd use units only in values of quantities, never as free parts of equations, otherwise these equations become dependent on the units used. My approach shown above is completely consistent in terms of units.
Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
#18 Posted: 8/4/2017 9:29:59 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Wrote

Any recommendations?For now, I will do copy and paste, but I would like to "automate" the solution since these numbers are used downstream.I'm not seeing a programmed solve program, unless there is an example I missed in the forum.


==========================================

You have missed nothing from the forum, has never bee exemplified.
I don't work with units, sorry for your work sheet.
As explained before => Gamma(x) is one that needs be vectorised locally
That does not exclude global vectorising for further use, downstream the document.

Cheers Esild ... Jean

SolveVectorise.PNG
#19 Posted: 8/4/2017 10:29:03 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

That' the best I can do.
1. Subscript should not exist [complicate matters].
2. What are "Units" doing in Engineering works ?

Jean

solve block 2 Examine Solutions.sm (7 KiB) downloaded 40 time(s).
#20 Posted: 8/5/2017 12:41:32 AM
Ruben Sidranski

Ruben Sidranski

17 likes in 463 posts.

Group: User

mkrask,
maybe it's my misunderstanding and lack of familiarity with the Maple plugin

Jean,
Finally, an example that shows WHY we need units the items in the array are in kW and need to be in BTU/hr for the empiracle equation to be valid. If I split the array like you did, it does what i need ... wondering why I cannot do it with vectorize Solve Block 2a.smz (4 KiB) downloaded 38 time(s).
  • New Posts New Posts
  • No New Posts No New Posts