Hi,

I cannot thank you enough for the help.

I figure that there is a way to write a bunch of lines on a DE plot using some sort of command? I figure that it has to do with the sequence feature of maple.  What I am trying to do is show a bunch of lines in a DE plot that does y(n)=-1 where n is an integer ranging from -5 to 5.  I'll link my worksheet. Differential_Equation_slope_field.mw

  Hello,

     I have a large module of procedures. A procedure thet I know takes a long time to execute the result. I would like to know at what stage of the procedure the calculation takes time. 

Here is the code I've wrote for third order polynomial .

I want to extend it for 100th order polynomial without writing the somehow same things 100 times !

what should I do ?

( I've done same thing in PHP with sprintf , there we could define variables like x{i} = ... , and so we could make a lot of variable with a for loop )

Thanks in advance .

Here is the code in *.mw format :

Test1.mw
 

with(plottools):with(plots): display(seq(seq(display(polygon([[i,j],[i,j+1],[i+1,j+1],[i+1,j]], color=`if``((j)::odd,ColorTools:-Color=magenta))), textplot([1+.5,j+.5,fprintf("%d",i*j)])),i=1..10), j=1..10),axes=none);

Hello Maple experts:

Using Maple 2018.2.1 on windows, with Physics package version 301.

There is some strange issue in parsing somewhere.  This works

r:='r'; u:='u'; t:='t';theta:='theta';k:='k';
k:=1/25;
pde := diff(u(r,theta,t),t)=k*(diff(u(r,theta,t),r$2) + 1/r*diff(u(r,theta,t),r)+1/r^2*diff(u(r,theta,t),theta$2));
bc1:= eval(diff(u(r,theta,t),r),r=1)=0;
bc2:= u(r,0,t)=0, u(r,Pi,t)=0;
ic := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
sol:=pdsolve([pde, bc1, bc2, ic], u(r, theta, t),HINT = boundedseries(r = [0]));

No problem with the above, pdsolve gives solution. Now if I change k from 1/25 to 1, it gives parsing error in map

k:=1;
pde := diff(u(r,theta,t),t)=k*(diff(u(r,theta,t),r$2) + 1/r*diff(u(r,theta,t),r)+1/r^2*diff(u(r,theta,t),theta$2));
bc1:= eval(diff(u(r,theta,t),r),r=1)=0; 
bc2:= u(r,0,t)=0, u(r,Pi,t)=0;
ic  := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
sol := pdsolve([pde, bc1, bc2, ic], u(r, theta, t),HINT = boundedseries(r = [0]));

Error, (in assuming) invalid input: map expects 2 or more arguments, but received 1

any K value less than one works. 

Now I change K from 1 to 2, the error goes away. But no solution. Which is OK.

Why does the map error shows up when k=1 only? Same error shows up if K is missing all togother, as in

r:='r'; u:='u'; t:='t';theta:='theta';k:='k';
pde := diff(u(r,theta,t),t)=diff(u(r,theta,t),r$2) + 1/r*diff(u(r,theta,t),r)+1/r^2*diff(u(r,theta,t),theta$2);
bc1:= eval(diff(u(r,theta,t),r),r=1)=0; 
bc2:= u(r,0,t)=0, u(r,Pi,t)=0;
ic  := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
sol := pdsolve([pde, bc1, bc2, ic], u(r, theta, t),HINT = boundedseries(r = [0]));

Error, (in assuming) invalid input: map expects 2 or more arguments, but received 1

Hi

I want to create a Heaviside function called MyStep() which will behave exactly like the built in Heaviside with one exception: for MyStep(0) it will return 1 instead of undefined.

MyStep(t) = 0, 0<t

MyStep(t) = 1, else

Thanks

Hi

I got an square matrix (70×70) from MATLAB (please download attached text file 1.txt). Following codes are used in MAPLE, but an unknown error is occurred. It seems that matrix is divided in two submatrices. Please hint me or run the codes to obtain fiirst three minimum real positive roots of the determinant. 

Thank you for taking your time

Maple codes

Digits:=150:

M:=parse(FileTools[Text][ReadFile]("1.txt")):

LinearAlgebra:-Determinant(M):

fsolve(%,P=0..0.5)*100;

Good morning everyone

I'm sorry if I'm asking a simple question, but I would like to create a maplet from a maple code I written.

My first problem is I would like to now how to create a table like in excel, to be fill up by the user of the maplet.

Basically, the maplet will consist of enter two parameters who define the number of colums and row needed for the table,

Then fill up this table, which gonna be plot, and then a button to fit this curve using an equation. And the starting value in order to fit the curve will also be set up by the user.

Thanks you in advance if you can help me to set up this table in the maplet.

cheers!

Hello,

in Maple 2017 (could not try 2018) the following is a bug;
 

with(LinearAlgebra);
U:= Matrix(4,4, shape=triangular[upper, unit], symbol =u);
-U;

It returns the wrong matrix, with +1 on the diagonal and -u_{ij} on the upper part. That is, it changes the sign only of the entries in the upper part. The bug is for any size, and also for lower triangular.

 

I have been trying to use fsolve to iteratively solve a free convection heat transfer problem for a heated vertical plate using units and the ThermoPhysicalData library.  I have received some help and have managed to get rid of the eror messages but I still can't seem to get it to solve.

Here is the code

I would appreciate help wih this.

Thanks

For example if a:=[3,3,1,5,7,8,5,4,4,4,4,3,9]

What is the simplest way, or alternate ways, to remove only one 4 for example from the dataset.

model.mw

after loop how to pick equations automatically solve system of equations??? 

Hi

I've just upgraded from v2017.3 to 2018. It worked OK until I installed latest service pack 2018.2.1 (server license provided by my University). Ever since I cannot use Maple. Anything I type I get Typesetting:-mparsed(...) error and the text/command I typed.

I’ve contacted our software tech support and they told me to change the typesetting level from advanced to standard and it did fix the problem.

But why does it happen in the first place? I’m running Windows 8.1 64bit. Out tech support told me to has something to do with 3D display issue on my machine and told me to bring my laptop on Monday to see if they can resolve the issue.

Anyone else have this problem? Why didn’t it happen with older Maple versions? What am I missing by using standard typesetting instead of the default advanced?

Thanks

how can I get a positive output in the following example

solve((hv*hw+lv*lw)/(lv+hv)*(lv+hv) = dw*(lv+hv), hv)

instead of -lv*(dw-lw)/(dw-hw) I would like lv*(dw-lw)/(hw-dw)