solve & roots with units

solve & roots with units - Messages

#1 Posted: 2023/6/3 08:48:30
Valery Ochkov

Valery Ochkov

64 likes in 670 posts.

Group: User

I can it solve in Mathcad - see please mcdx and pdf files
I try solve it (heat transfer) in SMath - see sm file
Help please
File not found.File not found.File not found.
Kerimov-Cvetkov-Velichko-1-7.png
#2 Posted: 2023/6/3 15:37:44
overlord

overlord

580 likes in 1,348 posts.

Group: Moderator

Wrote

I can it solve in Mathcad - see please mcdx and pdf files
I try solve it (heat transfer) in SMath - see sm file
Help please


SMath symbolic solving of integration is not capable as Mathcad.
Thus integrals can't be used with numeric solvers.
However, you can use a little trick to solve them.
Using maple or maxima plugins is the key.
Here is your problem solved with various ways, roots() and Findroot().

Regards

solvewithunits.sm (23.4 KiB) downloaded 795 time(s).
3 users liked this post
Valery Ochkov 2023/6/3 15:53:00, sergio 2023/6/4 13:03:00, NDTM Amarasekera 2023/6/4 03:15:00
#3 Posted: 2023/6/4 02:18:41
Alvaro

Alvaro

1,078 likes in 1,749 posts.

Group: User

Hi. The attached have an al_nleqsol wrapper for a Mathcad block with units, that can accept SMath's integrals too.

Click to enlarge

Kerimov-Cvetkov-Velichko-1-7 - nSol.sm (17.97 KiB) downloaded 787 time(s).

Best regards.
Alvaro.
2 users liked this post
sergio 2023/6/4 13:03:00, Valery Ochkov 2023/6/4 02:22:00
#4 Posted: 2023/6/4 10:00:56
overlord

overlord

580 likes in 1,348 posts.

Group: Moderator

Actually, you don't need to find Tc2 first hand.
Findroot() can solve with all three equations at once.
All important is the order of unknowns.

Regards

2023-06-04_15-58.png
2 users liked this post
Valery Ochkov 2023/6/4 12:03:00, sergio 2023/6/4 13:48:00
#5 Posted: 2023/6/5 00:08:21
overlord

overlord

580 likes in 1,348 posts.

Group: Moderator

Here is a paper-friendly version.
Also check how variable guess points are assigned.

Regards

2023-06-05_06-05.png
1 users liked this post
Valery Ochkov 2023/6/5 02:23:00
#6 Posted: 2023/6/5 02:24:02
Valery Ochkov

Valery Ochkov

64 likes in 670 posts.

Group: User

Wrote

Here is a paper-friendly version.
Also check how variable guess points are assigned.

Regards

2023-06-05_06-05.png


Thanks!
Send please the sm file.

#7 Posted: 2023/6/5 11:08:24
Jean Giraud

Jean Giraud

983 likes in 6,866 posts.

Group: User

Hello Valery,
In some application(s), you can associate FindRoot with epsilon.

Solve Diophantine Chemical Balance epsilon.sm (106.1 KiB) downloaded 787 time(s).
1 users liked this post
Valery Ochkov 2023/6/5 11:26:00
#8 Posted: 2023/6/5 17:39:04
overlord

overlord

580 likes in 1,348 posts.

Group: Moderator

Wrote

Thanks!
Send please the sm file.


solvewithunits.sm (14.13 KiB) downloaded 801 time(s).
#9 Posted: 2023/6/26 02:33:01
netsonicyxf

netsonicyxf

2 likes in 79 posts.

Group: User

Please help, why can't find root for this simple equation?
File not found.File not found.
root.sm (20.2 KiB) downloaded 757 time(s).
#10 Posted: 2023/6/26 05:54:06
Davide Carpi

Davide Carpi

1,435 likes in 2,913 posts.

Group: Moderator

Wrote

why can't find root for this simple equation?



Because there aren't unknows

2023-06-26 10_49_28-SMath Solver - [root.sm_].png
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
#11 Posted: 2023/6/26 08:45:01
Jean Giraud

Jean Giraud

983 likes in 6,866 posts.

Group: User

Wrote

Please help,
why can't find root for this simple equation ?


roots solves wrt scalar function(s)
In this example, both functions are scalar wrt 'x'.
Try ... roots[(sin(x)-Gamma(x),x,1]
'x not defined' ... Gamma(x) is not scalar wrt 'x'
Few more numerical functions are not scalar
erf(x), erf(x,n) ... erfc(x), erfc(x,n) ...

NLroots.PNG
#12 Posted: 2023/6/26 08:55:02
Jean Giraud

Jean Giraud

983 likes in 6,866 posts.

Group: User

Wrote

Try ... roots[(sin(x)-Gamma(x),x,1]


In that case and similar ones,
RootSecant algorithm will spits out 15 decimals, instantaneous.
Cheers ... Jean.

#13 Posted: 2023/6/26 10:06:59
netsonicyxf

netsonicyxf

2 likes in 79 posts.

Group: User

Wrote

Wrote

why can't find root for this simple equation?



Because there aren't unknows

2023-06-26 10_49_28-SMath Solver - [root.sm_].png



Thanks Davide.
Now I put "N" as unknown in the Eq. But it can't be solved either.
1.png

root.sm (19.41 KiB) downloaded 761 time(s).
#14 Posted: 2023/6/26 10:58:58
Davide Carpi

Davide Carpi

1,435 likes in 2,913 posts.

Group: Moderator

You have to play a little with the precision, especially if there are units involved; in this case, seems that 10^-3 mm is a stable target (10^-2 works too, but in some case it returns the guess value, meaning that is "close enough" for the algorithm - you can see in the plot that 0kN to 200kN is in +0.02/-0.02mm range)

2023-06-26 15_55_45-SMath Solver - [root_r1.sm].png

root_r1.sm (24.03 KiB) downloaded 789 time(s).
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
netsonicyxf 2023/6/26 12:18:00
  • New Posts
  • No New Posts