## what does 'rsolve', the Maple command do?...

Hi Mapleprimes,

We know that '' rsolve '' is a recurrence equation solver.  It is more than an expression simplifier.

Congratulations to the Maple computer algebra team for creating such a great computer tool.  simply want to know more.

rsolve_on_May_16_2017.pdf

Surely there are many steps to determine the values to place.

Regards,

Matt

## Maple symbolic toolbox for matlab...

I have an important question regarding the symbolic toolbox from maple for matlab

I have declared three symbolic variables:

syms t;

x = sym('x(t)', 'real');

y = sym('y(t)', 'real');

z = sym('z(t)', 'real');

f = x + y + z;

Now i want to  calculate the derivation of f regarding the variable t:

df = diff (f, t);

The result should be

df = dx(t) / dt + dy(t) / dt + dz(t) / dt

but  instead i get the result df = 0. If i use the command diff (f) i get the result

df = dx(t) / dt

Does anyone know how to fix this problem?

## Mean - Real Measurements unsorted...

Hi Guys,

got my last Problem solved :) Now i do have a new one...

restart; with(RealDomain); with(CodeGeneration); with(ExcelTools); with(plots);
dx1:=133;
dy1:=132;

n := 1;
for i to 256 do
for j to 256 do
r := evalf(sqrt((i-dx1)^2+(j-dy1)^2));
Ints := R0[i, j];
IntsR[n] := [r, Ints];
n := n+1
end do;
end do;
IntsR := [seq(IntsR[i], i = 1 .. n-1)]

The list IntsR consists of unsorted values doublets. Now it would be nice to get some sort of mean value of my Ints over r.  Didnt really find a solution for it until now... do you have any hint?

## deleted question?...

I saw a question earlier today regarding 3 particles in a gravitational field.  It has since disappeared.

I see no rhyme or reason why it was deleted by someone, whoever did.  It was Math related and it was Maple related hence no reason for it's deletion.  This happened just before Mapleprimes and Maplesoft went offline for a couple of hours.

Can someone explain?  Did anyone else see it?

## how to calculate probability of matrix in this cas...

i count the number among group
but when the list a large such as over 1000 records, the count will be over 30,000
use which denominator to find probability?
is there any functions in maple for this case?

with(LinearAlgebra):
correlationlist1 := [[1,2,3],[1,3,5]....]:
PAB := Matrix(50):
for ii from 1 to nops(correlationlist) do
for jj from 1 to nops(correlationlist[ii]) do
for kk from 1 to nops(correlationlist) do
for qq from 1 to nops(correlationlist[kk]) do
if ii <> kk then
#print("scan=",correlationlist2[kk],"kk=",kk,"qq=",qq,"row=",correlationlist[ii][jj],"column=",correlationlist[kk][qq]):
PAB[correlationlist[ii][jj],correlationlist[kk][qq]] := PAB[correlationlist[ii][jj],correlationlist[kk][qq]] + 1: # group to group relations
end if:
od:
od:
od:
od:

## How to use statistics correlate function with a li...

If there is a list

[[1,2],[2,2],[3,3]...

how to use correlate function?

assume [1,2] and [2,1] count as 2

when find correlation between 1and 2

if [1,2] and [2,1] are different

first position as key

means 1 relate with 2 count as 1 case for 1 only

[2,1] count as 1 case for 2 only

## Do maple have a Boolean simplify?...

If input a group of outputs which are binary numbers

can it simplify to give a logic that can output these outputs

## Incomplete worksheet...

I'm studying for an exam and earlier today I was working in my Maple sheet. I saved it as usual and no .bak file was in the directory afterwards. When I now later in the day tried to open it, I got the pop up message:
"There were problems during the loading process. your worksheet may be incomplete" and half of the content is missing.

I've searched on google and tried to solve it myself, unfortunately without any luck.

I would really appreciate any help.

Best regards

## Maybe a limitation to DocumentTools :-RunWorkshee...

Hi everybody,

I am used to using  the DocumentTools :-RunWorksheet command and I am confronted with the following problem :

• My worksheet is deemed to generate a plot in a png file
> plotsetup(png, plotoutput=MyFile, …)
> plot( …)

But MyFile is not created

• Thinking that the mw file launched by  DocumentTools :-RunWorksheet is seen as a procedure, I have thought that print(plot(…)) would fix the issue.
Unfortunately the png file is still not created (same behavior with a jpeg file)

I know that  DocumentTools :-RunWorksheet is “experimental” but I wonder if I am doing something wrong or if DocumentTools :-RunWorksheet simply does not support plots ?

## intial value is not converging...

Dear please check once it showing an error program.mw as intial value is not conververging

## Using the Poincare-Dulac Thm to develop an iterati...

proc:=changecoords([x,y],logrithmic);
dsolvediff(x,t):=Matrix([[ 1 , 0],
[ 0 , i ]])x+Vector[column]([ x^i ,
0]);
x:=y+h(y);
for i from 2 to 10;
solve;
end do;

endproc;

MakeIteration := proc( expr::algebraic, x::name )
local iteration := diff(x,t)=Matrix([[ 1 , 0 ],
[ 0 , i]])x+Vector[column]([ x^i ,
0]),x:=y+h(y);
for i from 2 to 10 ;
end do ;
unapply(iteration, x);
endproc;
Error, missing operator or ;

Hi , I would like to solve nonlinear equation by using maple,but I wander how I can write it since there is a change of veriable bysubstiuting the value of y in the diifferenitat . Iactually coudnot write the equation sine it contains differential ,matrix and vector . do Ihave to use dalembertian?

best regards

## How Do I Display Fast Parametric Surface with Smoo...

Dear MaplePrimes team,

For display a beautiful parametric surface with hole (or pierced surface), unfortunately unlike CAD (computer-aided design) tool, Maple requires a large grid number. Therefore the computing is too long and requires huge RAM only with grid of [500,500] on a modern computer.

Fig. 1: For a complex surface such an aircraft fuselage here, Needs high grid, therefore, the compilation is too long with modern computer (Intel Core i7, 2.4 GHz CPU, 16 GB RAM). If displays supplementary windows holes on fuselage, the software may bug when the RAM is full.

Question 1:

In matter of low cost computation, do you have better solution to create a fast surface with smooth hole from given a parametric surface equation S(x,y) and given any convex/concave-closed parametric curve C(t) that is projected on surface S(x,y)?

Fig. 2: These surfaces are been created in Maple 2016 from given initial parametric surface equation and curve projected on the surface. But requires huge grid, huge RAM, and long time to get smooth hole when free CAD tools are fast and low RAM.

Question 2:

Because I noted that the Maple’s view option renderers smooth edges of any surface. My question is: Is it possible to use the internal structure of display/view used by Maple to apply at free boundaries? I think it will be possible to create many class of view, for examples: cylindrical sector, spherical sector views or more rather to be limited with conventional cubic view. And thus, I think it will be possible to create own smooth hole or multiple-hole on surface and get smooth pierced surface.

I know that almost of technicians/engineers will recommend me to use a CAD tool to create surface with hole. But the objective here is to keep purely a mathematical mind and work with exact equations (analytic equations).

On a free CAD tool, even the complex holes are created very fast only with low RAM. Why not on Maple?

Best.

Guy.

## Nonlinear equations ...

Hi , I have question about poincare's theorem which needs to change the veriable and that transform the coordinates , but when I change the variable from x to y if doesn't defined any suggestion ? Please

## How do I solve non linear optimization problem and...

I want to maximize a total profit (TP) function which is dependent on five independent variables (E,W,T, theta, tp). All these five variables can have non negative values. The TP function is given below- ( first TP is directly copied from maple worksheet and than copied again as a picture for clear viewing).

TP = (p1*(Q-q)+p1*(1-theta)*(q-E)+s*E-c*Q-o-h*((1/6)*alpha*W^beta*a*p1^(-b)*tp^3/m-(1/2)*alpha*W^beta*a*p1^(-b)*tp^2+Q*tp)-(t1-tp)*h*((1/2*(-(2/3)*t1+m-(1/3)*tp))*W^beta*a*alpha*(t1-tp)*(-p1*(-1+theta))^(-b)+W*m)/m-(1/2)*h*(W+E)*(T-t1))/T-u*W

I am trying the maximize TP with respect to above five independent variables. I tried to solve  five equations ( representing first order partial derivative of TP with respect to each of the independent variables equated to zero) simultaneously by "solve" and "fsolve" command but both these commands fail to give any output. I have also tried three other commands in optimization package ( QPSolve, NLPSolve, Maximize) but all these three commands also doesn't give any output. I want to prove the concavity of TP function with respect to five independent variables, please guide how it can be done. ( I have computed the Hessian matrix but since five first order equations doesn't give output ( through fsolve command) so I am unable to compute Hessian at these first order optimiality condition solution.). The values of the paramters in the TP equation are -

[alpha = 50, beta = .7, c = 20, h = 4, m = .4, o = 10, p1 = 40, s = 10, u = 5, a = 15000, b = 2]

## How can i get rid out of this error?...

restart;
with(plots); with(DEtools);
&epsilon; := .1;
de1 := x[0](t)+&epsilon;*x[1](t);
ode2 := sin(t)-&epsilon;*t*sin(t);
MODEL := {ode1, ode2};
VARS := {x(t), y(t)};
DOMAIN := t = 0 .. 20;
RANGE := x = -3 .. 3, y = -3 .. 3; COLORS := [BLACK, BLUE];
IC1 := [x[0](0) = 0, x[1](0) = 0]; IC2 := [(D(x[0]))(0) = 1, (D(x[1]))(0) = 0];
DEplot(MODEL, VARS, DOMAIN, RANGE, [IC1, IC2], stepsize = .1, arrows = THIN, linecolor = COLORS);
Error, (in DEtools/DEplot/CheckInitial) the 'number' option must be specified before initial conditions

 3 4 5 6 7 8 9 Last Page 5 of 205
﻿