I am plotting a pair of lists of points using

pointplot(Listap, Listbp, symbol = point, symbolsize = 1, size = [1200, 1200])

How could is do this with plot so I can add colours? Along the lines Listap(i)^2+Listbp(i)^2 =R, R is in the range 0..1,then colour =R*256 or any other imaginative way of adding colour.
 

Download 096-Chayley_transform_for_MP_question.mw

i tried to solve a nonlinear ode with numerical method but maple can't solve it and this error occur:

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

my maple codes are attached below:

numeriacal_sol.mw

can any help me?

As I understand it, Maple will detect and use the available cores in a system, if the calculation is suitable for multi-core use.

As I am installing Maple on a multi-user cluster, using a scheduler to run maple scripts, I want to ensure the maple jobs only use the number of cores allocated to the job.  

Is it possible to set the number of cores used ? 

If I have misunderstood how Maple works (I am new to it), or if there is a section in the documentation which explains this, please point me in the right direction.  I haven't found this info so far.

if I use the command below, I will get the right answer.

and the answer is :

while I use the same command for the compliacated Ei function, it doesn't work any more, I don't know the reason.

 and it's reault is :

Dear Maple users

I was using the animation command to create an animation consisting of a running sine curve with a ball following the curve, as shown below. I wonder if it is possible to create it with the Explore command too? Or is this command more limited? I need to be able to add more animations, like I can with the display command ...

Regards, Erik.

Hello people in mapleprimes,

I want to modify the inequality of e assuming d.
But, e/(1-a) assuming(d) does not work well.
How can I do for this?

d:=(a>0,a<1,b>0,c>0);
                   0 < a, a < 1, 0 < b, 0 < c
e:=(1-a)*b>c;
                         c < (1 - a) b
e/(1-a)assuming(d);
                         `*`(c<(1-a)b,1/(1-a))

Thanks in advance.

inequality.mw

I have made an animate-plot with trace function enabled and I was wondering if there was something I could do so it would show the value of the variable when trace reached frame x and y and plot it in the plot or maybe somewhere else.

Thanks in advance.

Hello,

when trying to solve my equation I always get this error-message:

Error, (in DEtools/convertsys) invalid input: degree expects 1 or 2 arguments, but received 3

To solve, I use this one:

Anfangsbed := phi(0) = 0,  D(phi)(0) = 10;

Einst :=             range = 0..10,            #Integrationsgrenzen
                     relerr = 1e-3,            #Toleranz zur Bestimmung der Integrationsschrittweite
                     optimize = true,          #Optimierung der DGL
                     maxfun = 0,               #Anzahl der Integrationsschritte -- 0 = keine Obergrenze
                     output = operator,        #Lösungsfunktionausgabe: damit Definition als Funktion möglich 
                     stiff = true:

Lösung_Bewgl := dsolve({Bewegungsgl_mit_Para, Anfangsbed}, numeric, Einst);

The equation to solve is this one: http://imgur.com/a/SFxI8

What am I doing wrong?

Edit:

This is the worksheet, problem is the last point.

WS.mw

Hi,

Please help me to solve numerically and visualise the following time-delayed integro-differential equation (see also the file attached)

diff(x(t), t$2) = .8/(x(t)^3*exp((2*0.1e-1)*(int(1/x(t), t = 0 .. t-tau))))-1/x(t)^2 , 

with the initial conditions x(0) = 1, (D(x))(0) = 0, and tau=0.1 (variable parameter).

The expected solution is an oscillatory function.

Thanks for your time and advice!

IDE.mw

hi--how i can solve following equation?

thanks

Eq.mw


Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/Eq.mw .
 

Download Eq.mw

Hello,

i have a Problem with maple. I have written this simple procedure

myPi_1:=proc(r)
local i,l;
l:=0;
for i from 1 to floor(r) do
    if isprime(i) and (i mod 4 = 1) then
        l:=l+1;
    end if;
end do;
return l;
end proc:    

When i use it like this 

k := myPi_1(80);

it works

But when i triy to plot it

plot(myPi_1(t), t = 0 .. 500, numpoints = 500, thickness = 2, color = black);

i get this message
Warning, expecting only range variable t in expression myPi_1(t) to be plotted but found name myPi_1

And no plot

Thank you for any Help

I want to calculate the following integral numerically with required precision.

First, the functions are defined:

G1:=-0.9445379894;
f:= (x) -> 0.9/abs(x-0.4)^(1/3)+0.1/abs(x-0.6)^(1/2);
U1 := unapply(-exp(-x)*(evalf(Int(f(t)*exp(t), t = 0 .. x))+G1)/2-exp(x)*(evalf(Int(f(t)*exp(-t), t = 0 .. x))+G1)/2, x);
U:= unapply(-exp(x)/2*(evalf(Int(f(t)*exp(-t),t=0..x))+G1)+exp(-x)/2*(evalf(Int(f(t)*exp(t),t=0..x))+G1), x);

Next, I calculate the integral in numerical form:

evalf(Int(U1(x)^2+U(x)^2-2*f(x)*U(x), x=0..1, digits=4, method = _Gquad));

If I specify digits=4, Maple return the answer -0.4291

If I use digits=5 or larger, Maple return someting like this

Is it possible to increase precision of calculation?

Hey there!

this is the formula:

diff(ca(t), t) = -3.600000000*10^20*exp(-15098.13790/(340-20*ca(t)))*ca(t), ca(0) = 2

I can easily plot it by solving DE numerically, however I can't seem to find a way to determine what 't' is at ca(t)=0.2

thank you so much in advance!

Kind regards, 

Gerard 

Hello People in mapleprimes,

I have a question about collecting part of expression.
Please teach me an answer to this.

I want to modify a1 to b1.


e3.mw

Thanks in advance.

For a 2D parametric plot with only one parameter, in Maple we use plot([x,y,range]) but what about the case x and y are expressed by more than one parameter? How should I assign the range? For example I tried the following one

h := k[1]*k[3]+k[1]*k[4]+k[2]*k[4]:
x := k[2]/h:
y := (k[3]+k[4])/h:
A := seq(k[i] = 1 .. 2, i = 1 .. 4):
plot([x, y, A]);

But then I encountered the following error:
Error, (in plot) incorrect first argument [k[2]/(k[1]*k[3]+k[1]*k[4]+k[2]*k[4]), (k[3]+k[4])/(k[1]*k[3]+k[1]*k[4]+k[2]*k[4]), k[1] = 0.1e-2 .. 1000, k[2] = 0.1e-2 .. 1000, k[3] = 0.1e-2 .. 1000, k[4] = 0.1e-2 .. 1000]

For one parameter the above method doesn't make an error like the following:

x := t+1:
y := t^2:
A := t = 1 .. 2:
plot([x, y, A]);

So for more than one parameter, there should be another way to ask the plot. Does anyone know how to do that?