solve & roots with units

solve & roots with units - Сообщения

#1 Опубликовано: 03.06.2023 08:48:30
Valery Ochkov

Valery Ochkov

55 сообщений из 633 понравились пользователям.

Группа: 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
Файл не найден.Файл не найден.Файл не найден.
Kerimov-Cvetkov-Velichko-1-7.png
#2 Опубликовано: 03.06.2023 15:37:44
overlord

overlord

548 сообщений из 1332 понравились пользователям.

Группа: 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 КиБ) скачан 44 раз(а).
3 пользователям понравился этот пост
Valery Ochkov 03.06.2023 15:53:00, sergio 04.06.2023 13:03:00, NDTM Amarasekera 04.06.2023 03:15:00
#3 Опубликовано: 04.06.2023 02:18:41
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 сообщений из 1674 понравились пользователям.

Группа: 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 КиБ) скачан 49 раз(а).

Best regards.
Alvaro.
2 пользователям понравился этот пост
sergio 04.06.2023 13:03:00, Valery Ochkov 04.06.2023 02:22:00
#4 Опубликовано: 04.06.2023 10:00:56
overlord

overlord

548 сообщений из 1332 понравились пользователям.

Группа: 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 пользователям понравился этот пост
Valery Ochkov 04.06.2023 12:03:00, sergio 04.06.2023 13:48:00
#5 Опубликовано: 05.06.2023 00:08:21
overlord

overlord

548 сообщений из 1332 понравились пользователям.

Группа: Moderator

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

Regards

2023-06-05_06-05.png
1 пользователям понравился этот пост
Valery Ochkov 05.06.2023 02:23:00
#6 Опубликовано: 05.06.2023 02:24:02
Valery Ochkov

Valery Ochkov

55 сообщений из 633 понравились пользователям.

Группа: 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 Опубликовано: 05.06.2023 11:08:24
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

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

Solve Diophantine Chemical Balance epsilon.sm (106 КиБ) скачан 34 раз(а).
1 пользователям понравился этот пост
Valery Ochkov 05.06.2023 11:26:00
#8 Опубликовано: 05.06.2023 17:39:04
overlord

overlord

548 сообщений из 1332 понравились пользователям.

Группа: Moderator

Wrote

Thanks!
Send please the sm file.


solvewithunits.sm (14 КиБ) скачан 38 раз(а).
#9 Опубликовано: 26.06.2023 02:33:01
netsonicyxf

netsonicyxf

2 сообщений из 79 понравились пользователям.

Группа: User

Please help, why can't find root for this simple equation?
Файл не найден.Файл не найден.
root.sm (20 КиБ) скачан 26 раз(а).
#10 Опубликовано: 26.06.2023 05:54:06
Davide Carpi

Davide Carpi

1415 сообщений из 2872 понравились пользователям.

Группа: 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 Опубликовано: 26.06.2023 08:45:01
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: 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 Опубликовано: 26.06.2023 08:55:02
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: 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 Опубликовано: 26.06.2023 10:06:59
netsonicyxf

netsonicyxf

2 сообщений из 79 понравились пользователям.

Группа: 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 КиБ) скачан 21 раз(а).
#14 Опубликовано: 26.06.2023 10:58:58
Davide Carpi

Davide Carpi

1415 сообщений из 2872 понравились пользователям.

Группа: 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 КиБ) скачан 27 раз(а).
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
1 пользователям понравился этот пост
netsonicyxf 26.06.2023 12:18:00
  • Новые сообщения Новые сообщения
  • Нет новых сообщений Нет новых сообщений