Maple 2016 Questions and Posts

These are Posts and Questions associated with the product, Maple 2016

Hi Maple friends,

I am trying to run a script I have written in Maple 2016 directly from my Terminal on my MacAir. I have OSSierra 10.12.1 and Maple 2016.

My script needs to read a file in and then run, whilst writing to some files on the way. I need this, since I want my script to be run on a large computer on my institute, where I have no interface/Maple open and has to access data on the way.

So, I go to Terminal, and stand in the directory of where the file/scrips is (lets call it terminal_test.mpl OR terminal_test.txt) and then write:

./maple terminal_test.mpl

and get back

No such file or directory 

I also get this if I use .txt file AND if I write "Maple 2016" instead of ./maple. I have tried to write 

open "Maple 2016" terminal_test.mpl 

and then Maple opens in the right file, but does not run it. This is just to show, that my terminal can find the program. But this is not what I want. As I said, I would like it to:

start running a script from Terminal - read a file to be used in the script - on the way, write to 6-8 different files, that I can then access whilst the script is running. So my questions are: 

  1. How do I run the script/file in the first place, and am I in the right place?
  2. Where do I place the files to be read during the running of the script (Reading to script works fine already, so it should work also if I can just run the script/file)?
  3. Where do the files I create in my script end, so I can access them?

Thank you so much for any answer - I have not been able to find anything on the net so far that could help me, but I know you guys are really good on this.

Tomas.

 

Hello everyone.

Please I am trying to obtain series expansion of the expression below in u and v up to order 30 but encounter difficulties cum maple is slow to display solution. Can I get help on the code and what to do to optimize the displayed time of maple?

Thank you in anticipation of your quick and positive responses and suggestions.

convert(series(convert(series((y[n]+((-8 h u^2 v^2-4 u^3 sin(u) h+2 sin(2 u) h u^3+2 sin(2 v) h v^3-4 v^3 h sin(v)+2 v^3 h sin(2 u+v)+2 u^3 h sin(u-2 v)+2 u^3 h sin(u+2 v)-2 v^3 h sin(2 u-v)-u^3 h sin(2 u+2 v)-v^3 h sin(2 u+2 v)-u^3 h sin(2 u-2 v)+v^3 h sin(2 u-2 v)+4 h u^3 v^2 sin(2 u)+4 h u^2 v^3 sin(2 v)-4 h u^3 v^2 sin(u-v)+4 h u^2 v^3 sin(u-v)-4 h u^3 v^2 sin(u+v)-4 h u^2 v^3 sin(u+v)+4 h u^2 v^2 cos(u)+4 h u^2 v^2 cos(2 u)+4 h u^2 v^2 cos(2 v)+4 h u^2 v^2 cos(v)-4 h u^3 v cos(2 u-v)-2 h u^2 v^2 cos(2 u-v)+2 h u v^3 cos(2 u-v)+4 h u^3 v cos(2 u+v)-2 h u^2 v^2 cos(2 u+v)-2 h u v^3 cos(2 u+v)+2 h u^3 v cos(u-2 v)-2 h u^2 v^2 cos(u-2 v)-4 h u v^3 cos(u-2 v)-2 h u^3 v cos(u+2 v)-2 h u^2 v^2 cos(u+2 v)+4 h u v^3 cos(u+2 v)+4 h u^3 v cos(u-v)+4 h u v^3 cos(u-v)-4 h u^3 v cos(u+v)-4 h u v^3 cos(u+v)+4 u sin(u) v^2 h-2 sin(2 u) h u v^2-2 sin(2 v) h u^2 v+4 v h sin(v) u^2-2 v u^2 h sin(2 u+v)-2 u v^2 h sin(u-2 v)-2 u v^2 h sin(u+2 v)+2 v u^2 h sin(2 u-v)+v u^2 h sin(2 u+2 v)+u v^2 h sin(2 u+2 v)-v u^2 h sin(2 u-2 v)+u v^2 h sin(2 u-2 v)) f[n])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-8 h u^2 v^2+8 u^3 sin(u) h-4 sin(2 u) h u^3-4 sin(2 v) h v^3+8 v^3 h sin(v)-4 v^3 h sin(2 u+v)-4 u^3 h sin(u-2 v)-4 u^3 h sin(u+2 v)+4 v^3 h sin(2 u-v)+2 u^3 h sin(2 u+2 v)+2 v^3 h sin(2 u+2 v)+2 u^3 h sin(2 u-2 v)-2 v^3 h sin(2 u-2 v)+8 h u^3 v^2 sin(u)+8 h u^2 v^3 sin(v)-4 h u^3 v^2 sin(2 u+v)+4 h u^2 v^3 sin(2 u+v)+4 h u^3 v^2 sin(u-2 v)+4 h u^2 v^3 sin(u-2 v)+4 h u^3 v^2 sin(u+2 v)-4 h u^2 v^3 sin(u+2 v)-4 h u^3 v^2 sin(2 u-v)-4 h u^2 v^3 sin(2 u-v)+8 h u^2 v^2 cos(u)+8 h u^2 v^2 cos(v)+4 h u^3 v cos(2 u-v)-4 h u^2 v^2 cos(2 u-v)-8 h u v^3 cos(2 u-v)-4 h u^3 v cos(2 u+v)-4 h u^2 v^2 cos(2 u+v)+8 h u v^3 cos(2 u+v)-8 h u^3 v cos(u-2 v)-4 h u^2 v^2 cos(u-2 v)+4 h u v^3 cos(u-2 v)+8 h u^3 v cos(u+2 v)-4 h u^2 v^2 cos(u+2 v)-4 h u v^3 cos(u+2 v)-2 h u^3 v cos(2 u+2 v)+4 h u^2 v^2 cos(2 u+2 v)-2 h u v^3 cos(2 u+2 v)+2 h u^3 v cos(2 u-2 v)+4 h u^2 v^2 cos(2 u-2 v)+2 h u v^3 cos(2 u-2 v)-8 u sin(u) v^2 h+4 sin(2 u) h u v^2+4 sin(2 v) h u^2 v-8 v h sin(v) u^2+4 v u^2 h sin(2 u+v)+4 u v^2 h sin(u-2 v)+4 u v^2 h sin(u+2 v)-4 v u^2 h sin(2 u-v)-2 v u^2 h sin(2 u+2 v)-2 u v^2 h sin(2 u+2 v)+2 v u^2 h sin(2 u-2 v)-2 u v^2 h sin(2 u-2 v)) f[n+1])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-8 h u^2 v^2-4 u^3 sin(u) h+2 sin(2 u) h u^3+2 sin(2 v) h v^3-4 v^3 h sin(v)+2 v^3 h sin(2 u+v)+2 u^3 h sin(u-2 v)+2 u^3 h sin(u+2 v)-2 v^3 h sin(2 u-v)-u^3 h sin(2 u+2 v)-v^3 h sin(2 u+2 v)-u^3 h sin(2 u-2 v)+v^3 h sin(2 u-2 v)+4 h u^3 v^2 sin(2 u)+4 h u^2 v^3 sin(2 v)-4 h u^3 v^2 sin(u-v)+4 h u^2 v^3 sin(u-v)-4 h u^3 v^2 sin(u+v)-4 h u^2 v^3 sin(u+v)+4 h u^2 v^2 cos(u)+4 h u^2 v^2 cos(2 u)+4 h u^2 v^2 cos(2 v)+4 h u^2 v^2 cos(v)-4 h u^3 v cos(2 u-v)-2 h u^2 v^2 cos(2 u-v)+2 h u v^3 cos(2 u-v)+4 h u^3 v cos(2 u+v)-2 h u^2 v^2 cos(2 u+v)-2 h u v^3 cos(2 u+v)+2 h u^3 v cos(u-2 v)-2 h u^2 v^2 cos(u-2 v)-4 h u v^3 cos(u-2 v)-2 h u^3 v cos(u+2 v)-2 h u^2 v^2 cos(u+2 v)+4 h u v^3 cos(u+2 v)+4 h u^3 v cos(u-v)+4 h u v^3 cos(u-v)-4 h u^3 v cos(u+v)-4 h u v^3 cos(u+v)+4 u sin(u) v^2 h-2 sin(2 u) h u v^2-2 sin(2 v) h u^2 v+4 v h sin(v) u^2-2 v u^2 h sin(2 u+v)-2 u v^2 h sin(u-2 v)-2 u v^2 h sin(u+2 v)+2 v u^2 h sin(2 u-v)+v u^2 h sin(2 u+2 v)+u v^2 h sin(2 u+2 v)-v u^2 h sin(2 u-2 v)+u v^2 h sin(2 u-2 v)) f[n+2])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((-6 u^2 h^2-6 v^2 h^2-4 cos(2 u) h^2 u^2 v^2-4 cos(2 v) h^2 u^2 v^2+8 v^2 u^2 h^2 cos(u-v)+8 v^2 u^2 h^2 cos(u+v)+8 u sin(u) v^2 h^2+8 sin(2 u) h^2 u v^2+8 sin(2 v) h^2 u^2 v+8 v h^2 sin(v) u^2+4 v u^2 h^2 sin(2 u+v)-4 u v^2 h^2 sin(2 u+v)+8 v u^2 h^2 sin(u-v)-8 sin(u-v) h^2 u v^2-8 v u^2 h^2 sin(u+v)-8 sin(u+v) h^2 u v^2+4 v u^2 h^2 sin(u-2 v)+4 u v^2 h^2 sin(u-2 v)-4 v u^2 h^2 sin(u+2 v)+4 u v^2 h^2 sin(u+2 v)-4 v u^2 h^2 sin(2 u-v)-4 u v^2 h^2 sin(2 u-v)-4 u v h^2 cos(2 u-v)+4 u v h^2 cos(2 u+v)-4 u v h^2 cos(u-2 v)+4 u v h^2 cos(u+2 v)-2 u v h^2 cos(2 u+2 v)+2 u v h^2 cos(2 u-2 v)+8 cos(u-v) h^2 u v-8 cos(u+v) h^2 u v-8 v^2 u^2 h^2+8 h^2 cos(u) u^2-2 cos(2 u) h^2 u^2+6 cos(2 u) h^2 v^2+6 cos(2 v) h^2 u^2-2 cos(2 v) h^2 v^2+8 h^2 cos(v) v^2-4 v^2 h^2 cos(2 u-v)-4 v^2 h^2 cos(2 u+v)-4 u^2 h^2 cos(u-2 v)-4 u^2 h^2 cos(u+2 v)+u^2 h^2 cos(2 u+2 v)+v^2 h^2 cos(2 u+2 v)+u^2 h^2 cos(2 u-2 v)+v^2 h^2 cos(2 u-2 v)) g[n])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)+((6 u^2 h^2+6 v^2 h^2+4 cos(2 u) h^2 u^2 v^2+4 cos(2 v) h^2 u^2 v^2-8 v^2 u^2 h^2 cos(u-v)-8 v^2 u^2 h^2 cos(u+v)-8 u sin(u) v^2 h^2-8 sin(2 u) h^2 u v^2-8 sin(2 v) h^2 u^2 v-8 v h^2 sin(v) u^2-4 v u^2 h^2 sin(2 u+v)+4 u v^2 h^2 sin(2 u+v)-8 v u^2 h^2 sin(u-v)+8 sin(u-v) h^2 u v^2+8 v u^2 h^2 sin(u+v)+8 sin(u+v) h^2 u v^2-4 v u^2 h^2 sin(u-2 v)-4 u v^2 h^2 sin(u-2 v)+4 v u^2 h^2 sin(u+2 v)-4 u v^2 h^2 sin(u+2 v)+4 v u^2 h^2 sin(2 u-v)+4 u v^2 h^2 sin(2 u-v)+4 u v h^2 cos(2 u-v)-4 u v h^2 cos(2 u+v)+4 u v h^2 cos(u-2 v)-4 u v h^2 cos(u+2 v)+2 u v h^2 cos(2 u+2 v)-2 u v h^2 cos(2 u-2 v)-8 cos(u-v) h^2 u v+8 cos(u+v) h^2 u v+8 v^2 u^2 h^2-8 h^2 cos(u) u^2+2 cos(2 u) h^2 u^2-6 cos(2 u) h^2 v^2-6 cos(2 v) h^2 u^2+2 cos(2 v) h^2 v^2-8 h^2 cos(v) v^2+4 v^2 h^2 cos(2 u-v)+4 v^2 h^2 cos(2 u+v)+4 u^2 h^2 cos(u-2 v)+4 u^2 h^2 cos(u+2 v)-u^2 h^2 cos(2 u+2 v)-v^2 h^2 cos(2 u+2 v)-u^2 h^2 cos(2 u-2 v)-v^2 h^2 cos(2 u-2 v)) g[n+2])/(-12 u^2 v^2+4 sin(u) u^3 v^2+4 sin(2 u) u^3 v^2+4 sin(2 v) u^2 v^3+4 sin(v) u^2 v^3-2 sin(2 u+v) u^3 v^2+2 sin(2 u+v) u^2 v^3-4 sin(u-v) u^3 v^2+4 sin(u-v) u^2 v^3-4 sin(u+v) u^3 v^2-4 sin(u+v) u^2 v^3+2 sin(u-2 v) u^3 v^2+2 sin(u-2 v) u^2 v^3+2 sin(u+2 v) u^3 v^2-2 sin(u+2 v) u^2 v^3-2 sin(2 u-v) u^3 v^2-2 sin(2 u-v) u^2 v^3+8 cos(u) u^2 v^2+4 cos(2 u) u^2 v^2+4 cos(2 v) u^2 v^2+8 cos(v) u^2 v^2-2 cos(2 u-v) u^3 v-4 cos(2 u-v) u^2 v^2-2 cos(2 u-v) u v^3+2 cos(2 u+v) u^3 v-4 cos(2 u+v) u^2 v^2+2 cos(2 u+v) u v^3-2 cos(u-2 v) u^3 v-4 cos(u-2 v) u^2 v^2-2 cos(u-2 v) u v^3+2 cos(u+2 v) u^3 v-4 cos(u+2 v) u^2 v^2+2 cos(u+2 v) u v^3-cos(2 u+2 v) u^3 v+2 cos(2 u+2 v) u^2 v^2-cos(2 u+2 v) u v^3+cos(2 u-2 v) u^3 v+2 cos(2 u-2 v) u^2 v^2+cos(2 u-2 v) u v^3+4 cos(u-v) u^3 v+4 cos(u-v) u v^3-4 cos(u+v) u^3 v-4 cos(u+v) u v^3)),u=0,32),polynom),v=0,32),polynom);

Hi, 

Can we tell maple to give output in a form which contain any specific term or expression. For example- Can I tell maple to give output which contain (x^2+2*x+3) expression in the output. I used various combinations in collect command but it doesn't work. It is shown in enclosed maple worksheet.

I have some questions about output produced by maple. I have attached a maple worksheet in which I asked the questions in comments. Please help mein in this problem because it become very tedious to type each expression by hand in desired form. If maple produce it by some means it would be great help.

desired_form_in_output.mw

Thanks and Regards,

Nilesh

 

 

Dear Community,

Entering a long line in 1D mode breaks at the right window bounday. How could I prevent this? I would like the entered line stay just as one, single, long line extending over the window boundary and I would like to use the bottom window scroll bar to view it. BTW it also does not appear, so how could I invoke it? I use Maple2016. I have neither found any information about this in the help, nor in the documentation.

Also is there a possibility to split the IDE screen vertically, as e.g. in Visual Studio? Looking at long programs this feature is very useful.

tx in advance,

best regards

Andras

I am computing some finite points in the plane. Then I'm using point style plot to plot them. My points have order and I want to conncet the i-th point to the (i+1)-th point with a line segment. I searched the help but I didn't figure out how to do it. Is it possible to do this in Maple?

Is there any method or a command to use or type so Maple start running something or evaluating and if it took more than a specific amount of time, then itself automatically stops itself without needing me to interrupt the evaluation?

For example let's say I wrote a proc, for some inputs it may compute the output fast but in some it may even don't answer after hours. Now I am writing a for and I don't know for which steps it can answer in less than a minute. So I like to add something so that Maple starts the FOR but in steps that it takes more than a minute, it automatically stops evaluation at that step and jump to the next step.

I have a set of points. I want to make a linear spline of that points and plot the resulting function:

function:=CurveFitting[Spline](points_x,points_y,x,degree=1):

However I need the function to define constant y endpoint values where x is out of points_x range (i.e. first-derivation is zero where x is out of points_x range). Endpoint options does not work for me.

             constant1   x < xmin
function :=  ..
             constant2   otherwise

How can I define spline function like this?

I'm trying to obtain integral of Planck radiation law in Maple:

with this command:

f := (2*h*(c)^(2))/((x)^(5))*(1)/(exp((h*c)/(x*k*T))-1);
int(f,x=0..infinity);

but I get some terrible limit that cannot be solved instead of the correct result:

How to obtain correct integral?

planck.mw

I am coding a big module to solving my project : analyze function in math, but when I compile my module maple return "Error," but it doesn't tell me what error happened.

I check the maple help and it said:" If no msgString is given, error raises the most recently occurring exception" but I have no exception before. 

This is the pic of that error.

Thank for your help.

Hi, my dear friend,

i am solving 9 ODE with boundary conditionsNigam.mw

Eq1 := 2.*F1(eta)+diff(H1(eta), eta) = 0

2.*F1(eta)+diff(H1(eta), eta) = 0

(1)

Eq2 := F1(eta)^2-G1(eta)^2+(diff(F1(eta), eta))*H1(eta)-(diff(F1(eta), eta, eta)) = 0

F1(eta)^2-G1(eta)^2+(diff(F1(eta), eta))*H1(eta)-(diff(diff(F1(eta), eta), eta)) = 0

(2)

Eq3 := 2*F1(eta)*G1(eta)+H1(eta)*(diff(G1(eta), eta))-(diff(G1(eta), eta, eta)) = 0

2*F1(eta)*G1(eta)+H1(eta)*(diff(G1(eta), eta))-(diff(diff(G1(eta), eta), eta)) = 0

(3)

Eq4 := 4*F1(eta)*F3(eta)+H3(eta)*(diff(F1(eta), eta))+H1(eta)*(diff(F3(eta), eta))-2*G1(eta)*G3(eta)-2.*F1(eta)^2-1.5*H1(eta)-(diff(F3(eta), eta, eta)) = 0

4*F1(eta)*F3(eta)+H3(eta)*(diff(F1(eta), eta))+H1(eta)*(diff(F3(eta), eta))-2*G1(eta)*G3(eta)-2.*F1(eta)^2-1.5*H1(eta)-(diff(diff(F3(eta), eta), eta)) = 0

(4)

Eq5 := 2*F3(eta)*G1(eta)+4*F1(eta)*G3(eta)+H3(eta)*(diff(G1(eta), eta))-H1(eta)*(diff(G3(eta), eta))-2*F1(eta)*G1(eta)-1.5*H1(eta)*(diff(G1(eta), eta))-(diff(G3(eta), eta, eta)) = 0

2*F3(eta)*G1(eta)+4*F1(eta)*G3(eta)+H3(eta)*(diff(G1(eta), eta))-H1(eta)*(diff(G3(eta), eta))-2*F1(eta)*G1(eta)-1.5*H1(eta)*(diff(G1(eta), eta))-(diff(diff(G3(eta), eta), eta)) = 0

(5)

Eq6 := 4.*F3(eta)+diff(H3(eta), eta) = 0

4.*F3(eta)+diff(H3(eta), eta) = 0

(6)

Eq7 := 6*F1(eta)*F5(eta)-6*F1(eta)*F3(eta)+3.*F3(eta)^2+H1(eta)*(diff(F5(eta), eta))+H3(eta)*(diff(F3(eta), eta))+H5(eta)*(diff(F1(eta), eta))-1.5*(H1(eta)*(diff(F3(eta), eta))+H3(eta)*(diff(F1(eta), eta)))-G3(eta)^2-2*G1(eta)*G5(eta)-(diff(F5(eta), eta, eta)) = 0

6*F1(eta)*F5(eta)-6*F1(eta)*F3(eta)+3.*F3(eta)^2+H1(eta)*(diff(F5(eta), eta))+H3(eta)*(diff(F3(eta), eta))+H5(eta)*(diff(F1(eta), eta))-1.5*H1(eta)*(diff(F3(eta), eta))-1.5*H3(eta)*(diff(F1(eta), eta))-G3(eta)^2-2*G1(eta)*G5(eta)-(diff(diff(F5(eta), eta), eta)) = 0

(7)

Eq8 := 6*G5(eta)*F1(eta)+2*G1(eta)*F5(eta)+4*G3(eta)*F3(eta)-4*F1(eta)*G3(eta)-2*F3(eta)*G1(eta)+H1(eta)*(diff(G5(eta), eta))-1.5*(H1(eta)*(diff(G3(eta), eta))+H3(eta)*(diff(G1(eta), eta)))+H3(eta)*(diff(G3(eta), eta))+H5(eta)*(diff(G1(eta), eta))-(diff(G5(eta), eta, eta)) = 0

6*G5(eta)*F1(eta)+2*G1(eta)*F5(eta)+4*G3(eta)*F3(eta)-4*F1(eta)*G3(eta)-2*F3(eta)*G1(eta)+H1(eta)*(diff(G5(eta), eta))-1.5*H1(eta)*(diff(G3(eta), eta))-1.5*H3(eta)*(diff(G1(eta), eta))+H3(eta)*(diff(G3(eta), eta))+H5(eta)*(diff(G1(eta), eta))-(diff(diff(G5(eta), eta), eta)) = 0

(8)

Eq9 := 6.*F5(eta)+F3(eta)+diff(H5(eta), eta) = 0

6.*F5(eta)+F3(eta)+diff(H5(eta), eta) = 0

(9)

bcs1 := F1(0) = 0, F3(0) = 0, F5(0) = 0

F1(0) = 0, F3(0) = 0, F5(0) = 0

(10)

bcs2 := G1(0) = 1, G3(0) = 0, G5(0) = 0

G1(0) = 1, G3(0) = 0, G5(0) = 0

(11)

bcs3 := H1(0) = 0, H3(0) = 0, H5(0) = 0

H1(0) = 0, H3(0) = 0, H5(0) = 0

(12)

bcs4 := F1(10) = 0, F3(10) = 0, F5(10) = 0

F1(10) = 0, F3(10) = 0, F5(10) = 0

(13)

bcs5 := G1(10) = 0, G3(10) = 0, G5(10) = 0

G1(10) = 0, G3(10) = 0, G5(10) = 0

(14)

R := dsolve(eval({Eq1, Eq2, Eq3, Eq4, Eq5, Eq6, Eq7, Eq8, Eq9, bcs1, bcs2, bcs3, bcs4, bcs5}), [F1(eta), F3(eta), F5(eta), G1(eta), G3(eta), G5(eta), H1(eta), H3(eta), H5(eta)], numeric, output = listprocedure)

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

 

``


Maple Worksheet - Error

Failed to load the worksheet 

Download Nigam.mwNigam.mw

then i got this error

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

i dont know where i need to change.. could you help me..

Dear Maple users,

I have very interesting problem with evaluating of only symbolic equations.

The problem is: (In document mode)
When I try to evaluate all numeric values for example typing 1+2 and press Enter it successfuly evaluates 3.
But when I try to evaluate symbolic values for example x + y nothing happens.

I tried lots of things to solve this, but no luck.
Shortly if the equation contains only numbers it evaluates successfuly.
And if the equation contains one or more symbolic variables like (x, y, z, variable1, test1), it does nothing.

What could be the problem?
 

(Because of this problem I cannot use the document mode, so I'm using the worksheet mode. Worksheet mode works good, but sometimes it is not calculating like document mode)

 

 

System is Windows 7 x64
Maple 2016
8 x 3.50GHz Xeon CPU
128 GB RAM
Windows Language is Turkish, (I tried with also USA English, but no luck either)
Keyboard is Turkish TR, (I tried with also US English, but no luck either)

Hello everyone.

Please can I meet with Computational or/and Numerical anlysts that have worked or working on the algorihms particularly (Runge Kutta Nystrom, Block multistep methods including hybrid and Block Boundaru Value methods) for the solution of both IVP and BVP.

I will appreciante if I can learn from them and possibly collaborate with them. Thank you in anticipation of your positive response.

Find and classify all the critical point(s) for 

f(x,y)=x3 + 3xy2 - 15x + y3 - 15y .

I need a help from someone who knows the GRTensor commands.

If you create a tensor of rank 2 using the grdef ("F {(a) (b)}"). Until then, okay!

However, I would like to define each of the 16 tensor components and I do not know how to do.

Could someone help me?

Thank you very much!

Please I am having problem with this code particularly the last subroutine

#subroutine 1

restart;
Digits:=30:

f:=proc(n)
    -25*y[n]+12*cos(x[n]):
end proc:

#subroutine 2

e1:=y[n+4] = -y[n]+2*y[n+2]+((1/15)*h^2+(2/945)*h^2*u^2+(1/56700)*h^2*u^4-(1/415800)*h^2*u^6-(167/833976000)*h^2*u^8-(2633/245188944000)*h^2*u^10-(2671/5557616064000)*h^2*u^12-(257857/13304932857216000)*h^2*u^14-(3073333/4215002729166028800)*h^2*u^16)*f(n)+((16/15)*h^2-(8/945)*h^2*u^2-(1/14175)*h^2*u^4+(1/103950)*h^2*u^6+(167/208494000)*h^2*u^8+(2633/61297236000)*h^2*u^10+(2671/1389404016000)*h^2*u^12+(257857/3326233214304000)*h^2*u^14+(3073333/1053750682291507200)*h^2*u^16)*f(n+1)+((26/15)*h^2+(4/315)*h^2*u^2+(1/9450)*h^2*u^4-(1/69300)*h^2*u^6-(167/138996000)*h^2*u^8-(2633/40864824000)*h^2*u^10-(2671/926269344000)*h^2*u^12-(257857/2217488809536000)*h^2*u^14-(3073333/702500454861004800)*h^2*u^16)*f(n+2)+((16/15)*h^2-(8/945)*h^2*u^2-(1/14175)*h^2*u^4+(1/103950)*h^2*u^6+(167/208494000)*h^2*u^8+(2633/61297236000)*h^2*u^10+(2671/1389404016000)*h^2*u^12+(257857/3326233214304000)*h^2*u^14+(3073333/1053750682291507200)*h^2*u^16)*f(n+3)+((1/15)*h^2+(2/945)*h^2*u^2+(1/56700)*h^2*u^4-(1/415800)*h^2*u^6-(167/833976000)*h^2*u^8-(2633/245188944000)*h^2*u^10-(2671/5557616064000)*h^2*u^12-(257857/13304932857216000)*h^2*u^14-(3073333/4215002729166028800)*h^2*u^16)*f(n+4):

e2:=y[n+3] = -y[n+1]+2*y[n+2]+(-(1/240)*h^2-(31/60480)*h^2*u^2-(67/1814400)*h^2*u^4-(109/53222400)*h^2*u^6-(18127/186810624000)*h^2*u^8-(64931/15692092416000)*h^2*u^10-(9701/59281238016000)*h^2*u^12-(20832397/3406062811447296000)*h^2*u^14-(11349439/51876956666658816000)*h^2*u^16)*f(n)+((1/10)*h^2+(31/15120)*h^2*u^2+(67/453600)*h^2*u^4+(109/13305600)*h^2*u^6+(18127/46702656000)*h^2*u^8+(64931/3923023104000)*h^2*u^10+(9701/14820309504000)*h^2*u^12+(20832397/851515702861824000)*h^2*u^14+(11349439/12969239166664704000)*h^2*u^16)*f(n+1)+((97/120)*h^2-(31/10080)*h^2*u^2-(67/302400)*h^2*u^4-(109/8870400)*h^2*u^6-(18127/31135104000)*h^2*u^8-(64931/2615348736000)*h^2*u^10-(9701/9880206336000)*h^2*u^12-(20832397/567677135241216000)*h^2*u^14-(11349439/8646159444443136000)*h^2*u^16)*f(n+2)+((1/10)*h^2+(31/15120)*h^2*u^2+(67/453600)*h^2*u^4+(109/13305600)*h^2*u^6+(18127/46702656000)*h^2*u^8+(64931/3923023104000)*h^2*u^10+(9701/14820309504000)*h^2*u^12+(20832397/851515702861824000)*h^2*u^14+(11349439/12969239166664704000)*h^2*u^16)*f(n+3)+(-(1/240)*h^2-(31/60480)*h^2*u^2-(67/1814400)*h^2*u^4-(109/53222400)*h^2*u^6-(18127/186810624000)*h^2*u^8-(64931/15692092416000)*h^2*u^10-(9701/59281238016000)*h^2*u^12-(20832397/3406062811447296000)*h^2*u^14-(11349439/51876956666658816000)*h^2*u^16)*f(n+4):

e3:=h*delta[n] = (-149/42-(16/245)*u^2-(1324/169785)*u^4-(559246/695269575)*u^6-(14310311/175207932900)*u^8-(170868550903/20641246574949000)*u^10)*y[n]+(128/21+(32/245)*u^2+(2648/169785)*u^4+(1118492/695269575)*u^6+(14310311/87603966450)*u^8+(170868550903/10320623287474500)*u^10)*y[n+1]+(-107/42-(16/245)*u^2-(1324/169785)*u^4-(559246/695269575)*u^6-(14310311/175207932900)*u^8-(170868550903/20641246574949000)*u^10)*y[n+2]+(-(67/1260)*h^2+(1241/198450)*h^2*u^2+(277961/366735600)*h^2*u^4+(26460409/333729396000)*h^2*u^6+(1363374533/168199615584000)*h^2*u^8+(16323847966961/19815596711951040000)*h^2*u^10)*f(n)+((188/105)*h^2+(5078/99225)*h^2*u^2+(556159/91683900)*h^2*u^4+(51834031/83432349000)*h^2*u^6+(67782373/1078202664000)*h^2*u^8+(1854079193287/291405833999280000)*h^2*u^10)*f(n+1)+((31/90)*h^2+(341/33075)*h^2*u^2+(79361/61122600)*h^2*u^4+(23456627/166864698000)*h^2*u^6+(1228061399/84099807792000)*h^2*u^8+(14797833720283/9907798355975520000)*h^2*u^10)*f(n+2)+(-(4/105)*h^2-(46/14175)*h^2*u^2-(809/1871100)*h^2*u^4-(27827/567567000)*h^2*u^6-(637171/122594472000)*h^2*u^8-(33500737/62523180720000)*h^2*u^10)*f(n+3)+((1/252)*h^2+(23/28350)*h^2*u^2+(809/7484400)*h^2*u^4+(27827/2270268000)*h^2*u^6+(637171/490377888000)*h^2*u^8+(33500737/250092722880000)*h^2*u^10)*f(n+4):

e4:=y[3] = -y[1]+2*y[2]+(-(1/240)*h^2-(31/60480)*h^2*u^2-(67/1814400)*h^2*u^4-(109/53222400)*h^2*u^6-(18127/186810624000)*h^2*u^8-(64931/15692092416000)*h^2*u^10-(9701/59281238016000)*h^2*u^12-(20832397/3406062811447296000)*h^2*u^14-(11349439/51876956666658816000)*h^2*u^16)*f(0)+((1/10)*h^2+(31/15120)*h^2*u^2+(67/453600)*h^2*u^4+(109/13305600)*h^2*u^6+(18127/46702656000)*h^2*u^8+(64931/3923023104000)*h^2*u^10+(9701/14820309504000)*h^2*u^12+(20832397/851515702861824000)*h^2*u^14+(11349439/12969239166664704000)*h^2*u^16)*f(1)+((97/120)*h^2-(31/10080)*h^2*u^2-(67/302400)*h^2*u^4-(109/8870400)*h^2*u^6-(18127/31135104000)*h^2*u^8-(64931/2615348736000)*h^2*u^10-(9701/9880206336000)*h^2*u^12-(20832397/567677135241216000)*h^2*u^14-(11349439/8646159444443136000)*h^2*u^16)*f(2)+((1/10)*h^2+(31/15120)*h^2*u^2+(67/453600)*h^2*u^4+(109/13305600)*h^2*u^6+(18127/46702656000)*h^2*u^8+(64931/3923023104000)*h^2*u^10+(9701/14820309504000)*h^2*u^12+(20832397/851515702861824000)*h^2*u^14+(11349439/12969239166664704000)*h^2*u^16)*f(3)+(-(1/240)*h^2-(31/60480)*h^2*u^2-(67/1814400)*h^2*u^4-(109/53222400)*h^2*u^6-(18127/186810624000)*h^2*u^8-(64931/15692092416000)*h^2*u^10-(9701/59281238016000)*h^2*u^12-(20832397/3406062811447296000)*h^2*u^14-(11349439/51876956666658816000)*h^2*u^16)*f(4):

e5:=h*delta[0] = (-149/42-(16/245)*u^2-(1324/169785)*u^4-(559246/695269575)*u^6-(14310311/175207932900)*u^8-(170868550903/20641246574949000)*u^10)*y[0]+(128/21+(32/245)*u^2+(2648/169785)*u^4+(1118492/695269575)*u^6+(14310311/87603966450)*u^8+(170868550903/10320623287474500)*u^10)*y[1]+(-107/42-(16/245)*u^2-(1324/169785)*u^4-(559246/695269575)*u^6-(14310311/175207932900)*u^8-(170868550903/20641246574949000)*u^10)*y[2]+(-(67/1260)*h^2+(1241/198450)*h^2*u^2+(277961/366735600)*h^2*u^4+(26460409/333729396000)*h^2*u^6+(1363374533/168199615584000)*h^2*u^8+(16323847966961/19815596711951040000)*h^2*u^10)*f(0)+((188/105)*h^2+(5078/99225)*h^2*u^2+(556159/91683900)*h^2*u^4+(51834031/83432349000)*h^2*u^6+(67782373/1078202664000)*h^2*u^8+(1854079193287/291405833999280000)*h^2*u^10)*f(1)+((31/90)*h^2+(341/33075)*h^2*u^2+(79361/61122600)*h^2*u^4+(23456627/166864698000)*h^2*u^6+(1228061399/84099807792000)*h^2*u^8+(14797833720283/9907798355975520000)*h^2*u^10)*f(2)+(-(4/105)*h^2-(46/14175)*h^2*u^2-(809/1871100)*h^2*u^4-(27827/567567000)*h^2*u^6-(637171/122594472000)*h^2*u^8-(33500737/62523180720000)*h^2*u^10)*f(3)+((1/252)*h^2+(23/28350)*h^2*u^2+(809/7484400)*h^2*u^4+(27827/2270268000)*h^2*u^6+(637171/490377888000)*h^2*u^8+(33500737/250092722880000)*h^2*u^10)*f(4):

 

#subroutine 3

 


inx:=0:
ind:=0:
iny:=1:
h:=Pi/4.0:
n:=0:
omega:=5:
u:=omega*h:
N:=solve(h*p = 500*Pi/2, p):

c:=1:
for j from 0 to 5 do
    t[j]:=inx+j*h:
end do:
#e||(1..6);
vars:=y[n+1],y[n+2],y[n+3],delta[n],y[n+4]:

printf("%6s%15s%15s%15s\n",
    "h","Num.y","Ex.y","Error y");
for k from 1 to N do

    par1:=x[0]=t[0],x[1]=t[1],x[2]=t[2],x[3]=t[3],x[4]=t[4],x[5]=t[5]:
    par2:=y[n]=iny,delta[n]=ind:
   

res:=eval(<vars>, fsolve(eval({e||(1..5)},[par1,par2]), {vars}));

    for i from 1 to 5 do
        exy:=eval(0.5*cos(5*c*h)+0.5*cos(c*h)):
        printf("%6.5f%17.9f%15.9f%13.5g\n",
        h*c,res[i],exy,abs(res[i]-exy)):
        
        c:=c+1:
    end do:
    iny:=res[5]:
    inx:=t[5]:
    for j from 0 to 5 do
        t[j]:=inx + j*h:
    end do:
end do:

 

 

 

 

1 2 3 4 5 6 7 Last Page 2 of 38