How do you move plot labels individually? Just click-hold-and-drag does not work.

"A Virtual 3D Solar System": http://www.maplesoft.com/applications/view.aspx?SID=4484&view=html#mapleautobookmark5

How do you color the orbits individually?

I am also not sure how to program a check in the code to determine that dsolve solved the differential equation.

For example

ode:=diff(y(x),x)-a*(x^n-x)*y(x)^3-y(x)^2=0;
sol:=dsolve(ode,y(x));

In this case, sol is () as it could not solve it. When I tried odetest, I get an error

odetest(sol,ode);

Error, invalid input: odetest uses a 2nd argument, ODE, which is missing

What is the correct way to check dsolve was successful that will work for all cases? I am looking for programmable method, no GUI use.

In Maple help on int, it says

"If Maple cannot find a closed form expression for the integral, the function call itself is returned"

But then, what is the correct way to check for this in the code? How do I know that the result returned is the call I made? (it would have been easier if these functions throw an error, or set some status code that one can check for success or failure).

Hello friends;

I uploaded the file. Firstly i solved the 4. order ode. I converted it to trigonometric functions from exponantial.

I have 4 equation and 4 unkowns. I thought , it is sufficient to solve this. As you see , i can not solve the system .

May you help me please where i am doing wrong. 

Thanks.

odev4.mw

(PS:i downloaded the file again correctly)

How do you plot different kinds of data simultaneously in one graph? E.g. surfdata and spacecurve?

I am trying to do a substitution as shown in the attached document. I know variants of this question have been asked before but dont quiet get what to do. It is problem with algsubs and how it handles denominators I think. Can get substiturion to work for simple fractions but more complicated ones fail. Would appreciate any guidance here.

Download UHG5_substitution.mw

Dear Maple experts,

I want to solve the following problem with Maple. I appreciate if you can help.

n:=1+(2/Pi)*evalf(Int((x*0.187549975)/(x²-10.69²), x=0..infinity));

Hi,

How do you plot a list of values/an array into a surface?

E.g.

t = [seq(0+i*(2*evalf(Pi)*(1/10)), i = 0 .. 10)]

x:=[3,4,2*evalf(Pi), 7.83]:

y:=[2.5,4.3,6,2*evalf(Pi)+2]:

z:=[-2,0.3,1.5,evalf(Pi)]:

plot3d([x[i]*cos(t), x[i]*sin(t), y[i]*cos(t-z[i])],-10..10,-10..10);

Hi,

  Is there any benchmark of the performance of maple on windows 7 64 bit vs linux, especially in solving generalized eigenvalue problem?

Thank you very much

Hello,

I use maple for the determination of the symbolic relations between the geometric relations on multibody system.

I would like to create a dll function for the functions that I have determined symbolically.

I have already a paper presenting : 
- the transformation Maple --> code C with the package with(CodeGeneration)
- and the transformation code C --> another software.

I would like to directly create a dll function from my maple function with using code C. Is it possible and how ?

Thanks a lot for your help.

Hi all!

I do a small calculation and get a system of 6
nonlinear equations.
And "n" is the degree of the equation is float.

Here are the calculations that lead to the system.

restart;
 with(DirectSearch):
 B:=1: 
 q:=1: 
 l:=1: 
 n:=4.7:
 V:=0.05:
 N:=1200:
 
 
 kappa:=Vector(N+1,[]):
 theta:=Vector(N+1,[]):
 u:=Vector(N,[]):
 M:=Vector(N,[]):
 Z:=Vector(N,[]):
 
 M_F:=q*(6*l*(z-l)-z^2/2):
 M_1:=piecewise((z<l), l-z, 0):
 M_2:=piecewise((z<2*l), 2*l-z, 0):
 M_3:=piecewise((z<3*l), 3*l-z, 0):
 M_4:=piecewise((z<4*l), 4*l-z, 0):
 M_5:=piecewise((z<5*l), 5*l-z, 0):
 M_6:=6*l-z:
 M_finish:=(X_1,X_2,X_3,X_4,X_5,X_6,z)->M_1*X_1+M_2*X_2+M_3*X_3+M_4*X_4+M_5*X_5+M_6*X_6+M_F:
 
 
 kappa_old:=0:
 theta_old:=0:
 u_old:=0:
 M_old:=0:
 
 
 step:=6*l/N:
 u[1]:=0:
 kappa[1]:=0:
 theta[1]:=0:
 
 
 
 
 for i from 2 to N do
 
 z:=i*step:
 kappa_new:=kappa_old+B/V*(M_finish(X_1,X_2,X_3,X_4,X_5,X_6,z))^n*step:
 
 theta_new:=theta_old+1/2*(kappa_old+kappa_new)*step:
 
 u_new:=u_old+1/2*(theta_old+theta_new)*step:
 
 Z[i]:=z:
 kappa[i]:=kappa_new:
 theta[i]:=theta_new:
 u[i]:=u_new:
 kappa_old:=kappa_new:
 theta_old:=theta_new:
 u_old:=u_new:
 
 end do:
 
 So,my system:

 u[N/6]=0;
 u[N/3]=0;
 u[N/2]=0;
 u[2*N/3]=0;
 u[5*N/6]=0;
 u[N]=0;

I want to ask advice on how to solve the system.
I wanted to use Newton's method, but I don't know the initial values X_1..X_6.

Tried to set the values X_1..X_6 and to minimize the functional
Fl:=(X_1,X_2,X_3,X_4,X_5,X_6)->(u[N/6])^2+(u[N/3])^2+(u[N/2])^2+(u[2*N/3])^2+(u[5*N/6])^2+(u[N])^2:

with the help with(DirectSearch):
GlobalOptima(Fl);
But I don't know what to do next

Please, advise me how to solve the system! I would be grateful for examples!

Hi there,
I have the following set of equation which I want to solve using Maple's dsolve command:

d[V(t)*C(t)]/dt = G - K *C(t)
dV(t)/dt = alpha - beta

where V is the volume, C is the concentration, and t is the time variable. G, K, alpha and beta are constant parameters of the problem.

The solution for V(t) is easy to find operating the second equation:
V(t) = V_0 + (alpha - beta)*t

But solving for C(t) is a bit harder.

I would like Maple to solve the system, but the result I get does not really make sense to me. My attempt is this: Maple_artificial_kidney.mw

Any thoughts about how to introduce the equations successfully?

Thank you,
jon

Hello everybody,

I am working on a somewhat complicated analytical problem, which includes the evaluation of the following (finite and singularity free) double integral (the stated numerical parameters may vary):

There are no unevaluated variables inside (all the e's are real exp()s). However, when putting it through evalf(), Maple returns it (after calculating for some time) unevaluated. I have no idea why, as a simple box sum should succeed.

Any input would be appreciated. Best regards,

  K. Reuther

P.S.: Raw expression for further deliberations

Int(exp(18.1818*(Int((0.579160e-1*sqrt(x)*Ei(1., 0.500000e-4*x)+(1072.23*(.999950-1.*exp(-0.500000e-4*x)))/sqrt(x))/sqrt(x), x = 1. .. eta))-9.10000*eta)/eta, eta = 1. .. 100.)

is it possible?

i want to prove the injective of new custom operator based on custom logic table according to previous post

and want to see why new custom operator for integers are unordered or sometimes ordered

Hello,

we have a "little" Multibody Dynamics Project in our University. We try to compute a system with 4 DOF with Lagrange. The problem is that in the end our dsolve give us an error. We checked the whole system like 5 times and searched only for the dsolve problem over 6 hours.

Error from dsolve:

"Error, (in DEtools/convertsys) numeric exception: division by zero"

This error shouldnt be possible because we have no divisions at all in our system or somekind of inifity though arctan or what ever.

Any help would be perfekt.

Thanks a lot

Wackeraf

MultibodyDynamics_Gruppe_aktuell_V3.mw