dorna01

15 Reputation

One Badge

3 years, 164 days

MaplePrimes Activity


These are questions asked by dorna01

Hi everybody

I try to save an array of equations in a "mpl file", and then read file and array by command read mpl file. But system error "bad id" occurs. Previously I saved this mpl file by "save" and read array by "read" command. I ran this code several times without any problem, suddenly error "bad id" occurred. When I save mpl file by the code edit region and then read it by "read" in the same or  a new worksheet, the error doesn't happen!
 What is wrong?
"Id" refers to the array index?

Thanks.
Maple 18 and Windows 10

Hi everybody
I have some problems with fsolve(complex equation). It results some answers (I expect answers in the range 10e6 to 10e11) but substitution them into the main equation leads to numbers of order 10e-8 to 10e8. I know fsolve solves equation numerically, so 10e-8 t0 10e-6 is acceptable, but what about 10e7? How can I handle this problem? I have an Array of this kind of equations to solve and then analyze answers.
How can I increase the speed of calculations? I try to do some parallelization (thanks dohashi for posts about parallel programming) but I couldn't do. I upload the code below.

Thanks.

EQ1 := 1.780876811*10^90*(-(1.857495893*10^(-32)*I)*(-(.9215096529*(-1.077177489*10^(-57)*omega^2+1.251444314*10^(-43)-7.423792254*10^(-74)*omega^4))*(1.042248387*10^(-7)*omega-3.773917830*10^(-22)*omega^3)+1.022012860*10^(-43)-9.365146438*10^(-58)*omega^2+1.290731820*10^(-74)*omega^4+8.072440803*10^(-47)*omega^2*(7.038725244*10^(-13)-9.109383000*10^(-28)*omega^2))*exp(-.9800000000*I-4.717786244*10^(-17)*omega^2)-(1.857495893*10^(-32)*I)*((.9215096529*(5.411991727*10^(-58)*omega^2-1.370413754*10^(-43)+1.063387455*10^(-73)*omega^4))*(1.042248387*10^(-7)*omega-3.773917830*10^(-22)*omega^3)-1.119171234*10^(-43)+3.850718130*10^(-58)*omega^2+1.279097989*10^(-74)*omega^4+1.703871878*10^(-48)*omega^2*(5.154059190*10^(-14)+3.036461000*10^(-28)*omega^2)+8.072440803*10^(-47)*omega^2*(7.038725244*10^(-13)+9.109383000*10^(-28)*omega^2))*exp(.9800000000*I-4.717786244*10^(-17)*omega^2)+2.054040475*10^(-31)*((1.936145393*10^(-59)+1.043762907*10^(-58)*I)*omega^2+4.297601656*10^(-46)-1.690952584*10^(-44)*I+(-1.159596547*10^(-75)+1.164619044*10^(-74)*I)*omega^4)*exp(-4.717786244*10^(-17)*omega^2)*(1.042248387*10^(-7)*omega-3.773917830*10^(-22)*omega^3)+(2.799879047*10^(-71)*I)*(-6.704964363*10^(-12)-3.118737242*10^(-28)*omega^2)*omega*exp(.9800000000*I-1.090999486*10^(-14)*omega^2)-(2.799879047*10^(-71)*I)*(8.281232388*10^(-12)+2.177273887*10^(-28)*omega^2)*omega*exp(-.9800000000*I-1.090999486*10^(-14)*omega^2)+3.476335242*10^(-51)*((.1388433141*I)*(-2.893776471*10^(-25)-1.303697368*10^(-38)*omega^2+7.808106616*10^(-55)*omega^4)+4.959435112*10^(-25)-3.098806468*10^(-39)*omega^2-3.391707726*10^(-55)*omega^4-1.314961283*10^(-30)*(-2.854029409*10^(-11)+1.827522021*10^(-27)*omega^2)*omega^2)*exp(-4.717786244*10^(-17)*omega^2)-2.814230381*10^(-37)*(9.949004410*10^(-35)*(-6.832852706*10^(-13)-1.621609260*10^(-14)*I-(2.889900216*10^(-30)*I)*omega^2-(.9082907587*I)*(8.002616800*10^(-12)-1.954522389*10^(-30)*omega^2)-(.4487255373*I)*(9.612550267*10^(-12)+9.109383000*10^(-28)*omega^2)+4.081082866*10^(-29)*omega^2)*exp(-1.090999486*10^(-14)*omega^2)-1.995292057*10^(-54)*omega)*omega)/omega^2

Test_MaplePrime971127.mw

Hello

I have problem to calculate a numerical integration. I guess, I can't use integration cammands correctly or singularity is the problem.
The equation has a singularity at y=0 (exp(1E-9*I*x/y)). The limits are semi_infinite.
integration:

f1 := (-(1.671130827*10^(-21)*I)*((0.6132469529e-1-.5253537767*I)*x^3+(.2894208344-2.479398027*I)*y*x^2+(0.1107017247e-1+(-.2377238011-0.2774956673e-1*I)*y^2)*x+(-.1167302837+.9999999999*I)*y^3+(0.2912028092e-1-.2494663781*I)*y)*exp(-1.732070784*10^(-15)*x*y)*((-0.6132469529e-1-.5253537767*I)*x^3+(.3115611798+2.669068985*I)*y*x^2+(-0.1107017247e-1+(-.3699321954+0.4318229002e-1*I)*y^2)*x+(.1455017102+1.246477826*I)*y^3+(0.2512356422e-1+.2152274756*I)*y)*exp(1.732070784*10^(-15)*x*y)*exp(-2.178669712*10^(-15)*y^2-3.534838336*10^(-16)*x^2+(1.000000000*10^(-9)*I)*x/y)/x)*y:

evalf(Int(f1,y=0..infinity));



Thank you
Test4_MaplePrimes.mw

 




 

Hi
I am confused by solve and SolveEqutions and their options!
Solving same equation results different solutions.
Which  solution is correct and why?
How i can be sure one of the solutions is correct about the other equations?
I attached my code.

Thanks.

Test2_MaplePrimes.mw

Hi all

How I use "solve" or "fsolve" for this equation ?

M2 := evalf[4](Matrix(4, 4, {(1, 1) = BesselJ(0, 0.5e-1*sqrt(0.1111111111e-16*omega^2-25.00027778)), (1, 2) = -BesselJ(0, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25)), (1, 3) = -BesselY(0, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25)), (1, 4) = 0, (2, 1) = (0.1111111111e-16*I)*omega*(1-25000000000000/omega^2)*BesselJ(1, 0.5e-1*sqrt(0.1111111111e-16*omega^2-25.00027778))/sqrt(0.1111111111e-16*omega^2-25.00027778), (2, 2) = -(0.4444444444e-16*I)*BesselJ(1, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (2, 3) = (0.4444444444e-16*I)*BesselY(1, 0.5e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (2, 4) = 0, (3, 1) = 0, (3, 2) = BesselJ(0, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25)), (3, 3) = BesselY(0, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25)), (3, 4) = -BesselY(0, 0.60e-1*sqrt(0.1111111111e-16*omega^2-25)), (4, 1) = 0, (4, 2) = (0.4444444444e-16*I)*BesselJ(1, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (4, 3) = (0.4444444444e-16*I)*BesselY(1, 0.60e-1*sqrt(0.4444444445e-16*omega^2-25))/sqrt(0.4444444445e-16*omega^2-25), (4, 4) = -(0.1111111111e-16*I)*omega*BesselY(1, 0.60e-1*sqrt(0.1111111111e-16*omega^2-25))/sqrt(0.1111111111e-16*omega^2-25)})):


with(LinearAlgebra):
DETM2 := Determinant(M2):
solve(DETM2 = 0, omega);


Error, (in solve) cannot solve for an unknown function with other operations in its arguments

Is this Error because of combination of bessel functions? if I use asymptatic forms, does it work?

Thanks

1 2 Page 1 of 2