gaurav_rs

65 Reputation

4 Badges

7 years, 239 days

MaplePrimes Activity


These are replies submitted by gaurav_rs

@Carl Love thanks for the comment. My digits environment is set to 50 and I have one symbol "x" is involved in some calculations . Sorry for the confusion , 2 hours is the total time (of generating the matrix and solving them).

If there is a symbol then for no of terms = 60 is takes 3-4 hours depending on the the coefficient in the matrix A.

It is good to know that LinearSOlve cant be parallelized.

Thanks .

@Carl Love n is 400 and coefficients are floats. In my pc it takes 2  hours and so. I was thinking the posibility of parallel computation for AX=B, using solve and hybrid command both.

I am looking for a simple script file

to schedule my script file on a server having 24 cores(Linux machine)

for parallel computation of

sol:=LinearSolve(A, B, method = 'hybrid/solve');

regards,

@vv Thanks for the solution. I was expecting real factor though. Thanks!

@Preben Alsholm tahnks for the reply. I am sorry for the mess. By mistake, I pasted the wrong poluynomial.I have modified my question, now.

 

Thanks &regards,

grv

@vv I have used the command but maple is unable give the solution. It is taking a lot of time! any guess why is so?

@tomleslie I rechecked ,A11:= ImportVector("F:/Maplework/abc.txt", source=delimited, delimiter= " "), imports whole expression (123*y+567; 834*y+546 ...) in maple 18,  where abc.txt is a vector contains string of type( 123*y^n+567*y^(n-1)+..) . But for the the same in Maple 2017 i have to specify , A11:= ImportVector("F:/Maplework/abc.txt", source=delimited, delimiter= " ", datatype=string) and then I have to use A22:=  map(parse,A11) to convert it back in some usable form.

I dont know how this feature in Maple 2017.2 makes it better than Maple 18. It just took some time to figure it out.

Export command is same  for both Maple 18 and 2017.2.

 

regards,

 

@tomleslie Thankyou so much. It works!

 

Thnaks and regards,

 

@tomleslie thanks it works for floating point data. But is there a way to tell the system that data type = anything.

see if text file contains matrix with algebraic entries like 1234*x+5467 , datatype=float wont work . By default it should take datatype= anything but during import if i use

MatA:=ImportMatrix("/home/15_degree_3izto1_fixed/Mat.txt", source = delimited) ; 

Imported matrix contains float[8] only. i.e it only imports "1234" and NOT the value  "1234*x+5467" that it must be importing.

 

thanks and regards,

@Kitonum thanks a lot. It solves my problem. Thanks again for your time and help. :-)

 

 

@taro Thanks for the reply I have appended the file.

@Kitonum Thanks a lot for helping me out. There is fraction involved so direct use of solve explicit doesnot work for me. Can you suggest some modirfication in my work sheet? I tried separately for numer and denom , it worked but the roots were very combursome.

I am attaching the work sheet here. Thanks.
 

Polynomial Example prob

 

case 1

 

 

``

restart;

p:=(delta^2-1)*(2*(1+nu)*Omega*I_in/(F*a*s^2))/((1+nu)*((delta^2-1))/k+3*p1/q1*(delta^2+3*q1/(q1+6*p1*k_b)));

2*(delta^2-1)*(1+nu)*Omega*I_in/(F*a*s^2*((1+nu)*(delta^2-1)/k+3*p1*(delta^2+3*q1/(6*k_b*p1+q1))/q1))

(1.1.1)

p1:=1+(tau1+tau2)*s+tau1*tau2*s^2;

1+(tau1+tau2)*s+tau1*tau2*s^2

(1.1.2)

q1:=2*(tau1*g1+tau2*g2)*s+2*tau1*tau2*(g1+g2)*s^2;

2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2

(1.1.3)

p;

2*(delta^2-1)*(1+nu)*Omega*I_in/(F*a*s^2*((1+nu)*(delta^2-1)/k+3*(1+(tau1+tau2)*s+tau1*tau2*s^2)*(delta^2+3*(2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2)/(6*k_b*(1+(tau1+tau2)*s+tau1*tau2*s^2)+2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2))/(2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2)))

(1.1.4)

L:=[solve(p, s, explicit)]:

mul(s-p, p=L);   # Factorization in symbolic form

#eval(%, [a=1, b=2, c=-1, d=-4, e=-3]);  # Replace symbols with numbers (in symbolic form)

#evalf(%);  # Numerical factorization

 

(s+(g1*tau1+g2*tau2)/(tau1*tau2*(g1+g2)))*(s-(1/2)*(-tau1*g1-tau2*g2-3*k_b*tau1-3*k_b*tau2+(g1^2*tau1^2+2*g1*g2*tau1*tau2+6*g1*k_b*tau1^2-6*g1*k_b*tau1*tau2+g2^2*tau2^2-6*g2*k_b*tau1*tau2+6*g2*k_b*tau2^2+9*k_b^2*tau1^2-18*k_b^2*tau1*tau2+9*k_b^2*tau2^2)^(1/2))/(tau1*tau2*(g1+g2+3*k_b)))*(s+(1/2)*(tau1*g1+tau2*g2+3*k_b*tau1+3*k_b*tau2+(g1^2*tau1^2+2*g1*g2*tau1*tau2+6*g1*k_b*tau1^2-6*g1*k_b*tau1*tau2+g2^2*tau2^2-6*g2*k_b*tau1*tau2+6*g2*k_b*tau2^2+9*k_b^2*tau1^2-18*k_b^2*tau1*tau2+9*k_b^2*tau2^2)^(1/2))/(tau1*tau2*(g1+g2+3*k_b)))

(1.1.5)

factor(p);

4*(delta-1)*(delta+1)*(1+nu)*Omega*I_in*k*(g1*s*tau1*tau2+g2*s*tau1*tau2+g1*tau1+g2*tau2)*(g1*s^2*tau1*tau2+g2*s^2*tau1*tau2+3*k_b*s^2*tau1*tau2+g1*s*tau1+g2*s*tau2+3*k_b*s*tau1+3*k_b*s*tau2+3*k_b)/(s*F*a*(2*delta^2*g1^2*nu*s^4*tau1^2*tau2^2+4*delta^2*g1*g2*nu*s^4*tau1^2*tau2^2+6*delta^2*g1*k_b*nu*s^4*tau1^2*tau2^2+2*delta^2*g2^2*nu*s^4*tau1^2*tau2^2+6*delta^2*g2*k_b*nu*s^4*tau1^2*tau2^2+2*delta^2*g1^2*s^4*tau1^2*tau2^2+4*delta^2*g1*g2*s^4*tau1^2*tau2^2+3*delta^2*g1*k*s^4*tau1^2*tau2^2+6*delta^2*g1*k_b*s^4*tau1^2*tau2^2+2*delta^2*g2^2*s^4*tau1^2*tau2^2+3*delta^2*g2*k*s^4*tau1^2*tau2^2+6*delta^2*g2*k_b*s^4*tau1^2*tau2^2+9*delta^2*k*k_b*s^4*tau1^2*tau2^2+4*delta^2*g1^2*nu*s^3*tau1^2*tau2+4*delta^2*g1*g2*nu*s^3*tau1^2*tau2+4*delta^2*g1*g2*nu*s^3*tau1*tau2^2+12*delta^2*g1*k_b*nu*s^3*tau1^2*tau2+6*delta^2*g1*k_b*nu*s^3*tau1*tau2^2+4*delta^2*g2^2*nu*s^3*tau1*tau2^2+6*delta^2*g2*k_b*nu*s^3*tau1^2*tau2+12*delta^2*g2*k_b*nu*s^3*tau1*tau2^2-2*g1^2*nu*s^4*tau1^2*tau2^2-4*g1*g2*nu*s^4*tau1^2*tau2^2-6*g1*k_b*nu*s^4*tau1^2*tau2^2-2*g2^2*nu*s^4*tau1^2*tau2^2-6*g2*k_b*nu*s^4*tau1^2*tau2^2+4*delta^2*g1^2*s^3*tau1^2*tau2+4*delta^2*g1*g2*s^3*tau1^2*tau2+4*delta^2*g1*g2*s^3*tau1*tau2^2+6*delta^2*g1*k*s^3*tau1^2*tau2+3*delta^2*g1*k*s^3*tau1*tau2^2+12*delta^2*g1*k_b*s^3*tau1^2*tau2+6*delta^2*g1*k_b*s^3*tau1*tau2^2+4*delta^2*g2^2*s^3*tau1*tau2^2+3*delta^2*g2*k*s^3*tau1^2*tau2+6*delta^2*g2*k*s^3*tau1*tau2^2+6*delta^2*g2*k_b*s^3*tau1^2*tau2+12*delta^2*g2*k_b*s^3*tau1*tau2^2+18*delta^2*k*k_b*s^3*tau1^2*tau2+18*delta^2*k*k_b*s^3*tau1*tau2^2-2*g1^2*s^4*tau1^2*tau2^2-4*g1*g2*s^4*tau1^2*tau2^2+9*g1*k*s^4*tau1^2*tau2^2-6*g1*k_b*s^4*tau1^2*tau2^2-2*g2^2*s^4*tau1^2*tau2^2+9*g2*k*s^4*tau1^2*tau2^2-6*g2*k_b*s^4*tau1^2*tau2^2+2*delta^2*g1^2*nu*s^2*tau1^2+4*delta^2*g1*g2*nu*s^2*tau1*tau2+6*delta^2*g1*k_b*nu*s^2*tau1^2+12*delta^2*g1*k_b*nu*s^2*tau1*tau2+2*delta^2*g2^2*nu*s^2*tau2^2+12*delta^2*g2*k_b*nu*s^2*tau1*tau2+6*delta^2*g2*k_b*nu*s^2*tau2^2-4*g1^2*nu*s^3*tau1^2*tau2-4*g1*g2*nu*s^3*tau1^2*tau2-4*g1*g2*nu*s^3*tau1*tau2^2-12*g1*k_b*nu*s^3*tau1^2*tau2-6*g1*k_b*nu*s^3*tau1*tau2^2-4*g2^2*nu*s^3*tau1*tau2^2-6*g2*k_b*nu*s^3*tau1^2*tau2-12*g2*k_b*nu*s^3*tau1*tau2^2+2*delta^2*g1^2*s^2*tau1^2+4*delta^2*g1*g2*s^2*tau1*tau2+3*delta^2*g1*k*s^2*tau1^2+6*delta^2*g1*k*s^2*tau1*tau2+6*delta^2*g1*k_b*s^2*tau1^2+12*delta^2*g1*k_b*s^2*tau1*tau2+2*delta^2*g2^2*s^2*tau2^2+6*delta^2*g2*k*s^2*tau1*tau2+3*delta^2*g2*k*s^2*tau2^2+12*delta^2*g2*k_b*s^2*tau1*tau2+6*delta^2*g2*k_b*s^2*tau2^2+9*delta^2*k*k_b*s^2*tau1^2+36*delta^2*k*k_b*s^2*tau1*tau2+9*delta^2*k*k_b*s^2*tau2^2-4*g1^2*s^3*tau1^2*tau2-4*g1*g2*s^3*tau1^2*tau2-4*g1*g2*s^3*tau1*tau2^2+18*g1*k*s^3*tau1^2*tau2+9*g1*k*s^3*tau1*tau2^2-12*g1*k_b*s^3*tau1^2*tau2-6*g1*k_b*s^3*tau1*tau2^2-4*g2^2*s^3*tau1*tau2^2+9*g2*k*s^3*tau1^2*tau2+18*g2*k*s^3*tau1*tau2^2-6*g2*k_b*s^3*tau1^2*tau2-12*g2*k_b*s^3*tau1*tau2^2+6*delta^2*g1*k_b*nu*s*tau1+6*delta^2*g2*k_b*nu*s*tau2-2*g1^2*nu*s^2*tau1^2-4*g1*g2*nu*s^2*tau1*tau2-6*g1*k_b*nu*s^2*tau1^2-12*g1*k_b*nu*s^2*tau1*tau2-2*g2^2*nu*s^2*tau2^2-12*g2*k_b*nu*s^2*tau1*tau2-6*g2*k_b*nu*s^2*tau2^2+3*delta^2*g1*k*s*tau1+6*delta^2*g1*k_b*s*tau1+3*delta^2*g2*k*s*tau2+6*delta^2*g2*k_b*s*tau2+18*delta^2*k*k_b*s*tau1+18*delta^2*k*k_b*s*tau2-2*g1^2*s^2*tau1^2-4*g1*g2*s^2*tau1*tau2+9*g1*k*s^2*tau1^2+18*g1*k*s^2*tau1*tau2-6*g1*k_b*s^2*tau1^2-12*g1*k_b*s^2*tau1*tau2-2*g2^2*s^2*tau2^2+18*g2*k*s^2*tau1*tau2+9*g2*k*s^2*tau2^2-12*g2*k_b*s^2*tau1*tau2-6*g2*k_b*s^2*tau2^2-6*g1*k_b*nu*s*tau1-6*g2*k_b*nu*s*tau2+9*delta^2*k*k_b+9*g1*k*s*tau1-6*g1*k_b*s*tau1+9*g2*k*s*tau2-6*g2*k_b*s*tau2))

(1.1.6)

 

 

 

 

simplify(p,size);

2*(delta^2-1)*(1+nu)*Omega*I_in/(F*a*s^2*((1+nu)*(delta^2-1)/k+3*(1+(tau1+tau2)*s+tau1*tau2*s^2)*(delta^2+3*(2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2)/(6*k_b*(1+(tau1+tau2)*s+tau1*tau2*s^2)+2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2))/(2*(g1*tau1+g2*tau2)*s+2*tau1*tau2*(g1+g2)*s^2)))

(1.1.7)

 

#assume(k>0,F>0,tau1>0,tau2>0,k_b>0,g1>0,g2>0,Omega>0,I_in>0,g1>g2,tau1>tau2,nu<1,nu>0,k>k_b,delta>1);

 

 

 

 

#input

g1:=0.9167*10^(9);

916700000.0

(1.1.8)

``

g2:=0.0768*10^(9);

76800000.00

(1.1.9)

tau1:=1.047*10^(5);

104700.000

(1.1.10)

tau2:=678.43;

678.43

(1.1.11)

Di:=1*10^(-16);

1/10000000000000000

(1.1.12)

Omega:=3/96500;

3/96500

(1.1.13)

k_b:=2*(g1+g2)*(1+.33)/3/(1-2*.33);

2590892157.

(1.1.14)

delta:=1.1;

1.1

(1.1.15)

I_in:=0.176;a:=1*10^(-6);F:=96500;

.176

 

1/1000000

 

96500

(1.1.16)

nu:=0.22;

.22

(1.1.17)

k:=47.619*10^9;

0.4761900000e11

(1.1.18)

p;

0.2905282826e-4/(s^2*(0.5380205380e-11+3*(1+105378.430*s+71031621.00*s^2)*(1.21+3*(0.1920611868e15*s+0.1411398309e18*s^2)/(0.1554535294e11+0.1830206074e16*s+0.1245351449e19*s^2))/(0.1920611868e15*s+0.1411398309e18*s^2)))

(1.1.19)

simplify(p,size);## desired output

12386.14110*(s+0.1360786573e-2)*(s+0.1461086742e-2)*(s+0.8543437688e-5)/(s*(s+HFloat(0.0014738955790471226))*(s+HFloat(0.001438868727786942))*(s+HFloat(9.543679813919816e-6))*(s+HFloat(6.762535285302987e-6)))

(1.1.20)

``

 

 

 

NULL

 

 

 


 

Download polynomial_prime.mw

@ThU I think you are right. For the cross check I have included the worksheet. Is it possible that the factor causes some error? How to cross check the solution (inv laplace) in such a case. Can you verify the poles ? Thanks,


 

################################################Inverting the result############################################

restart

Digits := 30;

30

(1)

``

A11 := -2.77777777777777777777777777778*10^(-29)*(9.28328070200010721495878047583*10^115*s^21+1.88205134700002323105272125304*10^112*s^20+3.91011721800004757051577536443*10^108*s^19+7.37738589300009183336634204689*10^104*s^18+1.23421857700001666775270370158*10^101*s^17+2.17979818780000995534264567186*10^461*s^150+4.45581406430002019209018913007*10^460*s^149+4.52276392550002033970405022044*10^459*s^148+7.49532505000022439494280487829*10^29+3.03937939810001356642550101498*10^458*s^147+1.52116811486000674025578681202*10^457*s^146+6.04628843990002660639763049503*10^455*s^145+1.98861721540000869166374199996*10^454*s^144+2.15309248170000861165557599155*10^411*s^122+3.62467906190001452212704247088*10^413*s^123+5.83511145030002342265060275116*10^415*s^124+8.97429079890003608578695079299*10^417*s^125+1.31553106925000530292559392169*10^420*s^126+1.83682535080000742174509610723*10^422*s^127+2.43748537550000987597492079021*10^424*s^128+3.06930822420001247217316521474*10^426*s^129+3.65976264080001491833551153549*10^428*s^130+4.12238295030001686376283829074*10^430*s^131+4.37755990140001797143262195171*10^432*s^132+4.36846009910001800633150896055*10^434*s^133+4.08585361910001691080126374422*10^436*s^134+3.56881489420001483605284193611*10^438*s^135+4.55318344420002073496747591314*10^312*s^83+2.72421757870001228130820857468*10^315*s^84+1.58024096250000708089641797370*10^318*s^85+8.93668093200003973695625187472*10^320*s^86+4.91317943050002170299277211070*10^323*s^87+2.62833174310001153492217322931*10^326*s^88+1.37016854370000596521884425070*10^329*s^89+6.93571055330003001211522080175*10^331*s^90+3.41273834680001468798550008504*10^334*s^91+1.63515193670000699125382287171*10^337*s^92+7.61591553560003235948852876455*10^339*s^93+3.44568355100001456206860453459*10^342*s^94+1.51461409410000636993220266389*10^345*s^95+6.47450043680002707986410347817*10^347*s^96+2.68724045360001118566024827147*10^350*s^97+1.08304244430000448825418791895*10^353*s^98+4.23599295330001748978758982077*10^355*s^99+1.61039436940000661708284840187*10^358*s^100+5.93580044990002429964836316009*10^360*s^101+2.12121960120000865871837300386*10^363*s^102+7.35986459290002992872945618326*10^365*s^103+2.47289139030001003128096144474*10^368*s^104+8.06085958040003259128397622795*10^370*s^105+2.54291906400001026024227195584*10^373*s^106+7.77739398330003128592291167182*10^375*s^107+2.30035521880000923610991587148*10^378*s^108+1.28354910500002512537605430497*10^35*s+2.16705950700002835629958642934*10^97*s^16+3.47822447500004506473903241377*10^93*s^15+4.67324092800006664429945402921*10^89*s^14+6.43589568000009130971303829866*10^85*s^13+7.57161964000011531888982670506*10^81*s^12+9.03002596000013346236450885806*10^77*s^11+9.54872022000014056852523327277*10^73*s^10+8.67745291000013363072006755635*10^69*s^9+6.92583409000011351951208550834*10^65*s^8+4.92739847000008511521837873418*10^61*s^7+3.51543375200005544957594509933*10^57*s^6+1.88278633900003074748887167538*10^53*s^5+7.55893919000014110021983105376*10^48*s^4+2.59767743600005144509222050811*10^44*s^3+6.88991849000013971717890352067*10^39*s^2+7.42533258700003408000299440322*10^309*s^82+1.17959287260000545336176808898*10^307*s^81+1.82212280300000849618663659604*10^304*s^80+2.73419886090001288846709325794*10^301*s^79+3.99972430820001903769089781830*10^298*s^78+5.70593556600002738276148632547*10^295*s^77+7.89931881200003835308432906323*10^292*s^76+1.06686955250000523105737212606*10^290*s^75+1.40889535250000694779222150066*10^287*s^74+1.80090539400000898608018179344*10^284*s^73+2.23658267220001131758646052921*10^281*s^72+2.71995483640001387993070391206*10^278*s^71+3.20928995460001657511405560610*10^275*s^70+3.69305884950001927272848975992*10^272*s^69+4.12495236460002181849976462563*10^269*s^68+4.49080132100002404777684680127*10^266*s^67+4.75373733300002580253988030717*10^263*s^66+4.95579513400002694959563915289*10^260*s^65+4.94278678700002739710778596481*10^257*s^64+4.83513520300002710664514775213*10^254*s^63+4.57438003200002609857208223865*10^251*s^62+4.26049678100002444974272353015*10^248*s^61+3.80244848900002228385097637212*10^245*s^60+3.32500520510001975611244362139*10^242*s^59+2.82977689630001703490893192329*10^239*s^58+2.33809497470001428340523835048*10^236*s^57+1.87592579740001164389842102630*10^233*s^56+1.46262735740000922688753161405*10^230*s^55+1.11044640610000710578114135945*10^227*s^54+8.18497419400005317062195156771*10^223*s^53+5.92500366300003864845736471578*10^220*s^52+4.10578860000002728246973954547*10^217*s^51+2.71032376800001869863761245154*10^214*s^50+1.78606627650001243905641095174*10^211*s^49+1.14059961900000802943378048675*10^208*s^48+6.97432114200005027629200260824*10^204*s^47+4.18130053600003052634403487099*10^201*s^46+2.43008078000001796648381585828*10^198*s^45+1.34718223730001024618583331216*10^195*s^44+7.34545282300005659738172472762*10^191*s^43+3.92219836300003026766904525089*10^188*s^42+1.94597758530001566435172077803*10^185*s^41+9.72173990600007841280215460624*10^181*s^40+4.60813624500003794723013721289*10^178*s^39+2.11127966000001774416130583591*10^175*s^38+9.32835075900008012391412379634*10^171*s^37+3.93827278100003491670636482573*10^168*s^36+1.62704233150001467522894935245*10^165*s^35+6.45943728700005944490085111741*10^161*s^34+2.51035787000002318996215775019*10^158*s^33+9.03103651500008705566961396355*10^154*s^32+3.26940643400003142229304518949*10^151*s^31+1.09759258840001089503997875561*10^148*s^30+3.55353083100003625333184772243*10^144*s^29+1.14777607920001156490353790688*10^141*s^28+3.33389114700003532836692736802*10^137*s^27+9.59580790600010322074170849868*10^133*s^26+2.64027489800002880741905453280*10^130*s^25+6.58295992600007668678998229860*10^126*s^24+1.64482553400001944229508494451*10^123*s^23+4.06500951300004686596373775530*10^119*s^22+6.58420753790002638640829172275*10^380*s^109+1.82194789150000729143183992837*10^383*s^110+4.87633159620001947888919393678*10^385*s^111+1.25957308120000502803295451832*10^388*s^112+3.14193814020001253324728898256*10^390*s^113+7.56540445970003015045176593336*10^392*s^114+1.75575288680000699527569771705*10^395*s^115+3.92686606920001564207115925321*10^397*s^116+8.45361411880003368524426334162*10^399*s^117+1.75170732900000698066907320027*10^402*s^118+3.48828727280001390908365014511*10^404*s^119+6.67120289230002662259959511905*10^406*s^120+1.22414851162000489024803226308*10^409*s^121+2.90025936640001211188816111900*10^440*s^136+2.18244894330000916080452775972*10^442*s^137+1.51409576770000638682993148519*10^444*s^138+9.62248985820004080491234341778*10^445*s^139+5.56382216630002372525779521630*10^447*s^140+2.90286276620001245084614961212*10^449*s^141+1.35296670027000583893011284349*10^451*s^142+5.56521322380002416799855195495*10^452*s^143)*(339.*s+.500000000000000000000000000000)/(s*(7.79437500000000000000000000000*10^317*s^95+1.74238750000000000000000000000*10^315*s^94+3.78725000000000000000000000000*10^312*s^93+8.00325000000000000000000000000*10^309*s^92+1.64450000000000000000000000000*10^307*s^91+3.28257500000000000000000000000*10^304*s^90+6.38550000000000000000000000000*10^301*s^89+5.87412500000000000000000000000*10^325*s^98+2.0262500*10^7*s+4.84500000000000000000000000000*10^69*s^16+6.80000000000000000000000000000*10^65*s^15+1.00250000000000000000000000000*10^62*s^14+1.37875000000000000000000000000*10^58*s^13+1.61250000000000000000000000000*10^54*s^12+1.79625000000000000000000000000*10^50*s^11+1.95250000000000000000000000000*10^46*s^10+1.67000000000000000000000000000*10^42*s^9+1.36000000000000000000000000000*10^38*s^8+1.01375000000000000000000000000*10^34*s^7+6.35000000000000000000000000000*10^29*s^6+3.4287500000000000000000000*10^25*s^5+1.538750000000000000000*10^21*s^4+5.2000000000000000*10^16*s^3+1.362500000000*10^12*s^2+3.13050000000000000000000000000*10^282*s^82+4.89387500000000000000000000000*10^279*s^81+7.45887500000000000000000000000*10^276*s^80+1.10663750000000000000000000000*10^274*s^79+1.60045000000000000000000000000*10^271*s^78+2.24975000000000000000000000000*10^268*s^77+3.07875000000000000000000000000*10^265*s^76+4.10512500000000000000000000000*10^262*s^75+5.32950000000000000000000000000*10^259*s^74+6.74212500000000000000000000000*10^256*s^73+8.30675000000000000000000000000*10^253*s^72+9.94175000000000000000000000000*10^250*s^71+1.15908750000000000000000000000*10^248*s^70+1.32112500000000000000000000000*10^245*s^69+1.45937500000000000000000000000*10^242*s^68+1.57050000000000000000000000000*10^239*s^67+1.65137500000000000000000000000*10^236*s^66+1.68062500000000000000000000000*10^233*s^65+1.67075000000000000000000000000*10^230*s^64+1.61387500000000000000000000000*10^227*s^63+1.51687500000000000000000000000*10^224*s^62+1.38875000000000000000000000000*10^221*s^61+1.23487500000000000000000000000*10^218*s^60+1.07025000000000000000000000000*10^215*s^59+8.98625000000000000000000000000*10^211*s^58+7.32750000000000000000000000000*10^208*s^57+5.85412500000000000000000000000*10^205*s^56+4.51575000000000000000000000000*10^202*s^55+3.40350000000000000000000000000*10^199*s^54+2.47362500000000000000000000000*10^196*s^53+1.75512500000000000000000000000*10^193*s^52+1.21175000000000000000000000000*10^190*s^51+8.03625000000000000000000000000*10^186*s^50+5.22000000000000000000000000000*10^183*s^49+3.29275000000000000000000000000*10^180*s^48+2.00825000000000000000000000000*10^177*s^47+1.18962500000000000000000000000*10^174*s^46+6.80500000000000000000000000000*10^170*s^45+3.78000000000000000000000000000*10^167*s^44+2.03025000000000000000000000000*10^164*s^43+1.05737500000000000000000000000*10^161*s^42+5.30000000000000000000000000000*10^157*s^41+2.60112500000000000000000000000*10^154*s^40+1.22712500000000000000000000000*10^151*s^39+5.53875000000000000000000000000*10^147*s^38+2.44387500000000000000000000000*10^144*s^37+1.03412500000000000000000000000*10^141*s^36+4.24250000000000000000000000000*10^137*s^35+1.68325000000000000000000000000*10^134*s^34+6.28625000000000000000000000000*10^130*s^33+2.28250000000000000000000000000*10^127*s^32+8.04000000000000000000000000000*10^123*s^31+2.71375000000000000000000000000*10^120*s^30+8.81500000000000000000000000000*10^116*s^29+2.70625000000000000000000000000*10^113*s^28+8.09500000000000000000000000000*10^109*s^27+2.28750000000000000000000000000*10^106*s^26+6.05000000000000000000000000000*10^102*s^25+1.59375000000000000000000000000*10^99*s^24+3.87875000000000000000000000000*10^95*s^23+9.27375000000000000000000000000*10^91*s^22+2.05750000000000000000000000000*10^88*s^21+4.07500000000000000000000000000*10^84*s^20+8.69625000000000000000000000000*10^80*s^19+1.61125000000000000000000000000*10^77*s^18+2.78625000000000000000000000000*10^73*s^17+6.09521250000000000000000000000*10^435*s^150+6.23350000000000000000000000000*10^434*s^149+4.22026250000000000000000000000*10^433*s^148+2.12785750000000000000000000000*10^432*s^147+8.52270000000000000000000000000*10^430*s^146+2.82440000000000000000000000000*10^429*s^145+7.96603750000000000000000000000*10^427*s^144+1.95178750000000000000000000000*10^426*s^143+4.21967500000000000000000000000*10^424*s^142+8.15070000000000000000000000000*10^422*s^141+1.42078125000000000000000000000*10^421*s^140+2.25338750000000000000000000000*10^419*s^139+3.27485000000000000000000000000*10^417*s^138+4.38636250000000000000000000000*10^415*s^137+5.44115000000000000000000000000*10^413*s^136+6.28072500000000000000000000000*10^411*s^135+6.77112500000000000000000000000*10^409*s^134+6.84078750000000000000000000000*10^407*s^133+6.49688750000000000000000000000*10^405*s^132+5.81567500000000000000000000000*10^403*s^131+4.91967500000000000000000000000*10^401*s^130+3.93996250000000000000000000000*10^399*s^129+2.99457500000000000000000000000*10^397*s^128+2.16400000000000000000000000000*10^395*s^127+1.48870000000000000000000000000*10^393*s^126+9.76702500000000000000000000000*10^390*s^125+6.12057500000000000000000000000*10^388*s^124+3.66785000000000000000000000000*10^386*s^123+2.10382500000000000000000000000*10^384*s^122+1.15718750000000000000000000000*10^382*s^121+6.10702500000000000000000000000*10^379*s^120+3.09496250000000000000000000000*10^377*s^119+1.50807500000000000000000000000*10^375*s^118+7.06730000000000000000000000000*10^372*s^117+3.19063750000000000000000000000*10^370*s^116+1.38778750000000000000000000000*10^368*s^115+5.82162500000000000000000000000*10^365*s^114+2.35592500000000000000000000000*10^363*s^113+9.20875000000000000000000000000*10^360*s^112+3.47593750000000000000000000000*10^358*s^111+1.26826250000000000000000000000*10^356*s^110+4.47647500000000000000000000000*10^353*s^109+1.52853750000000000000000000000*10^351*s^108+5.05230000000000000000000000000*10^348*s^107+1.61761250000000000000000000000*10^346*s^106+5.01626250000000000000000000000*10^343*s^105+1.50781250000000000000000000000*10^341*s^104+4.39270000000000000000000000000*10^338*s^103+1.24156250000000000000000000000*10^336*s^102+3.40548750000000000000000000000*10^333*s^101+9.06287500000000000000000000000*10^330*s^100+2.34108750000000000000000000000*10^328*s^99+1.43112500000000000000000000000*10^323*s^97+3.38603750000000000000000000000*10^320*s^96+2.95956250000000000000000000000*10^436*s^151+1.94535000000000000000000000000*10^285*s^83+1.17815000000000000000000000000*10^288*s^84+6.94362500000000000000000000000*10^290*s^85+3.98737500000000000000000000000*10^293*s^86+2.22566250000000000000000000000*10^296*s^87+1.20895000000000000000000000000*10^299*s^88+241.)):

``

A22 := factor(A11)

 

 

A22;

-0.693563083567816540884333863873e-1*(s^2+0.690931567083030562466492330164e-1*s+0.131419321908717817695986946802e-2)*(s^2+0.445612375757305814885907278258e-1*s+0.532301883634578204728036976352e-3)*(s^2+0.402062167667597088685587695446e-1*s+0.941585548325071550814140274653e-3)*(s^2+0.178146770743279479521546129578e-1*s+0.121430671829720613437700920041e-3)*(s^2+0.154359208565186521100458530794e-1*s+0.492005561434134764997390981824e-3)*(s^2+0.122918188745204333840354014128e-1*s+0.388300607745741227960683367216e-3)*(s^2+0.863140300873429921542456514714e-2*s+0.291745299251779255545001336604e-4)*(s^2+0.828039196507119580433769223220e-2*s+0.727497407245863812834187062561e-4)*(s^2+0.607436481046816728921279089128e-2*s+0.134186075347974245047612580380e-4)*(s^2+0.458436271424978357269227653428e-2*s+0.232926971678378773189431455766e-4)*(s^2+0.414343160903409582986969835712e-2*s+0.455962146303873008721781979537e-5)*(s^2+0.374405065332046635942464241096e-2*s+0.537483327064821956031752207744e-5)*(s^2+0.308164983295178642323181283804e-2*s+0.103733277114317462278221084800e-4)*(s^2+0.270802943601296759070390950912e-2*s+0.211692821958284031281117385160e-5)*(s^2+0.240487571179809311656093343448e-2*s+0.605020952058039325763244776641e-5)*(s^2+0.206479986964489326004678652256e-2*s+0.155476075193134071076241124565e-5)*(s^2+0.200840509898854997796302077352e-2*s+0.260946246260649583255193167511e-5)*(s^2+0.141228534442833303380857872480e-2*s+0.529884449346263671140458121116e-6)*(s^2+0.130712000477243093043988841410e-2*s+0.787853383708358871188211960532e-6)*(s^2+0.110177784064827936300275131990e-2*s+0.911242008117255852116992485023e-6)*(s^2+0.103187163101568514827753008784e-2*s+0.420305261905851372350608750704e-6)*(s^2+0.758315114625615756665156688672e-3*s+0.367667823547168692615446769931e-5)*(s^2+0.652292747049027721358424399256e-3*s+0.191381382285983212184664005136e-6)*(s^2+0.589019742399326438479733091708e-3*s+0.888354199408741547272716704796e-7)*(s^2+0.474011157834195062414679226940e-3*s+0.142134066686469366028973021414e-5)*(s^2+0.368704925499072007417610012226e-3*s+0.375283073644486365069107515211e-7)*(s^2+0.354065768414870835700456577106e-3*s+0.616994004710938357085510045428e-7)*(s^2+0.270874499590682457664184231874e-3*s+0.183627458369838140059179136241e-6)*(s^2+0.214257667169997892608885644276e-3*s+0.728296772582286004789761726844e-6)*(s^2+0.111790156467102191961969547699e-3*s+0.900560187553512296433072442719e-8)*(s^2+0.878122149954993762670378130746e-4*s+0.136152537358349528106243131608e-6)*(s^2+0.487870895223957494236428631748e-4*s+0.150633572915609877940132524299e-7)*(s^2+0.479344509868562361714905661004e-4*s+0.124381399298249577539121437149e-8)*(s^2+0.454745009556528016814010264950e-4*s+0.337342095880535738721894957562e-6)*(s^2-0.143495148631070416925092154015e-4*s+0.145327562110866754431115890171e-8)*(s^2-0.595822818080922880035242180570e-4*s+0.414101242962173843976655082214e-7)*(s^2-0.833725981684672214728149735802e-4*s+0.677211652338524488486045625395e-8)*(s^2-0.140387906073740752175949930393e-3*s+0.121373780663882262421998834705e-6)*(s^2-0.164996447035474522287137906342e-3*s+0.175759498261178224611587486501e-7)*(s^2-0.247000566237266392859399547006e-3*s+0.343450333978302545393470550672e-7)*(s^2-0.388855559426644543582883187630e-3*s+0.851288915863237525845855765743e-4)*(s^2-0.407053084377187413443204281254e-3*s+0.700825174589984777007770109154e-7)*(s+0.147492625368731563421828908555e-2)*(s^2-0.449673129088862465755109736440e-3*s+0.276805843066021138055480463079e-6)*(s^2-0.492699905118873316191938075884e-3*s+0.136388449252720773695296741225e-6)*(s^2-0.508808923376905529293029940652e-3*s+0.235668119381973081023203743495e-5)*(s^2-0.635141956505546615915396057152e-3*s+0.101906509778082892272747094459e-4)*(s^2-0.639735148104664748376377176506e-3*s+0.167334961281348149434227325357e-6)*(s^2-0.645279185772462312455444792726e-3*s+0.871332153627322180794720432979e-6)*(s^2-0.797634112705862241855550904846e-3*s+0.302452874138562695037670122153e-6)*(s^2-0.100042741323876254586674347312e-2*s+0.384718986726557567257741475414e-6)*(s^2-0.127039171105278636358307660493e-2*s+0.856696382581232762421213410957e-6)*(s^2-0.134771455852246605471767534215e-2*s+0.152420065081772472995156639087e-5)*(s^2-0.135770495872597390776930400277e-2*s+0.693417317477572748835432804749e-6)*(s^2-0.147244208677084907527616482611e-2*s+0.318793877550154042055406020451e-4)*(s^2-0.153212330954335657228430928450e-2*s+0.621260453203199004409122409386e-5)*(s^2-0.180354707053551114625973226922e-2*s+0.123996494870363933341841185121e-5)*(s^2-0.226074180278301586802878916464e-2*s+0.200454601566268590151369041903e-5)*(s^2-0.255599096743279306868943964272e-2*s+0.433862397451416380599892043761e-5)*(s^2-0.256211712686845918538017685252e-2*s+0.383410354327731205170010774555e-5)*(s^2-0.274564275487206542973516981530e-2*s+0.269302824761514647913892187093e-3)*(s^2-0.277244148166490505296776415940e-2*s+0.332706140963591415674374268505e-5)*(s^2-0.318752315819316191143469743584e-2*s+0.141765601219933423390422468084e-4)*(s^2-0.342282367304190991427403153066e-2*s+0.587417025759086227616500020290e-5)*(s^2-0.422971493182656200457720461336e-2*s+0.892951315804218011876289186209e-5)*(s^2-0.502856921561931027038854770516e-2*s+0.145924824583205710191246918636e-4)*(s^2-0.548467870154574502272103959628e-2*s+0.211381072866210365252746274844e-4)*(s^2-0.638386364361623612213700895240e-2*s+0.517102155120860738075309431215e-4)*(s^2-0.645854859703869051257346796042e-2*s+0.749308097208108689055801583941e-4)*(s^2-0.649579807834511122650589458614e-2*s+0.339323827955044515308100453611e-4)*(s^2-0.654857360078241600344582441518e-2*s+0.134903499583845910348322827849e-3)*(s+0.837509110705437027400183697928e-2)*(s+0.570259649669054252586463601053e-2)*(s+0.461771610357815283698142126909e-2)*(s+0.251077960036616416536656156893e-2)*(s+0.151796375059262599725980332730e-2)*(s+0.680456226667619818186507188592e-3)*(s+0.453629356701283760204795474245e-3)*(s+0.888577679798658000810594931464e-4)*(s+0.486721115037638963814063528697e-4)*(s+0.906329104271389581148284599671e-5)/(s*(s+0.327199108333099771478749074871e-1)*(s+0.203982055992343495590359155544e-1)*(s+0.117203099114620452924498515142e-1)*(s+0.727724549273518634231758482172e-2)*(s+0.622378218328195009869414856533e-2)*(s+0.288127931670050768303247450186e-2)*(s+0.201481086727237827074768415542e-2)*(s+0.259539877492647059098500228636e-3)*(s+0.211077748555797101608965607637e-3)*(s^2+0.515272961446762610202390218234e-1*s+0.958005764745424876600420647259e-3)*(s^2+0.405775310495930755768862053512e-1*s+0.476071731934613073587966338905e-3)*(s^2+0.263482304181727232535676854454e-1*s+0.679361216984965963948442975010e-3)*(s^2+0.112145772760791757855663118179e-1*s+0.107291036304092743306497293536e-3)*(s^2+0.106202304706343470547740389143e-1*s+0.487376200975008322068911203741e-4)*(s^2+0.953583263879931643445864438276e-2*s+0.403406969820213867662154460582e-3)*(s^2+0.652073727510838719951868378750e-2*s+0.117424119137545983511813369364e-4)*(s^2+0.636436278018703360412650434524e-2*s+0.182965863324330176600352710174e-4)*(s^2+0.454679030228354312447429261774e-2*s+0.778538966829913874322032209827e-5)*(s^2+0.363754879182606712168497329910e-2*s+0.529307355349755917454806278848e-5)*(s^2+0.294231716552996158632852575292e-2*s+0.241613912597391551349826209978e-4)*(s^2+0.276823092226373369594752636476e-2*s+0.106519108078968728422598708276e-4)*(s^2+0.260054320868396596243553512984e-2*s+0.286952277944613834344715048866e-5)*(s^2+0.243436527887358781187297993612e-2*s+0.151395629874495033476127245732e-5)*(s^2+0.242266746049997457824783378956e-2*s+0.201977445474704810893240246492e-5)*(s^2+0.185015267427942152001107961551e-2*s+0.941588507425162328452973364213e-6)*(s^2+0.169959735063261519091207336642e-2*s+0.264393789822114683593633156653e-5)*(s^2+0.127395019453773898281197207010e-2*s+0.421516592735940709480960132730e-6)*(s^2+0.108337606204351774210796814265e-2*s+0.816845875805346520906451698300e-6)*(s^2+0.966576530708079141254764545124e-3*s+0.282522690824929579698621523829e-3)*(s^2+0.871367686850150062443453042450e-3*s+0.330272413216146763394816198219e-6)*(s^2+0.775095698376549888941168763130e-3*s+0.165491226876338932099065532217e-6)*(s^2+0.773730272079490808965980902662e-3*s+0.152937140814993481927865445238e-6)*(s^2+0.682858055497174995749984398578e-3*s+0.705733975784660784051742955934e-6)*(s^2+0.519242830757910854832065304022e-3*s+0.317486247464500149215068029065e-5)*(s^2+0.395680382276048187254371908406e-3*s+0.124309170247910787575704014302e-6)*(s^2+0.347306986787839178223651559796e-3*s+0.337794468860758931015389046332e-6)*(s^2+0.300144123555496752643532340738e-3*s+0.210447792396506568995912493866e-5)*(s^2+0.210670876193828775650173034732e-3*s+0.140124381121470846895171120113e-7)*(s^2+0.186543128870002039199167906637e-3*s+0.849541463783836153340742478982e-7)*(s^2+0.164773101202727044990866209700e-3*s+0.211806381269268705665027611464e-6)*(s^2+0.115814026509466548020963292792e-3*s+0.372159348991986046910525061976e-8)*(s^2+0.790880549351762526113451020348e-4*s+0.645151836828898420663567738871e-5)*(s^2+0.438909089665102207812892266520e-4*s+0.334239261786490683553926355384e-7)*(s^2+0.392560690185253493079135255564e-4*s+0.174484775685372522405857081203e-8)*(s^2+0.352288075108914233668365221628e-4*s+0.103052736908851840621371966759e-5)*(s^2+0.155138316313735005695917967101e-4*s+0.481615386206920122711519556568e-9)*(s^2+0.277926036832391670874077517660e-6*s+0.316328087909844624572582764600e-8)*(s^2-0.651938424745803791628973788668e-4*s+0.144048985727731667159737768521e-7)*(s^2-0.755831461778435025233585158950e-4*s+0.708307658079019776065400385789e-8)*(s^2-0.175083063865879442021914355036e-3*s+0.105717675432312765725429311797e-6)*(s^2-0.184042746719725490062092959317e-3*s+0.211965137906444366068538697192e-7)*(s^2-0.299876171251979294935888714966e-3*s+0.498345303302121648926309714803e-7)*(s^2-0.343735786597367887222319836340e-3*s+0.356579824221405011980538479103e-6)*(s^2-0.379860637105873760210437970772e-3*s+0.724381948060352094676977694395e-7)*(s^2-0.535666434828203831405955866642e-3*s+0.525188674301704871442528628581e-6)*(s^2-0.551133268187637182461072080386e-3*s+0.133068726581012206067659502243e-6)*(s^2-0.648575827222084524375195244912e-3*s+0.239314952926844263894328081459e-6)*(s^2-0.736025161451458934694970235510e-3*s+0.114491885180938867279521240805e-5)*(s^2-0.778104300170502904975489664418e-3*s+0.256602553144173887923712269023e-6)*(s^2-0.880626440145307921819806664240e-3*s+0.144543153789098237182089345377e-4)*(s^2-0.103694198766105586671285756594e-2*s+0.440961771260250111695541477088e-6)*(s^2-0.128671629618371960575145002060e-2*s+0.680045389357376068966881695951e-6)*(s^2-0.144837352460465503944939937875e-2*s+0.715610392126514230868226912530e-5)*(s^2-0.151034680370896270414643307353e-2*s+0.101811776681834254125889617726e-5)*(s^2-0.159414078832314322615491286170e-2*s+0.249461059665038882065840827024e-5)*(s^2-0.178229210078585556190110727452e-2*s+0.141223614577947654317661019137e-5)*(s^2-0.179661448756953310116515895548e-2*s+0.159252798171102452477086166458e-5)*(s^2-0.242437818360575825278198645532e-2*s+0.245489333681429654745375904028e-5)*(s^2-0.254216928551367204107976005770e-2*s+0.145339491868644238479646321832e-4)*(s^2-0.300083023419162116629908146102e-2*s+0.379483865245216408115092615133e-5)*(s^2-0.339011022175993286952281662050e-2*s+0.470692845528993044693095365529e-4)*(s^2-0.354290338259809915794162704464e-2*s+0.575660059046203359502329998691e-5)*(s^2-0.391449610121541612324538142388e-2*s+0.779849627228849728718000391774e-5)*(s^2-0.395677824997954322083273223546e-2*s+0.101789875199088902744720140668e-3)*(s^2-0.463598443534434655525820487296e-2*s+0.113900686300167982280577577700e-4)*(s^2-0.498460071584213869172913912234e-2*s+0.160942605285866098849824148553e-3)*(s^2-0.547402264885373050977610875272e-2*s+0.174952713150723653605004984747e-4)*(s^2-0.644301339648887717502257504032e-2*s+0.282579353081156743433377874428e-4)*(s^2-0.676325753340368709367946301756e-2*s+0.789051136103869728436570417444e-4)*(s^2-0.706629995983822066227466813632e-2*s+0.466521151934957339485952693151e-4))

(2)

B11 := s*A22:

limit(B11, s = 0);

-0.431957414130948847103665564667e-1

(3)

with(inttrans);

[addtable, fourier, fouriercos, fouriersin, hankel, hilbert, invfourier, invhilbert, invlaplace, invmellin, laplace, mellin, savetable]

(4)

A33 := invlaplace(A22, s, t);

-0.465328823534752336947210931839e-7052*exp(-0.138963018416195835437038758830e-6*t)*(0.125933373860287822863803603363e7052*sin(0.562428801572070889775634078280e-4*t)+0.181822133813527135134040007223e7052*cos(0.562428801572070889775634078280e-4*t))+0.115612599157439491876754013728e-6866*exp(-0.932715644350010195995839533185e-4*t)*(0.161426154696394963566436986954e6866*sin(0.276142285147043476597455689710e-3*t)+0.436048948718899161108266375706e6866*cos(0.276142285147043476597455689710e-3*t))+0.196441467511478272604222722919e-6911*exp(0.377915730889217512616792579475e-4*t)*(0.890888259999922840037343752018e6910*sin(0.751988935041924871842839117242e-4*t)-0.459672234395989277350254256296e6910*cos(0.751988935041924871842839117242e-4*t))-0.183641610857742862333860636453e-6853*exp(0.338162876670184354683973150878e-2*t)*(0.481914892182810028372503314727e6851*sin(0.821399418642365131256105615302e-2*t)+0.203412086636475536502736307191e6852*cos(0.821399418642365131256105615302e-2*t))-0.442448778854946400086440037534e-6687*exp(0.121218909180287912639099322766e-2*t)*(0.298585882864655422819093302262e6686*sin(0.992718964525412715258424963906e-3*t)+0.802813537263461154704463389204e6685*cos(0.992718964525412715258424963906e-3*t))-0.583983870817200603541986321622e-6846*exp(0.353314997991911033113733406816e-2*t)*(0.140582305956723035754600851566e6845*sin(0.584542268898437525237425576733e-2*t)+0.464709269248615200323146484736e6844*cos(0.584542268898437525237425576733e-2*t))+0.375217371761822173262921668211e-6796*(-0.128284353558850778015238597256e6795*sin(0.199190269357738053881850544784e-2*t)+0.104629248182207837190909039180e6795*cos(0.199190269357738053881850544784e-2*t))*exp(0.195724805060770806162269071194e-2*t)+0.347011430020072792068496428395e-6973*(-0.175899015040020695647587361346e6972*sin(0.195108673254526612812659666623e-1*t)+0.104792070791784808379210579497e6972*cos(0.195108673254526612812659666623e-1*t))*exp(-0.476791631939965821722932219138e-2*t)-0.655907144636718138364700627570e-6745*exp(0.231799221767217327762910243648e-2*t)*(0.151891586334128119458074141540e6744*sin(0.245295346650278677612587596305e-2*t)+0.328479343612986724826704732111e6743*cos(0.245295346650278677612587596305e-2*t))-0.126678052777596414987758947130e-6764*exp(0.150041511709581058314954073051e-2*t)*(0.280489916296499292866497044281e6763*sin(0.124241423399868093168088734897e-2*t)+0.173091947394844971607902683821e6763*cos(0.124241423399868093168088734897e-2*t))+0.189970393522561458584351519294e-1*exp(-0.259539877492647059098500228636e-3*t)-0.194030538399994019036463824482e-1*exp(-0.211077748555797101608965607637e-3*t)-0.431957414130948847103665564667e-1+0.376025656526708185556820882893e-6939*exp(-0.196280345092626746539567627782e-4*t)*(0.100445010137063666680170841801e6939*sin(0.368725916929759465291825220478e-4*t)+0.189991992816462570234537167768e6939*cos(0.368725916929759465291825220478e-4*t))+0.711800269960535952138333020849e-6716*exp(-0.387547849188274944470584381565e-3*t)*(0.185171634049074292534790072123e6716*sin(0.123684645230849215397707110510e-3*t)+0.121121636157859239656357115171e6715*cos(0.123684645230849215397707110510e-3*t))+0.105805408765519797788330144737e-6914*exp(-0.202887655247965377884431026756e-1*t)*(0.510327842415991264204077015266e6914*sin(0.802731121948318023289556159656e-2*t)+0.459799914681424942528713540040e6914*cos(0.802731121948318023289556159656e-2*t))-0.478926432221900385732309542472e-6672*exp(-0.121133373024998728912391689478e-2*t)*(0.346055478134218793256628724082e6670*sin(0.743266472206098596933151331364e-3*t)+0.151833179312147703904248044644e6671*cos(0.743266472206098596933151331364e-3*t))+0.895493472173564770068303304334e-6788*exp(-0.531011523531717352738701945715e-2*t)*(0.656565659607569564424039931395e6786*sin(0.453214036467906309504323100432e-2*t)+0.281426450000747564266686221381e6787*cos(0.453214036467906309504323100432e-2*t))+0.133095385394816210180106410006e-6697*(-0.163664076449470508251841330224e6697*sin(0.452790059022145133352416279325e-3*t)+0.765076499959901445138875303378e6695*cos(0.452790059022145133352416279325e-3*t))*exp(-0.823865506013635224954331048500e-4*t)-0.774498542212907693394354889996e-6732*exp(0.797070394161571613077456430850e-3*t)*(0.165491946584947056726794408819e6731*sin(0.136355762012520236806234112476e-2*t)+0.216862419561360094143971887520e6731*cos(0.136355762012520236806234112476e-2*t))-0.214817830859973951024677161239e-6783*exp(0.322150669824443858751128752016e-2*t)*(0.371949445463108481940753209221e6782*sin(0.422845478884212383220287619608e-2*t)+0.287459711060602586354118041455e6781*cos(0.422845478884212383220287619608e-2*t))+0.380463661292665137755626626702e-6693*exp(-0.121718263943679390593648996806e-2*t)*(0.196178792509697198465352049232e6690*sin(0.180063102824065806212784884287e-3*t)-0.305527326551990629276240039983e6692*cos(0.180063102824065806212784884287e-3*t))+0.567614466027992815154599606243e-7027*exp(-0.395440274675881263056725510174e-4*t)*(0.509497787578785142539273570780e7026*sin(0.253967608922488909635583094550e-2*t)-0.563123916986653886256775806005e7026*cos(0.253967608922488909635583094550e-2*t))+0.212922638157864881483893128258e-6921*exp(-0.176144037554457116834182610814e-4*t)*(0.990780251785206727632658726480e6920*sin(0.101499610928754724277465725345e-2*t)-0.166840145046963610608747832497e6921*cos(0.101499610928754724277465725345e-2*t))-0.730954401551775857193901428583e-2*exp(-0.288127931670050768303247450186e-2*t)+0.877365567437703290483391772782e-2*exp(-0.622378218328195009869414856533e-2*t)-0.192234586963556052672845501912e-1*exp(-0.727724549273518634231758482172e-2*t)+0.129778567328976027966191784702e-6793*(-0.302630503103848161102650468571e6793*sin(0.540081807977106529387559092208e-4*t)+0.313058295719231692081633470668e6792*cos(0.540081807977106529387559092208e-4*t))*exp(-0.105335438096914387825086517366e-3*t)+0.734685669989477356527254599134e-6673*(-0.102538801273583512853067627728e6672*sin(0.125615756158433244816163802512e-3*t)+0.146626535388194737817681274664e6672*cos(0.125615756158433244816163802512e-3*t))*exp(-0.636975097268869491405986035050e-3*t)-0.636814909806545452373785274680e-1*exp(-0.117203099114620452924498515142e-1*t)+0.170536301197510358374831497924e-6798*exp(-0.219454544832551103906446133260e-4*t)*(0.462928352713640138927112381429e6797*sin(0.181500201669784506416497637928e-3*t)-0.218446017033888781075920272618e6798*cos(0.181500201669784506416497637928e-3*t))-0.196962313064969211659359120147e-6847*exp(-0.560728863803958789278315590895e-2*t)*(0.108766156760562297671520115779e6847*sin(0.870915325584553522729174271845e-2*t)+0.741091837646775409790898217255e6844*cos(0.870915325584553522729174271845e-2*t))+0.122571940249693951692873809868e-6739*exp(-0.326036863755419359975934189375e-2*t)*(0.219778352220308063372776710002e6739*sin(0.105470766613674033677334168272e-2*t)+0.179336479596878443368135753412e6739*cos(0.105470766613674033677334168272e-2*t))-0.994568221018990475991431059689e-6877*exp(-0.131741152090863616267838427227e-1*t)*(0.115767779883951980493478830581e6876*sin(0.224900846028352097487802045445e-1*t)+0.204556597699069898544449664541e6876*cos(0.224900846028352097487802045445e-1*t))+0.338459600489047827611453286286e-6820*exp(-0.259621415378955427416032652011e-3*t)*(0.675289782884221353194837663666e6819*sin(0.176279868258449449592057242000e-2*t)+0.397424832175805517822295386143e6819*cos(0.176279868258449449592057242000e-2*t))+0.335235519568688661308218774325e-6781*exp(-0.579070132547332740104816463960e-4*t)*(0.119406186455079321119857266833e6781*sin(0.191930014806440679569193999744e-4*t)-0.694698255526976313420822331737e6780*cos(0.191930014806440679569193999744e-4*t))+0.106202631408024735781216117860e-6690*exp(-0.435683843425075031221726521225e-3*t)*(-0.148316185803233937328416116085e6690*sin(0.374769264740988955984673801098e-3*t)+0.285307169729421150716057130688e6689*cos(0.374769264740988955984673801098e-3*t))+0.131173019592744008304813633486e-6804*(0.362017453150085387980866978017e6803*sin(0.664801268981783827376842741540e-2*t)-0.112769781667953295070655552841e6804*cos(0.664801268981783827376842741540e-2*t))*exp(0.169505511087996643476140831025e-2*t)+0.138825911190695126462789355761e-6792*(-0.962766729100988290460398772606e6791*sin(0.285837505113928509912871051594e-2*t)+0.621558942984099008488085670777e6791*cos(0.285837505113928509912871051594e-2*t))*exp(-0.318218139009351680206325217262e-2*t)-0.183460740170340015043805001426e-6898*exp(-0.257636480723381305101195109117e-1*t)*(0.249134671600592069629333575510e6897*sin(0.171534312238144531666331663082e-1*t)+0.270329712044712947945434765764e6898*cos(0.171534312238144531666331663082e-1*t))-0.106409426438456551448699786403e-6712*exp(0.891146050392927780950553637260e-3*t)*(0.108589369821586604987199887430e6712*sin(0.786190093201740849285171761724e-3*t)+0.671641560171187131782952955580e6712*cos(0.786190093201740849285171761724e-3*t))+0.323196151039921180214420518592e-6787*exp(-0.173653493393919589111825779898e-3*t)*(0.191391358796106563794128987559e6787*sin(0.554652082924824917347131356596e-3*t)+0.112514271064527622458077098335e6785*cos(0.554652082924824917347131356596e-3*t))+0.379460981697736523681715609776e-6730*exp(0.171867893298683943611159918170e-3*t)*(0.128333630055260457403649274359e6730*sin(0.571875206207156005575848516015e-3*t)-0.952376833604220513976455285371e6729*cos(0.571875206207156005575848516015e-3*t))-0.199350935303058840447650708027e-6782*exp(-0.197840191138024093627185954203e-3*t)*(0.137966527171553921858378121402e6782*sin(0.291836305175317213247866016319e-3*t)+0.195868660510416153730600120833e6782*cos(0.291836305175317213247866016319e-3*t))-0.608292428721856902090094122098e-6690*exp(-0.541688031021758871053984071325e-3*t)*(0.362140980881692749960106103898e6688*sin(0.723477679581835143618449535385e-3*t)+0.852576401135495282696204354729e6688*cos(0.723477679581835143618449535385e-3*t))+0.207123895759922563931714023915e-6765*exp(0.643358148091859802875725010300e-3*t)*(0.600122140273358818177475617014e6764*sin(0.515883400238066062596921115158e-3*t)+0.119730297245325607261843869105e6764*cos(0.515883400238066062596921115158e-3*t))-0.740947804486277854558168103080e-6865*exp(0.149938085625989647467944357483e-3*t)*(0.452152907029867389480112529989e6864*sin(0.165387728713546222272556914648e-3*t)+0.151401652584589191010423504934e6864*cos(0.165387728713546222272556914648e-3*t))-0.449599904433123605209267435269e-6817*exp(-0.150072061777748376321766170369e-3*t)*(0.100004541417280052667607644537e6817*sin(0.144289857586693922303396366188e-2*t)+0.469290473489852973913681415232e6816*cos(0.144289857586693922303396366188e-2*t))+0.687385062139039275884459104095e-6817*exp(-0.925076337139710760005539807755e-3*t)*(-0.241383974933548623736730128039e6816*sin(0.292954395579480123325740384080e-3*t)+0.207033841290379288726702569592e6816*cos(0.292954395579480123325740384080e-3*t))+0.870812929298647273391982975124e-1*exp(-0.203982055992343495590359155544e-1*t)-0.304292523787531212012561152247e-1*exp(-0.327199108333099771478749074871e-1*t)-0.141298181047584185335147048196e-2*exp(-0.201481086727237827074768415542e-2*t)-0.165393049052276862871581907891e-6793*exp(0.273701132442686525488805437636e-2*t)*(0.335790921554669652707903109537e6792*sin(0.316291642713358177522262716386e-2*t)+0.157589767950023203665278782809e6792*cos(0.316291642713358177522262716386e-2*t))+0.296708078930811604562907244449e-6695*exp(0.898307243784766550582579477740e-3*t)*(0.133481437019902474074433233630e6695*sin(0.886325040532445141100036657228e-3*t)+0.248342739859722970746331370654e6695*cos(0.886325040532445141100036657228e-3*t))-0.366623399305903918507206325254e-6764*exp(-0.386865136039745404482990451331e-3*t)*(0.358455376640523510020055388456e6764*sin(0.572058330237639604369260820230e-4*t)+0.653177615612162872520412971375e6763*cos(0.572058330237639604369260820230e-4*t))+0.257963370380733958947569858454e-6750*exp(0.440313220072653960909903332120e-3*t)*(0.143421991119057490865849885476e6750*sin(0.377629972951553122863753172171e-2*t)-0.229484540572315927100228670119e6749*cos(0.377629972951553122863753172171e-2*t))-0.308585585637123347529816465726e-6896*exp(0.249230035792106934586456956117e-2*t)*(0.231240038178433203188642467529e6895*sin(0.124390933838351984807684989838e-1*t)+0.645669135856110825090782148994e6895*cos(0.124390933838351984807684989838e-1*t))+0.912329277019821782648157386517e-6772*(-0.234387367242129729587078170775e6771*sin(0.469010486920027069585549735172e-2*t)+0.890300786659702229798347002419e6770*cos(0.469010486920027069585549735172e-2*t))*exp(-0.147115858276498079316426287646e-2*t)+0.193427312805841902040733690704e-6783*(-0.292943975678631249241453926559e6783*sin(0.673389962764010379761313979080e-3*t)+0.387796863932489982793948468863e6783*cos(0.673389962764010379761313979080e-3*t))*exp(0.267833217414101915702977933321e-3*t)+0.888866611415573770499597333800e-6713*exp(-0.849798675316307595456036683210e-3*t)*(0.354959171798970259141881543216e6711*sin(0.138628283898048585048305148978e-2*t)+0.207027780643951752245139491169e6712*cos(0.138628283898048585048305148978e-2*t))+0.112394809123040471313075938529e-6803*exp(0.324287913611042262187597622456e-3*t)*(0.126573517763348755223723525389e6803*sin(0.366268073974024460839372819190e-3*t)-0.214770573937809221808800783021e6803*cos(0.366268073974024460839372819190e-3*t))+0.629900330253865364805028551219e-7066*(-0.164255810099915347150065865204e7066*sin(0.205291413175336897770887167458e-4*t)+0.374647317756558842949340901118e7065*cos(0.205291413175336897770887167458e-4*t))*exp(-0.775691581568675028479589835505e-5*t)+0.357682406135697270039185479600e-6872*(-0.357560587662903164229211054511e6871*sin(0.115509044228679628236078892750e-3*t)+0.201416450471095448436790864552e6872*cos(0.115509044228679628236078892750e-3*t))*exp(0.325969212372901895814486894334e-4*t)+0.517857767931107670610523499813e-6718*exp(-0.130027160434198298121776756492e-2*t)*(0.153221597520403377038914392590e6717*sin(0.108573317826621841162979956234e-2*t)+0.434933673838959970897968477097e6717*cos(0.108573317826621841162979956234e-2*t))+0.268958319292336939467579623307e-6760*(-0.149627792678311941415080534550e6760*sin(0.767567713491985703938581033442e-3*t)+0.783398969598724398086191915749e6759*cos(0.767567713491985703938581033442e-3*t))*exp(-0.341429027748587497874992199289e-3*t)+0.533369074281957447224765547921e-6797*(-0.999466283864185659873060779669e6795*sin(0.239022502612652985490613817457e-3*t)+0.171100920582304724901462076997e6796*cos(0.239022502612652985490613817457e-3*t))*exp(0.275566634093818591230536040193e-3*t)+0.328057603809364545774838483459e-6713*exp(-0.181877439591303356084248664955e-2*t)*(0.131975402625298030554746029997e6712*sin(0.140894756831783453851755930591e-2*t)-0.116940616496294379005380877413e6712*cos(0.140894756831783453851755930591e-2*t))-0.129445013195254498241007475703e-6805*exp(0.920213733598627450310464796585e-4*t)*(0.927341695160661812130157409182e6804*sin(0.112821011498785897468839618775e-3*t)+0.466191761223814566952654916533e6804*cos(0.112821011498785897468839618775e-3*t))-0.200145106591611286269549711548e-6785*exp(-0.227339515114177156223714630887e-2*t)*(0.659215520198470957515103169623e6784*sin(0.161773426589913721154707968552e-2*t)+0.139170922681311091611295175230e6785*cos(0.161773426589913721154707968552e-2*t))+0.108341591216643367243739126000e-6738*(-0.125765945820385935040738089600e6738*sin(0.359420269570486945275824439408e-2*t)+0.546453622110725621662863811369e6737*cos(0.359420269570486945275824439408e-2*t))*exp(0.127108464275683602053988002885e-2*t)-0.170108591004670668831052327439e-6719*exp(0.177145169129904957897081352232e-2*t)*(0.204244208107659533628859850362e6718*sin(0.161819637122809367311344418563e-2*t)+0.480070252110055702957589926609e6717*cos(0.161819637122809367311344418563e-2*t))+0.141121750853159957765204568238e-6700*exp(0.518470993830527933356428782970e-3*t)*(0.747516040967906497147191571214e6698*sin(0.414909146460565240408744304047e-3*t)-0.101769610282339977731721690117e6700*cos(0.414909146460565240408744304047e-3*t))+0.739246153195778787275304448165e-6737*(-0.343269703192105751459614482627e6734*sin(0.100473160208932272090482284363e-2*t)+0.472310573128384832779329885885e6736*cos(0.100473160208932272090482284363e-2*t))*exp(0.368012580725729467347485117755e-3*t)+0.236577334845383687117450027788e-6764*exp(0.875415319329397210109571775180e-4*t)*(0.211668575789531525291725245962e6763*sin(0.313136001793385074892057514500e-3*t)+0.477674755949407599895135597838e6763*cos(0.313136001793385074892057514500e-3*t))+0.555511657985690309998376220680e-6720*(-0.327070376466714259420443396264e6718*sin(0.257520046881232422810843356558e-2*t)+0.179531918224826984504820203752e6719*cos(0.257520046881232422810843356558e-2*t))*exp(0.724186762302327519724699689375e-3*t)+0.113420467856928854375375582803e-7050*(0.111503453512799040640359858465e7050*sin(0.168014619386974972587456863096e-1*t)-0.129979990297876115955913608729e7050*cos(0.168014619386974972587456863096e-1*t))*exp(-0.483288265354039570627382272562e-3*t)+0.515615223343034383041454603438e-6765*(0.365015256334166460446016270951e6764*sin(0.190695225164174286981703786355e-3*t)-0.307145712695647409126217106644e6764*cos(0.190695225164174286981703786355e-3*t))*exp(0.189930318552936880105218985386e-3*t)+0.314462345098236640773066714835e-6859*(-0.465374419744439972668248082508e6858*sin(0.989322251186190133944993058442e-2*t)+0.401974916384439857042070327680e6858*cos(0.989322251186190133944993058442e-2*t))*exp(0.197838912498977161041636611773e-2*t)+0.285161685680084825491789771551e-6819*(-0.629975951521216465895204737319e6818*sin(0.324408658420543495447688816098e-3*t)+0.234877089180185601214456245140e6818*cos(0.324408658420543495447688816098e-3*t))*exp(0.389052150085251452487744832209e-3*t)+0.319550837900482420070431483887e-6724*exp(0.755173401854481352073216536765e-3*t)*(-0.913367242626541679831048579082e6723*sin(0.669201688543799597947682647905e-3*t)+0.256595116238063825103122911241e6723*cos(0.669201688543799597947682647905e-3*t))-0.591827335990771521368541406007e-6725*exp(-0.138411546113186684797376318238e-2*t)*(0.104298282889666604251490618420e6724*sin(0.295569538318017346310250253536e-2*t)+0.559506119415007181652797019101e6723*cos(0.295569538318017346310250253536e-2*t))

(5)

limit(A33, t = 0);

-0.693563083567816540884333863868e-1

(6)

limit(A33, t = 100);

-0.690665759941969127938176162995e-1

(7)

limit(A33, t = 1000);

-.264287267363129686538094956620

(8)

limit(A33, t = 10000);

-16197548572434.3018721820829842

(9)

plot(A33, t = 0 .. 800);

 

plot(A33, t = 800 .. 2000);

 

``


 

Download invlaplace_prime.mw

 

@AmusingYeti Thanks a lot for helping me out.  It really saves a considerable amount of time and memory. At present, I am not sure about the round of error cost. Once I have the final solution, I will get to know more about that.

Some more query about memory. I want to know more about the statement below :

"memory used=95.11GiB, alloc change=0.52GiB, gc time=67.67s". After computation, I see 389.MB memory at the right corner of the Maple worksheet. I guess it is RAM. What is the term memory used stand for? Is it virtual memory? what are gc-time and alloc change?

Regarding automation, some post suggests that restart should be used to clear any memory as unassigning the variable does not serve the purpose.

But even if I restart and append the same code again on the same worksheet, I don't see a memory clearance ( right corner of the worksheet still shows 389 MB memory).

 

@AmusingYeti Thanks for the suggestion. Will try to implement that.

thanks,

@AmusingYeti its 8 GB.

1 2 3 4 Page 2 of 4