hi. please help me about modeling population migration. i am writing a sample for drawing plot in maple, however maple has not any error, but working incorrectly. my program is:

> u := am+bm*m0*x*t+cm*n0*y*t;
> v := an+bn*m0*x*t+cn*n0*y*t;
> eq1 := diff*(x*t, t) = 2*A*(sinh(u)-x*t*cosh(u));
> eq2 := diff*(y*t, t) = 2*A*(sinh(v)-y*t*cosh(v));
> init := x*0 = X, y*0 = Y;
> sol := dsolve*({eq1, eq2, init}, {x(t), y(t)}, numeric);
> param :=A=1/(2) ,am=0 , bm=1.2*(e)^(-4) ,cm= 0.5*(e)^(-5) ,;
> m0 = 10000, an = 0, bn = -exp(-4), cn = 1.2*exp(-3), n0 = 1000;
> save u, v, param, " narm.sav";
> init := x*0 = 0, y*0 = 0.1e-1;
> eq1 := subs*(param, eq1);
> eq2 := subs*(param, eq2);
> sol := dsolve*({eq1, eq2, init}, {x(t), y(t)}, numeric);
> with*plots;
> odeplot*(sol, [x*t, y*t], 0 .. 60, numpoints = 300, view = [-1 .. 1, -1 .. 1]);
> odeplot*(sol*d, [t, x*t, y*t], 0 .. 60, numpoints = 300, orientation = [70, 55], colour = black, axes = normal);
> init := x*0 = -1, y*0 = 1;
> sold := dsolve*({eq1, eq2, init}, {x(t), y(t)}, numeric);
> odeplot*(sold, [x*t, y*t], 0 .. 60, numpoints = 300, view = [-1 .. 1, -1 .. 1], color = black);

I have two polynomials, say

x(t) = a[0] + a[1]t + a[2]t^2

y(t) = b[0] + b[2]t + b[3]t^2

with the following conditions, x(0) = 2 and y(0) = 1. The polynomial are related by

x'(t) = y(t)

y'(t) = x(t)

and x(0) = 2, y(0) = 1.

One can show that the solution is x(t) = 2 + t + t^2 + O(t^3) and y(t) = 1 + 2t + 0.5t^2 + O(t^3).

I am trying to write something so that it works on a larger system. I have a set of nonlinear DEs which i want to solve this way. Here is a sample code of a much simplifed problem.

See the uploaded file for details.

edit: link is reparied. I uploaded the wrong file initially

edit2: Some improvements were made, but Maple still not returning.

2-2_tests_01.mw

I've got an excel file and a maple file saved in the same folder and am trying to import the contents of the excel file using  the Exceltools[Import]("name of file.extension") command as given in the manual but it doesn't seem to work. I've also tried copying the entire path of the file in the argument but that doesn't seem to work either. What am I doing incorrectly?

Is it possible to use an option similar to range when using lsode method for dsolve? The ODEs I am trying to solve is stiff and will not work with the flag stiff=true or with method=rosenbrock unless i set Digits:=20. I want to avoid doing that as much as possible, since I believe wit will be very taxing, computationally. I have a very large systeom to solve. I found that method=lsode works with the default Digits=15. 

However I need to have the solutions in a given range stored for future access and manipulations. Using range gives me an error: 

Error, (in dsolve/numeric/an_args/lsode) lsode keyword was range, optional keyword must be one of 'ctrl', 'initial', 'itask', 'output', 'procedure', 'procvars', 'start', 'number', 'abserr', 'relerr', 'maxfun', 'minstep', 'maxstep', 'initstep', 'startinit', 'implicit', 'optimize', 'complex'

I cannot figure out how to use range or something similar with lsode. Anyone knows? 

Hello,

I am solving a large system of ODEs, using the following command,

> Sol := dsolve({seq(ode[j], j = 0 .. 21), seq(v[j](0) = 0, j = 1 .. 21), v[0](0) = 1}, [seq(v[j](t), j = 0 .. 21)]);

>

and then plot the quantities I want by something like

> plots:-odeplot(Sol, [t, v[3](t)+v[5](t)], t = 0 .. 1.5);

My problem is that, I do not know a priori which quantity I want to plot, and plotting using above method requires solving the ODEs each time separately, which takes a long time.

So I was curious if there is a scheme that I can solve my system for once and for all, and then plot any quantities that I would like to see.

Hi! I have to solve a very large system of ODEs for a set of functions that can be indexed with two integers, n and l, call them f[n,l](x), where n ranges from 1 to an order 1000 number and l ranges from 1 to an order 10 number. 

What is the best strategy in dsolve, keeping in mind I will only need to have quick access to the values of the first few such functions, say f[0,0](x) and f[0,1](x) for any given x in a given range. I dont need to store the rest, and trying to store all of them actually will give me a warning about length of output exceeds limit of 1000000. So, having maple output the result of dsolve in a procedural form seems to be what I want. However in doing so, for evey value of x i call the solving procedure maple will solve the system each time, thus taking a long time. 

My dilema is: to save storage space and memory I want to use dsolve output as procedures. To save time, when accessing the functions I need (only the first few in the large set of unknown functions that obey the ODE) for a given value of x, I would like maple to have those already stored instead of computing them with each invocation. So, what I need is a somehow split behavior of output. For some of my unknown functions to act as if i use range in dsolve, so that I have access at the values without further computation, and for the rest of the solutions (most of the unknown functions) just keep them in procedural form. How do I achieve this? 

Below is a schematic of my probem:

ODEs:={Set of about 10000 coupled ODEs for functions labeled f[n,l] n ranging from 1..1000 and l from 1..10}

dsolve(ODEs,numeric, output=???)#What options to put in dsolve such that I have quick access to f[0,0](x) and f[0,1](x) in a given range (say xmin=0.1 and xmax=10 for example) but not store the rest. 

guys, is there any possibility to obtain field equations of einstein-hilbert action?

best regards

i want to make a solution that makes sphere not using Sphereplot()

i mean it`s fomula for a sphere that x^2+y^2+z^2=r^2

but, as you know

if i do like this

plot3d(x^2+y^2+z^2,x=0..1,y=0..1)

Warning, expecting only range variables [x, y] in expression x^2+y^2+z^2 to be plotted but found name z

how to make a sphere !?

How to compute Eigenvalues and Eigenvectors of the matrix over finite field? Thank you very much

Dear all,

I want to plot for example cos(theta) from 0 to 2*Pi inside a circle at a radius R. The axis theta of the ploted function is at a radius R.

Is it possible ?

Thanks

Hello,

The idea: parameter "a" will have a new random value each 10 days.

The way I did it is working but it can get very long especially if I do it for a system of equations and for long time more than a year (365 days).

The code:

with(DEtools); with(plots);
n := 5;

for i to n do Ra[i] := RandomTools:-Generate(distribution(Uniform(0.1e-1, .5))); a[[i]] := Ra[i] end do;

b := 0.1e-2;

T := 10;

 eq := diff(L(t), t) = a*L(t)-b;

init[1] := L(0) = 100;
 sol[1] := dsolve({init[1], subs(a = a[[1]], eq)}, L(t), range = 0 .. T, numeric);

init[2] := L(T) = rhs(sol[1](T)[2]);

sol[2] := dsolve({init[2], subs(a = a[[2]], eq)}, L(t), range = T .. 2*T, numeric);

init[3] := L(2*T) = rhs(sol[2](2*T)[2]);
sol[3] := dsolve({init[3], subs(a = a[[3]], eq)}, L(t), range = 2*T .. 3*T, numeric);

p[1] := odeplot(sol[1], [t, L(t)], t = 0 .. T);

p[2] := odeplot(sol[2], [t, L(t)], t = T .. 2*T);

p[3] := odeplot(sol[3], [t, L(t)], t = 2*T .. 3*T);

p := display([p[1], p[2], p[3]]);
display(p);

Thank you

solve.mw

Hi all, I want to determine the roots of this figure that attached here. But the function has 5 parameters so the code doesn't work for it! Help me.

I have another question: The code that attached, determine the roots on horizontal axis, how could I find the values of root on vertical axis?In this figure I want to know the value of F(0) that cut the vertical axis?

Regards

Hello,

I'm trying to calculate and export a huge "Strategy-Matrix". The calculation works in a smaller version with less possible strategies. For this version I use a specialised simulation-PC with enough memory (that was my problem using my own PC). But now I've got this error massage ("Kernel connection has been lost. You should save this worksheet and restart Maple. Executing commands in Maple requires a connection to the Maple kernel. Firewalls have been known to cause problems with kernel connections in Maple. Please ensure that any firewall software is configured to allow Maple to create connections to the kernel. Consult the FAQ for more information.")

The firewall should not be the problem. Below the code. num_Strat is 63001. 

Thank you in advance for your help.

m_Gi:= Matrix(num_Strat, num_Strat):

for t from 1 by 1 to num_Strat
do
Digits:= 5:
v_h2:= evalf((-(24/17)* v_pEi[t]^2)*v_h1 + (v_pEi[t]*(-(128/51)*v_pKi[t]+980/51))*v_h1+ v_pEi[t]*((40/51)*v_pKj+(15/17)*v_pEj)-((110/51)*v_pKi[t]^2)*v_h1 + v_pKi[t]*((55/51)*v_pKj+1220/51*v_h1+(10/17)*v_pEj)):
m_Gi(..,t):= v_h2:
end do:
ExportMatrix(TestMatrixNK, m_Gi):

plot.mw

Hi all. Please help me to plot this. I don't know why it doesn't work correctly.

Regards

Hi, 

     I have a question regarding pdsolve, or Solve from the PDEtools package. I have a set of equations relating partial derivatives, and I'd like to isolate certain terms without explicitly known the functions. I can do this for a single equation, but not multiple ones. I'm curious if Maple can currently handle a system of eqns like these easily, since I will be increasing the number of eqns in the future. Here's the code 

Download PDESolveHelp.mw