## 115 Reputation

9 years, 114 days

## How can I extract an excellent approxima...

Maple

Dear All,

I have a polynomial in terms of time. I know that based on the nature of the problem, the real function that governs this problem is the sum of exponential functions with a negative power, for example, in the form of alpha[0]+add(alpha[i]*exp(-beta[i]*t), i=1.. 5).

Can you help me if there is a method that can be used to obtain these exponential functions using following polynomial?

The polynomial function is as follow:

f:=0.020399949322360296902872908942 + 0.0261353198432118595103693714851*t^3 + 0.0240968505875842806805439681431*t^4 + 0.0148456155621193706595799212802*t^5 + 0.0239969764160351203722354728376*t^2 + 0.0204278458408370651586217048716*t - 0.00450853634927256388740864146173*t^6 - 0.0355389767483113696513996149731*t^7 - 0.0766669789661906882315038416910*t^8 - 0.120843030849135239578151569663*t^9 - 0.153280689906711146639066606024*t^10 - 0.150288711858517536713273977277*t^11 - 0.0808171080937786380164380347445*t^12 + 0.0872390654213369913348407061899*t^13 + 0.373992140377042586618283139889*t^14 + 0.766807288928470485618700282187*t^15 + 1.19339994493571167326973251788*t^16 + 1.49476369302534328383069681700*t^17 + 1.41015598591182237492637420929*t^18 + 0.593451797299651247527539427688*t^19 - 1.31434443870999971750661332301*t^20;

Best wishes

## How can I use hardware floating-point co...

Maple

Dear All,

I want all calculations in Maple code to be done with double precision. Using Digits:=18 and evalf greatly affects the calculation speed. But using the evalhf performs calculations much faster. The problem is that the evalhf cannot work with expressions that contain symbolic terms. For example, the evalhf(1/3*x+1/5*exp(-x^2)) generates an error. How can I perform calculations on expression with double-precision floating point without significantly slowing down the calculation speed? Is there a method that can be used at the beginning of the code to use double precision (preferably hardware) in all computation steps without slowing down the computations, unlike software floating-point computations?
Best wishes

Digits:=18;
evalf(1/3.);
evalhf(1/3.);

str:=time[real]();
for i from 1 to 100000 do
evalf(sqrt(3)*i):
end do:
time[real]()-str;

str:=time[real]();
for i from 1 to 100000 do
evalhf(sqrt(3)*i):
end do:
time[real]()-str;

## How can I convert exponential functions ...

Maple

Dear All,

I have the following algebraic function. In vibrations, it is common to write the functions exp(alpha[i]*I*t) as separate terms of cos(alpha[i]*t) and sin(alpha[i]*t). While exp(-beta[i]*t) remains without converting to sinh and cosh. How can I find the solution given for T(t) as the sum of terms C1[i]*sin(alpha[i]*t)*exp(-beta[i]*t) and C2[i]*cos(alpha[i]*t)*exp(-beta[i]*t) with non-complex (real) coefficients C1[i] and C2[i].

Can anyone help me to achieve my goal in the following expression?

T(t):=(1.450761945*10^(-11) + (3.836655196*10^(-14))*I)*exp((-0.5000000000 + 222.6866468*I)*t) + (-3.770333746*10^(-9) + (2.179000257*10^(-6))*I)*exp((-0.5000000000 - 924.5904413*I)*t) + (-2.584086158*10^(-12) + (4.273321932*10^(-13))*I)*exp((-0.5000000000 + 326.7549627*I)*t) + (1.986287340*10^(-9) + (1.330623218*10^(-11))*I)*exp((-0.5000000000 - 74.63720909*I)*t) + (-5.980910367*10^(-12) + (5.816480027*10^(-11))*I)*exp((-0.5000000000 - 453.7574402*I)*t) + (1.450761945*10^(-11) - (3.836655196*10^(-14))*I)*exp((-0.5000000000 - 222.6866468*I)*t) + (8.923968224*10^(-10) - (8.844466162*10^(-9))*I)*exp((-0.5000000000 + 637.9999953*I)*t) - (1.217986141*10^(-10) + (4.431771836*10^(-13))*I)*exp((-0.5000000000 - 138.7904660*I)*t) + (-1.217986141*10^(-10) + (4.431771836*10^(-13))*I)*exp((-0.5000000000 + 138.7904660*I)*t) + (0.0002537882980 + 0.00002277791755*I)*exp((-0.5000000000 - 5.570928456*I)*t) - (3.770333746*10^(-9) + (2.179000257*10^(-6))*I)*exp((-0.5000000000 + 924.5904413*I)*t) + (-0.0001618723219 + 0.01288449595*I)*exp((-0.5000000000 - 1638.001654*I)*t) + (8.923968224*10^(-10) + (8.844466162*10^(-9))*I)*exp((-0.5000000000 - 637.9999953*I)*t) + (-1.153529195*10^(-7) + (1.908444485*10^(-9))*I)*exp((-0.5000000000 + 30.22171212*I)*t) - (1.153529195*10^(-7) + (1.908444485*10^(-9))*I)*exp((-0.5000000000 - 30.22171212*I)*t) - (0.0001618723219 + 0.01288449595*I)*exp((-0.5000000000 + 1638.001654*I)*t) + (1.986287340*10^(-9) - (1.330623218*10^(-11))*I)*exp((-0.5000000000 + 74.63720909*I)*t) + (0.0002537882980 - 0.00002277791755*I)*exp((-0.5000000000 + 5.570928456*I)*t) - (5.980910367*10^(-12) + (5.816480027*10^(-11))*I)*exp((-0.5000000000 + 453.7574402*I)*t) - (2.584086158*10^(-12) + (4.273321932*10^(-13))*I)*exp((-0.5000000000 - 326.7549627*I)*t);

## How can I extract the specific terms in ...

Maple

Dear All,

I have the given expression f as follows. I want to extract terms with sin, sinh, and exp in this expression. For example, for sin, cos, sinh and exp, I want the outputs to be -4*sin(x), -5*cos(x^3)*sin(y^2), 5*sinh(x^2) and 2*exp(y^2), respectively. Also, if I want to get terms having sin(x) and sin(y^2), I like to receive -4*sin(x) and -5*cos(x^3)*sin(y^2), respectively. And also, for terms including sinh(x^3), I want to get Void output. The above goals cannot be achieved with the following commands. Can anyone help me?

f := -4*sin(x) + 2*exp(y^2) + 5 - 5*cos(x^3)*sin(y^2) + 5*sinh(x^2);

indets(f, trig);

selectfun(f, exp);

selectremove(has, f, trig);

## How can I solve a set of linear differen...

Maple 2024

Deal All,

I have a system of linear differential equations with unknown functions T[1](t) to T[n](t). In the attached example, I considered the value of n equal to 10, but depending on the problem, the value of n may be higher. Maple is not able to solve this problem analytically with the ‘dsolve’.

Does anyone have an idea to analytically solve for such a set of linear differential equations?

Best wishes

Set_of_Linear_DEs.mw

 1 2 3 4 5 Page 1 of 5
﻿