solve & roots with units

solve & roots with units - Messaggi

#1 Pubblicato: 03/06/2023 08:48:30
Valery Ochkov

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: 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 non trovato.File non trovato.File non trovato.
Kerimov-Cvetkov-Velichko-1-7.png
#2 Pubblicato: 03/06/2023 15:37:44
overlord

overlord

600 apprezzamenti in 1.381 post.

Gruppo: 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) scaricato 1426 volta(e).
3 utenti hanno apprezzato questo post
Valery Ochkov 03/06/2023 15:53:00, sergio 04/06/2023 13:03:00, NDTM Amarasekera 04/06/2023 03:15:00
#3 Pubblicato: 04/06/2023 02:18:41
Alvaro

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: 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) scaricato 1458 volta(e).

Best regards.
Alvaro.
2 utenti hanno apprezzato questo post
sergio 04/06/2023 13:03:00, Valery Ochkov 04/06/2023 02:22:00
#4 Pubblicato: 04/06/2023 10:00:56
overlord

overlord

600 apprezzamenti in 1.381 post.

Gruppo: 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 utenti hanno apprezzato questo post
Valery Ochkov 04/06/2023 12:03:00, sergio 04/06/2023 13:48:00
#5 Pubblicato: 05/06/2023 00:08:21
overlord

overlord

600 apprezzamenti in 1.381 post.

Gruppo: Moderator

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

Regards

2023-06-05_06-05.png
1 utenti hanno apprezzato questo post
Valery Ochkov 05/06/2023 02:23:00
#6 Pubblicato: 05/06/2023 02:24:02
Valery Ochkov

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: 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 Pubblicato: 05/06/2023 11:08:24
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

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

Solve Diophantine Chemical Balance epsilon.sm (106,1 KiB) scaricato 1489 volta(e).
1 utenti hanno apprezzato questo post
Valery Ochkov 05/06/2023 11:26:00
#8 Pubblicato: 05/06/2023 17:39:04
overlord

overlord

600 apprezzamenti in 1.381 post.

Gruppo: Moderator

Wrote

Thanks!
Send please the sm file.


solvewithunits.sm (14,13 KiB) scaricato 1416 volta(e).
#9 Pubblicato: 26/06/2023 02:33:01
netsonicyxf

netsonicyxf

2 apprezzamenti in 79 post.

Gruppo: User

Please help, why can't find root for this simple equation?
File non trovato.File non trovato.
root.sm (20,2 KiB) scaricato 1394 volta(e).
#10 Pubblicato: 26/06/2023 05:54:06
Davide Carpi

Davide Carpi

1.443 apprezzamenti in 2.936 post.

Gruppo: 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 Pubblicato: 26/06/2023 08:45:01
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: 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 Pubblicato: 26/06/2023 08:55:02
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: 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 Pubblicato: 26/06/2023 10:06:59
netsonicyxf

netsonicyxf

2 apprezzamenti in 79 post.

Gruppo: 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) scaricato 1361 volta(e).
#14 Pubblicato: 26/06/2023 10:58:58
Davide Carpi

Davide Carpi

1.443 apprezzamenti in 2.936 post.

Gruppo: 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) scaricato 1430 volta(e).
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
1 utenti hanno apprezzato questo post
netsonicyxf 26/06/2023 12:18:00
  • Nuovi post
  • Nessun nuovo post