I want to solve numerically the PDE:

u_xx + u_yy= = u^{1/2}+(u_x)^2/(u)^{3/2}

My assumptions are that  |sqrt(2)u_x/u|<<1 (but I cannot neglect the first term since its in my first order approximation of another PDE.

So I tried solving by using pdsolve in maple, but to no cigar.

Here's the maple file:

 nonlinear.mw

Download nonlinear.mw

I created my own package and one procudure inside this package is to calculate the eigenvalues of a identity matrix, as shown below

Then, I use this package in another document, all other procedures work well except the procedure calculating egigenvalues, as shown below

How to solve this problem? Thank you.

Hi,

Another newbie question...

I have an excel file with lots of named ranges. I can use the ExcelTools WorkbookData command to get a list of all the available named ranges, and then I can assign the contents of each named range to a Maple variable with the Import command.

I'm thinking that there must be an elegant way to create a Maple variable for each string in the list of named ranges and assign it to the contents. Something like looping through each named range in the list, creating a variable of that name, and then assigning the contents of that named range. I hope I'm clear without repeating myself too much.

Am I explaining my need?

Thanks in advance. This looks like a good forum, and I hope to be able to contribute someday.

Mike McDermott

i was to design a GUI using embedded components, i started as shown in the screenshot i uploaded but i had alittle problem trying to code the `evaluate at params` bottton i wanted the botton to return the evaluation of the mathematical expression at the value of the various parameters [delta,gamma,alpha,G,C] into the mathcontainer infront of it i.e evaluating at these points which can be change as new vualues are type in the textfeilds. such as

this is the code i used and it poping out an error message

use DocumentTools in

module()
local f,q,delta,gamma,alpha,G,C;

Do(f=%MathContainer0);
Do(delta=value(%MathContainer1));
Do(gamma=value(%MathContainer2));
Do(alpha=value(%MathContainer3));
Do(G=value(%MathContainer4));
Do(C=value(%MathContainer5));



q:=eval(f,value,delta,gamma,alpha,G,C);
Do(%MathContainer6=q);




end module;
pls i need some one to correct me

FP3.mw

Hi all,

I write this code with some matricial calculations.

There are two important parameters N and Digits (in red in the code).

For N=8, I can do the calculation with Digits=10. For N=10, to get the results (without overflow), I should give Digits=30.

I want to do the calculation for N=100 and higher.

I do not know how to do this. I have a PC with 4G of RAM.

Thank you

While I am working with a problem , suddenly the message 'Kernel connection lost' appeared. After that I restarted maple and then I opened my work sheet again , but I lost all my equations and data in the work sheet. I dont know why? All my work is total empty. 

How can I recover the work sheet data? Solutions are appreciated.

I have attached the lost work sheet for your reference. 

Thanking you.

Hello All,

I looked through the Maple help on PDE systems and pdsolve and the physics problems that appear there. THere are a number of single-PDE cases with initial / boundary conditions; but I couldn't find PDE systems with ics/bcs.

Would you have a (simple) example of a PDE system with its initial / boundary conditions? I am attempting to build understanding of the syntax and different options of "pdsolve". Examples seem to be a great way to learn how to solve PDE systems. One can then pdsolve the PDE system without ics; then add them, try different options etc.

Thank you!

Hello everybody.

In the attached file, you find 6 equations. All of these parameters are known except "pd" and "qd". How can I find these two unknowns from 6 equations??? It should be pointed out that, "pd" and "qd" must contain "ud", "vd","wd" and "rd".

Thanks in advance.

D.mw

I need a help about running the following codes..

restart:
d:=1:
s:=0.8:
c1:=sqrt((s)/Pi/r^3)/2 * (-sin(r*(u-1)^2/4/(s)+Pi/4)+sin(r*(u+1)^2/4/(s))/sqrt(2)):

c22:=-((u+1)/(4*Pi*r*sqrt(y-1)))*(exp(-r*y*(u+1)^2/(8*s)))*(sin(r*y*(u+1)^2/8/(s*(y-1))+Pi/4)/2+sin(r*y*(u-1)^2/8/(s*(y-1))+Pi/4)):
c2:=Int(c22, y= 1..100):

c33:=(sin(r*y*(u+1)^2/8/(s*(y-1))+Pi/4)+sin(r*y*(u-1)^2/8/(s*(y-1))+Pi/4))*(((d*s)/(Pi*(r^2)*y*sqrt(y-1))*exp(-r*y*(u+1)^2/8/s))-((d^2*sqrt((2*s^3)/(Pi*(y-1)*y^3*r^5)))*exp(d*(u+1)+(2*d^2*s)/(y*r))*erfc((r*y*(u+1)+4*d*s)/(2*sqrt(2*s*r*y))))):
c3:=Int(c33, y= 1..100):

g1:=c1+c2+c3:
g2:= unapply(Int(g1, r= 1..infinity, epsilon=1e-4,digits=7),u):
plot(g2, 0..10, numpoints=100);

Thank you!

VIM.3.mwHi,How do I make a graph and error table for exact and by using this method i.e VIM 

In the below file, I have tried to create tree using Graphtheory and export it to file. But exporting works only if I read the procedure twice and that too after running the procedure.

test.mw

testproc.txt

I was solving thos differential equiation using dsolve:

diff(P(t),t)=b*P(t)+m*t

dsolve({diff(P(t),t)=b*P(t)+m*t,P(t0)=P0});

where t0=2000 and P0=100 and m=9/50

Then, I found b setting at t=2001, P=2*P0=200, I applied fsolve to do this and the result was:

b=-1.5979...

To this point, Maple computes correctly. The problem starts when I tried to plot the result in a single 2D-graph (P vs t) using the computed value for b, the expression involved is:

P(t)=.1126460951*t-0.7049523744e-1-1.090463096*10^1390*exp(-1.597924898*t)

I wrote

plot( .1126460951*t-0.7049523744e-1-1.090463096*10^1390*exp(-1.597924898*t) , t=2000..2001);

The result was a blank plot, only the axes appeared even if I used the y=A..B option in the plot command. Is there a way to get the desired plot without having to evaluate the function at several points and then using this set of points to generate the plot: plot({[a,b],[c,d],[f,g]...})? (that worked, of course).

Thanks and regards.

P.S. I have used Maple for over 15 years now and I've plotted expressions even more complicated than this one, I found this very strange. I guess the 10^1390 has something to do this time, but I'm not sure...

I would like to import the code from Pg 90 on in the following document into Maple. Any suggestions for doing so with minimal editing required after pasting? http://mspace.lib.umanitoba.ca/xmlui/bitstream/handle/1993/14907/Huang_Liji.pdf?sequence=1&isAllowed=y

thanks,