f := x^2*(y/x+sqrt(-7*y^2/x^2))/(y^2*(x/y+sqrt(-7*x^2/y^2)));
v := parametrization(f, x, y, t);

it can not parametrize.

i do not know which book teach group theory and algebraic curve

can we call this algebraic curve over finite field ?

how to represent a function as an algebraic curve equation for parametrization?

There is an heuristic-free algorithm, designed by Greg Reid to find structure constant of finite dimension Lie algebra of symmetries, admitted by differential equation. Details could be found in his paper:

http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=2316480

The main problem - I couldn't find implementation of it in Maple...

When I do simplify(LegendreP(n, 1, cos(t))), Maple gives me -sqrt(1-cos(t))*sqrt(cos(t)+1). Isn't it the same thing as -sin(t)? How can I have Maple further convert/simplify it to -sin(t)? (I tried simplify(%, trig) but it didn't work). I am new on Maple. Thanks in advance for anyone's help!

Trying to use the Explore() to parametrize a couple of plots simultaneously.

According to the help: "A plotting command to be explored can be a function call to plot, plot3d, a command in the plots package, or any user-defined procedure which returns a two- or three-dimensional plot structure or an Array of such plot structures."

However with the definitions

restart;
with(plots):
with(plottools):
Parab1 := ''plot(a*x^2, x = -1 .. 1, y = -3 .. 3)'';
Parab2 := ''plot(a*x^2+1, x = -1 .. 1, y = -3 .. 3)'';
L:=Array(1..2);
L[1]:=Parab1;L[2]:=Parab2;

I then find the following:

Explore(L[1], parameters = [a = -1.0 .. 1.0]);
gives the expected parabola, and the slider updates the parabola dynamically, as expected, and of course the same with L[2] in place of L[1]. However,

Explore(L, parameters = [a = -1.0 .. 1.0])

does not work, giving "Warning, expecting only range variable x in expression a*x^2 to be plotted but found name a". If I use display(L) in the command above, I receive two adjacent plots, as expected but they do not update as the parameter a is changed via the slider.

I'd guess I need to delay or promote variable evaluation somehow, but I can't seem to make it work. What am I missing?

Hello,

This is the system of equations in term of sin and cos. I have used the command "solve" in Maple but it yielded only 2 solutions. I've tried to use with(RealDomain): It yielded more solutions but most of them were wrong.

f1 := -8100+(-30+70*cos(t1)-40*cos(t2))^2+(-70*sin(t1)+40*sin(t2))^2

f2 := (-20-80*cos(t3))^2+(-15+70*cos(t1)+10*cos(t1+t))^2+(-70*sin(t1)-10*sin(t1+t)+80*sin(t3))^2-5625

f3 := (-20-80*cos(t3))^2+(15+40*cos(t2)+10*cos(t1+t))^2+(-40*sin(t2)-10*sin(t1+t)+80*sin(t3))^2-5625

f4 := 10*cos(t1+t)*(30-70*cos(t1)+40*cos(t2))-10*sin(t1+t)*(70*sin(t1)-40*sin(t2))

Anybody know how to solve this system of equations to get the full set of roots?

Thank you very much in advance.

for i from 0 to 3 do print("rawData[", i, "] = new double[] { ", hello(i+1), ",", data[i+1], "};") end do

wrong format after copy to notepad from the maple 15

do not know why it insert a empty line, and make "};" in another line

and there is unexpected " 

wrong format example:

"rawData[", 0, "] = new double[] { ", 73.25, ",", 0.1510425143,

"};"
"rawData[", 1, "] = new double[] { ", 73.15, ",", 0.3974080269,

"};"
"rawData[", 2, "] = new double[] { ", 72.85, ",", 0.4661517269,

"};"
"rawData[", 3, "] = new double[] { ", 73.25, ",", 0.3974080269,

"};"

expected format:
rawData[0] = new double[] { 25.0, 20.0 };
rawData[1] = new double[] { 27.0, 34.0 };

I have a expression: (x-2)/(x-2) and I want to retrieve all the variables of this expression.

I tried the following commands:

1. indets(expr, name);

2. indets(numer(expr),name);

but all i get is an empty set...it seems that maple will first simplify the expression and then figure out what variables are related to the given expression.

Is there any workaroud to achive what i want?

thanks

Hello,

I am trying to medel a mechanical system using parametric relations. Unfrotunately in some cases, the relations and equations become too lengthy. Maple calculates the relations but when I try to convert them to MATLAB language an error message appears. I have uploaded corresponding file in maple. As can be seen in this file, an error message appeard after trying to convert to MATLAB language. In some cases, for example C22, the relations are too length so that maple cannot show them. Do you have any idea to overcome this problem?

Thanks in advance.

P1.mw

Hello,

I have just updated maplesim to the last version (maplesim2015).

On the multibody examples of the maplesim2015 library, I didn't have any problems.

However, when i try to launch a simulation on a model i have made with a previous version (which worked before I use maplesim2015), i meet this problem :

May you help me to make the troubleshooting ? For the moment, I didn't manage to locate the error.

Do you have some ideas about this mistake ?

Thanks a lot for your help.

Hello, dear experts.
I have a question...
I solve the system of differential equations,where one of the initial conditions need to be chosen so thatcondition is metat the end of integration.
The task is not difficult, but I'm having trouble with the syntax.

1.I can't "pull"the desired function from the solution and find its value at a certain point.
I try to do so:
r_ravn:=s->subs(F,r(s));
evalf(r_ravn(s_end));
evalf(r_ravn(0));
but there is no result
2.In this case,instead of"for"it is better to use a while loop, but again the problem arises 1.
Tell me, please,how to implemen my program.

restart:
R:=0.3:
theta_min:=Pi/6:
theta_max:=Pi/2:
betta_max:=evalf(Pi/180*80);
p:=2*10^5:

theta0:=s->Pi/3/s_end*s+Pi/6:
r0:=s->R*sin(theta0(s)):
s_end:=evalf(R*(theta_max-theta_min)):

sol1:=solve({sin(betta_max)=c/r0(0)},{c});
const1:=0.1477211630;

betta0:=s->arcsin(const1/r0(s)):
betta:=s->arcsin(r(s)/r0(s)*sin(betta0(s))):
A:=s->cos(betta(s))/cos(betta0(s)):
T1:=s->rT1(s)/r(s):
T2:=s->T1(s)*tan(betta(s))^2:

step:=0.001:
delta:=0.001:
for i from 1 to 3000 do
r_min:=0.3-step:
rT1_n:=p*Pi*r_min^2/2/Pi/sin(theta_min):

sys := diff(rT1(s), s)-A(s)*T2(s)*cos(theta(s)),diff(theta(s), s)-A(s)/T1(s)*(p-T2(s)*sin(theta(s)/r(s))),diff(r(s),s)-A(s)*cos(theta(s)),diff(z(s),s)-A(s)*sin(theta(s));
fcns := {rT1(s),theta(s),r(s), z(s)};
F := dsolve({sys,rT1(0)=rT1_n, theta(0)=theta_min,r(0) = r_min, z(0) = 0}, fcns, numeric,output=listprocedure):
r_ravn:=s->subs(F,r(s)):
if abs(evalf(r_ravn(s_end))-R)=delta then break:
print(r_min):
end if:
end do:

r_ravn:=s->subs(F,r(s));
evalf(r_ravn(s_end));
evalf(r_ravn(0));
plot([r_ravn(s),r(s)],s=0..s_end);

Hello,

i did see that u can calculate minimial polynomials for for example \sqrt(2) or smth.

Suppose now u have a ground field K and an extension L=K(p,q) where p,q are some given rational functions in the indeterminates x,y, now i would like to calculate the minimal polynomial for lets say x+y or x over L which then obviously splits over the field K(x,y).

Thx in advance

(in dsolve/numeric/bvp) initial Newton iteration is not converging.

Dear Maple users

I am unsure how to handle events and their probabilities in Maple. Let's say I know that an event A has the probability say 0.3 and another event B has probability 0.8. I would like to make the following assignments:

P(A):=0.3:

P(B):=0.8

and maybe defining the conditional probability:

P(A|B):=0.55

but I am not allowed to do so in Maple because if will regarded as a function definition. My purpose is to make simple calculations with those probabilities for example:

P(C):=P(A)*P(B)  etc.

My problem is therefore more of a notational problem than a mathematical one. I hope someone can advice me on a proper setup. I am using 2D math notation, by the way. I could of course name the variables containing the probability using simple names like X1, X2, etc., but then I need constantly to remember what they really mean. The above assignments would be much better, because they are easier to handle mentally.

Regards,Erik

Trying a few integrals in Maple -- doesn't seem to handle them very well. Any workarounds please?

f1(p,b):= 1/(p^2 + b^2)^2 ;

f2(p,b):= exp(-p/b)^2;

T1pbm := Int(q*f1(q, b)*ln((p+q)^2+m^2), q = 0 .. infinity);

T2pbm := Int(q*f1(q, b)*ln((p-q)^2+m^2), q = 0 .. infinity);

with(IntegrationTools):

PT1 := evalf(Parts(T1pbm, q*f1(q, b), ln((p+q)^2+m^2)));

PT2 := evalf(Parts(T2pbm, q*f1(q, b), ln((p-q)^2+m^2)));

Ev := evalf(-pi*c*Int*(p*f1(p, b)*(PT1-PT2), p = 0 .. infinity))

The following toy procedure does not work (I give a value to a symbol). How to make it work ? Thanks.