Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hello!

How can I make MAPLE to create the solution of the following system?

https://math.stackexchange.com/questions/301068/how-do-you-find-a-corresponding-recurrence-relation-for-some-random-algorithm/301709

according to this link, how to parse or walk through the algorithm in maple to generate recurrence relation formula?

I have a differential equation involving several functions of the following form:

diff(h,z) = iAf + iBg,

where h, f and g are functions of the Cartesian coordinates x, y and R and the third coordinate corresponds to z = R for some fixed constant value R.  The derivative is then with respect to the coordinate z and A and B are constants, with i the usual imaginary unit.  Is there some way this equation could be solved explicitly with Maple?

Basically it spits out the subset of values for which a division by zero error will occur for the function you specify on  range you specify for each of it's arguments, but I get an ambigous error when ever exponentiation features in the function I specify, which of course dramatically reduces the application of the calculator. Division,addition,substraction and multiplication are currently the only available arithmetic operators availble for the function window that I know the error will not occur.

If some one can help it is much appriciated

 

DIVISION_BY_ZERO_CALCULATOR.mw

hello everyone,
   INGT.mw
 

L__d := 100:

L__b := 200:

L__c := L__d+(1/2)*L__b;

200

(1)

X := .3;

.3

(2)

Error, (in plot) procedure expected, as range contains no plotting variable

 

`ΔE__g` := 1.155*X+.37*X^2;

.3798

(3)

V__0 := .6*`ΔE__g`;

.22788

(4)

m__D := 0.67e-1*m[e];

0.67e-1*m[e]

(5)

m__B := (0.67e-1+0.83e-1*X)*m[e];

0.919e-1*m[e]

(6)

Sol := solve(2*cos(L__d*sqrt(E__x))+(m__D*sqrt((V__0-E__x)/E__x)/m__B-m__B*sqrt(E__x/(V__0-E__x))/m__D)*sin(L__d*sqrt(E__x))-(m__D*sqrt((V__0-E__x)/E__x)/m__B+m__B*sqrt(E__x/(V__0-E__x))/m__D)*sin(L__d*sqrt(E__x))*exp(-sqrt(V__0-E__x)*L__b) = 0, E__x);

0.1110897170e-2, 0.3531161505e-2, -.2585338615+0.9991335677e-27*I

(7)

 

0.1110897170e-2, 0.3531161505e-2, -.2585338615+0.9991335677e-27*I

(8)

E__1 := 0.1110897170e-2:

K__1 := sqrt(E__1);

0.3333012406e-1

(9)

K__2 := sqrt(V__0-E__1);

.4762027959

(10)

C := cosh((1/2)*K__2*L__b);

0.2399908351e21

(11)

beta := m__D*K__2/(m__B*K__1);

10.41631973

(12)

B := -beta*sinh((1/2)*K__2*L__b);

-0.2499821271e22

(13)

A := -B*sin(K__1*L__d)+C*cos(K__1*L__d);

-0.7112056933e21

(14)

h := proc (x) options operator, arrow; piecewise(x <= -L__c, A*exp(K__2*(x+L__c)), -L__c < x and x < -(1/2)*L__b, -B*sin(K__1*(x+(1/2)*L__b))+C*cos(K__1*(x+(1/2)*L__b)), abs(x) <= (1/2)*L__b, (1/2)*exp(K__2*x)+(1/2)*exp(-K__2*x), (1/2)*L__b < x and x < L__c, -B*sin(K__1*(x-(1/2)*L__b))+C*cos(K__1*(x-(1/2)*L__b)), L__c <= x, A*exp(K__2*(x-L__c))) end proc:

'h(x)' = h(x);

h(x) = piecewise(x <= -200, -7.112056933*10^20*exp(95.24055918+.4762027959*x), -200 < x and x < -100, 2.499821271*10^21*sin(3.333012406+0.3333012406e-1*x)+2.399908351*10^20*cos(3.333012406+0.3333012406e-1*x), abs(x) <= 100, (1/2)*exp(.4762027959*x)+(1/2)*exp(-.4762027959*x), 100 < x and x < 200, 2.499821271*10^21*sin(0.3333012406e-1*x-3.333012406)+2.399908351*10^20*cos(0.3333012406e-1*x-3.333012406), 200 <= x, -7.112056933*10^20*exp(-95.24055918+.4762027959*x))

(15)

L__y := 200:

L__z := 200:

P := proc (x, y, z) options operator, arrow; h(x)*cos(Pi*y/L__y)*cos(Pi*z/L__z) end proc:

'Psi(x, y, z)' = P(x, y, z);

Psi(x, y, z) = piecewise(x <= -200, -7.112056933*10^20*exp(95.24055918+.4762027959*x), -200 < x and x < -100, 2.499821271*10^21*sin(3.333012406+0.3333012406e-1*x)+2.399908351*10^20*cos(3.333012406+0.3333012406e-1*x), abs(x) <= 100, (1/2)*exp(.4762027959*x)+(1/2)*exp(-.4762027959*x), 100 < x and x < 200, 2.499821271*10^21*sin(0.3333012406e-1*x-3.333012406)+2.399908351*10^20*cos(0.3333012406e-1*x-3.333012406), 200 <= x, -7.112056933*10^20*exp(-95.24055918+.4762027959*x))*cos((1/200)*Pi*y)*cos((1/200)*Pi*z)

(16)

INGT := proc (x__i) `assuming`([evalf(int(int(int(P(x, y, z)^2*exp(-lambda*sqrt((x-x__i)^2+y^2+z^2)), x = -infinity .. infinity), y = -L__y .. L__y), z = -L__z .. L__z))], [0 < lambda]) end proc

evalf(INGT(2))

``

Warning,  computation interrupted

 

``


 

Download INGT.mw
 

L__d := 100:

L__b := 200:

L__c := L__d+(1/2)*L__b;

200

(1)

X := .3;

.3

(2)

Error, (in plot) procedure expected, as range contains no plotting variable

 

`&Delta;E__g` := 1.155*X+.37*X^2;

.3798

(3)

V__0 := .6*`&Delta;E__g`;

.22788

(4)

m__D := 0.67e-1*m[e];

0.67e-1*m[e]

(5)

m__B := (0.67e-1+0.83e-1*X)*m[e];

0.919e-1*m[e]

(6)

Sol := solve(2*cos(L__d*sqrt(E__x))+(m__D*sqrt((V__0-E__x)/E__x)/m__B-m__B*sqrt(E__x/(V__0-E__x))/m__D)*sin(L__d*sqrt(E__x))-(m__D*sqrt((V__0-E__x)/E__x)/m__B+m__B*sqrt(E__x/(V__0-E__x))/m__D)*sin(L__d*sqrt(E__x))*exp(-sqrt(V__0-E__x)*L__b) = 0, E__x);

0.1110897170e-2, 0.3531161505e-2, -.2585338615+0.9991335677e-27*I

(7)

 

0.1110897170e-2, 0.3531161505e-2, -.2585338615+0.9991335677e-27*I

(8)

E__1 := 0.1110897170e-2:

K__1 := sqrt(E__1);

0.3333012406e-1

(9)

K__2 := sqrt(V__0-E__1);

.4762027959

(10)

C := cosh((1/2)*K__2*L__b);

0.2399908351e21

(11)

beta := m__D*K__2/(m__B*K__1);

10.41631973

(12)

B := -beta*sinh((1/2)*K__2*L__b);

-0.2499821271e22

(13)

A := -B*sin(K__1*L__d)+C*cos(K__1*L__d);

-0.7112056933e21

(14)

h := proc (x) options operator, arrow; piecewise(x <= -L__c, A*exp(K__2*(x+L__c)), -L__c < x and x < -(1/2)*L__b, -B*sin(K__1*(x+(1/2)*L__b))+C*cos(K__1*(x+(1/2)*L__b)), abs(x) <= (1/2)*L__b, (1/2)*exp(K__2*x)+(1/2)*exp(-K__2*x), (1/2)*L__b < x and x < L__c, -B*sin(K__1*(x-(1/2)*L__b))+C*cos(K__1*(x-(1/2)*L__b)), L__c <= x, A*exp(K__2*(x-L__c))) end proc:

'h(x)' = h(x);

h(x) = piecewise(x <= -200, -7.112056933*10^20*exp(95.24055918+.4762027959*x), -200 < x and x < -100, 2.499821271*10^21*sin(3.333012406+0.3333012406e-1*x)+2.399908351*10^20*cos(3.333012406+0.3333012406e-1*x), abs(x) <= 100, (1/2)*exp(.4762027959*x)+(1/2)*exp(-.4762027959*x), 100 < x and x < 200, 2.499821271*10^21*sin(0.3333012406e-1*x-3.333012406)+2.399908351*10^20*cos(0.3333012406e-1*x-3.333012406), 200 <= x, -7.112056933*10^20*exp(-95.24055918+.4762027959*x))

(15)

L__y := 200:

L__z := 200:

P := proc (x, y, z) options operator, arrow; h(x)*cos(Pi*y/L__y)*cos(Pi*z/L__z) end proc:

'Psi(x, y, z)' = P(x, y, z);

Psi(x, y, z) = piecewise(x <= -200, -7.112056933*10^20*exp(95.24055918+.4762027959*x), -200 < x and x < -100, 2.499821271*10^21*sin(3.333012406+0.3333012406e-1*x)+2.399908351*10^20*cos(3.333012406+0.3333012406e-1*x), abs(x) <= 100, (1/2)*exp(.4762027959*x)+(1/2)*exp(-.4762027959*x), 100 < x and x < 200, 2.499821271*10^21*sin(0.3333012406e-1*x-3.333012406)+2.399908351*10^20*cos(0.3333012406e-1*x-3.333012406), 200 <= x, -7.112056933*10^20*exp(-95.24055918+.4762027959*x))*cos((1/200)*Pi*y)*cos((1/200)*Pi*z)

(16)

INGT := proc (x__i) `assuming`([evalf(int(int(int(P(x, y, z)^2*exp(-lambda*sqrt((x-x__i)^2+y^2+z^2)), x = -infinity .. infinity), y = -L__y .. L__y), z = -L__z .. L__z))], [0 < lambda]) end proc

evalf(INGT(2))

``

Warning,  computation interrupted

 

``


 

Download INGT.mw

 

I'm trying to calculate a triple integral complicated by a procedure that changes each time a variable xi, while the program takes a lot of time and it gives me the message "Warning, computation interrupted". If anyone can help me I will be very happy

Is there something I should be doing whenever I use simplify to avoid things like this, or should I stop using the "is" function all together?

 

interface(showassumed = 0):

 

sum(binomial(k+j, k), j = 0 .. n-k) = binomial(n+1, k+1)

(n-k+1)*binomial(n+1, k)/(k+1) = binomial(n+1, k+1)

(1)

#And we have:
is(sum(binomial(k+j, k), j = 0 .. n-k) = binomial(n+1, k+1))

FAIL

(2)

#And since:
is(simplify(convert(sum(binomial(k+j, k), j = 0 .. n-k) = binomial(n+1, k+1), 'factorial')))

true

(3)

is(sum(binomial(k+j, k), j = 0 .. n-k) = binomial(n+1, k+1)) = is(simplify(convert(sum(binomial(k+j, k), j = 0 .. n-k) = binomial(n+1, k+1), 'factorial')))


 

Download main.mw

Hello people :) 

As the captian says, im trying to remove an old task i've made.
But i get this:

Error in Get, invalid object [_XML_reply_data_get("reference" =
"_Maplets_reference_12","parameter" =
"value",_XML_content("Task,UserTasks,Nyops",&Entity "#xc3",&Enity
"#xa6","tning"))]

And i have no idea what it is, but it won't erase my task :'D

Thanks a bunch in advance! 

Have a great weekend you all
Best regards Lucas :)

Hi 

I was wondering someone could help me with something.

Basically I've got a function of f(x,z) that I have plotted which comes from an integration that was taken over a large number of small steps. I have been able to plot this function to determine a Potential Field as a function of x and z. (No Problems so far)

However the next step that I want to be able to do is compute a field form this potential. apologies for my lack of knowing how to input into here as I am fairly new to this. To do this I have:

Field = - Nabla * Potential 

Again no problem, I can plot a gradient field from it which is great using the gradplot command. I beleive this comes as a result of the following:

gradplot = [diff(f,x), diff(f,z)]

HOWEVER I want to be able to plot a Field Strength plot (In essence a Gradient Field Strength) where it takes the magnitude of the vectors at each point. Firstly am I correct in saying that this should be:

Magnitude Field = |F| = sqrt( Fx^2+Fz^2)      where Fx is diff(f,x) and Fz is diff(f,z)

 

However when I try to compute this, it takes an age to evaluate and doesn't seem to make much progress. I think this may come as a result of computing square roots of large numbers/large decimal places? I was wondering whether there is any command that anyones knows of that would calculate Gradient Field Strength/Magnitude that may solve this issue? Or even who to just compute this Square root?

A few seconds after calling up Help starts zucking araound and the whole computer then freezes. Ctrl-Alt-Delete doesn't work, hard reset required. Very funny. Am I alone?

Hi good afternoon , im looking for maple code of HPM to solve time dependent diffusion reaction , can anybody help me . I didnt get how to solve it by maple . 

Hello

I know that anames( 'user' ) recovers the contents of the workspace but I wonder whether it would be possible to recover only the contents of each m-file read within a worksheet.  

Many thanks

Ed

 

My maplet

mpt := Maplet(Window("Точка выше или ниже прямой",
[[Plotter[f]()], ["1st vector"],
[TextField[v1](), TextField[v2](), TextField[v5]()],
["2nd vector"],
[TextField[v3](), TextField[v4](), TextField[v6]()], ["Скаляр"],
[TextField[scalar]()], [Button("Add", Evaluate(f = 'work(1)'))],
[[ToggleButton['d2']("2D", 'value' = true, 'group' = 'BG1')],
[ToggleButton['d3']("3D", 'value' = false, 'group' = 'BG1')]]]),
ButtonGroup['BG1']())

 

Now I choose via toogle button what to plot 2D or 3D addition but I need to plot it in two different plotter, is there a way to do it?

 

Procedure 

with(Maplets);
with(Elements);
with(Tools);
[AddAttribute, AddContent, Get, ListBoxSplit, Parse, Print, Set, 

  SetTimeout, StartEngine, StopEngine]
with(plots);
with(LinearAlgebra);
with(Student[LinearAlgebra]);

with(Maplets[Elements]);

work:=proc(g)

     local x1, y1, x2, y2, scalar;z1, z2,d2,d3,u,v,U,V,r1,r2,r3,R1,R2,R3,V1:
 x1 := Get(v1::algebraic);
 y1 := Get(v2::algebraic);  

z1 := Get(v5::algebraic);  

x2 := Get(v3::algebraic);  

y2 := Get(v4::algebraic);  

z2 := Get(v6::algebraic);  

scalar := Get(scalar::algebraic);

 d2:= Get(d2);  

d3:= Get(d3);    

 if d2 = true then  u, v:=<x1, y1>, <x2, y2>;  

elif d3 = true then   u, v:=<x1, y1,z1>, <x2, y2,z2>;  

end if;  U:=arrow(u, color=red):  

V:=arrow(u, v, color=black):    

   r1:=u+v: r2:=u-v: r3:=ScalarMultiply(u,scalar):r4:=ScalarMultiply(v,scalar):  

R1:=arrow(r1, color=red):

 R2:=arrow(v, r2, color=black):  

R3:=arrow(r3, color=black):

 R4:=arrow(r4, color=black):  

V1:=arrow(v, color=blue):    

if g = 1 then return VectorSumPlot(u,v, show=1, caption="Треугольник", vectorcolors=[red,blue]);    

end if;  

end proc:

Repeat once again what do I need: I want to plot addition of vectors in 2D and 3D in 2 different plotters at once

how to convert a function to transition matrix?

how to convert a differential system to transition matrix?

 

Does global optimization toolbox has randomization option it gives same value every time

First 117 118 119 120 121 122 123 Last Page 119 of 368