Dear all

I have the following proc I, I need a help so that the proc run without error. 

I get Error, invalid =

Code_eigenvalues_eigenvectors.mw

Thank you for any help 

Determine the sert of point  Z in complex plane isuch that  l2z-il=Im(z+1-i)

Lets say I wanna write a procedure DrawSpaceCurve3d. 

DrawSpaceCurve3d:=proc(fnc::algebraic, vars:name,h::integer, xvalue::range=a{integer}...b..{integer})
plots:-spacecurve([vars, h, funk(vars, h)], ':-x' = 'xvalue')
end proc; 

How I try to run this procedure, then I get the following error. 

"Error, (in Plot:-SpaceCurve) parameter range in the form name=range is missing"

So my question is it possible to get Maple to accept that I write DrawSpaceCurve3D(x^2*y, 2, xvalue=[-5..5]) 

I know the y - curve is missing, but I would like to see if I can get Maple to accept this kind of input?
 

I am trying to solve three simultaneous PDE where the first two PDEs are 1D while the third is 2D. When using pdsolve with numeric option I am getting the following error

Error, (in pdsolve/numeric/process_PDEs) PDEs can only contain dependent variables with direct dependence on the independent variables of the problem, got {Tg(t, z, 0.6985e-1)}. Can someone please help me with this.

restart;
T_well := 10. + 0.026*z;

PDE_in := -0.493381*diff(Ti(t, z), z) + diff(Ti(t, z), t) = 0.000176303*(-Ti(t, z) + To(t, z));
PDE_out := 0.186546*diff(To(t, z), z) + diff(To(t, z), t) = 0.0397694*(Tg(t, z, 0.06985) - To(t, z)) + 0.0000666597*(Ti(t, z) - To(t, z));
PDE_g := 0.22828*10^7/3.5*diff(Tg(t, z, r), t) = diff(Tg(t, z, r), r)/r + diff(Tg(t, z, r), r, r) + diff(Tg(t, z, r), z, z);
PDE := {PDE_g, PDE_in, PDE_out};
IC := {Tg(0, z, r) = T_well, Ti(0, z) = T_well, To(0, z) = T_well}
BC := {3.5*D[2](Tg)(t, 2000, r) = -3/40, 3.5*D[2](Tg)(t, z, 0.06985) = 0.0397694*(Tg(t, z, 0.06985) - To(t, z)), Tg(t, 0, r) = 10, Tg(t, z, 50) = 10 + 0.026*z, To(t, 0) = 10, To(t, 2000) = Ti(t, 2000)}
pdsolve(PDE, IC, BC, numeric)

with(plots);
a := -Pi;
b := Pi;
implicitplot3d(1/2*(sin(2*x)*cos(y)*sin(z) + sin(2*y)*cos(z)*sin(x) + sin(2*z)*cos(x)*sin(y)) - 1/2*(cos(2*x)*cos(2*y) + cos(2*y)*cos(2*z) + cos(2*z)*cos(2*x)) + 0.15 = 0, x = a .. b, y = a .. b, z = 1/2*a .. 3/2*b, grid = [70, 70, 70], style = patchnogrid, shading = xy, lightmodel = light2);
 

Hi! Is it possible to get more clear mesh plot when the grid is such large? 

Hello,

I am working on a project to code a conjecture in Maple and I am stuck.  It is as follows:

"Someone had a conjecture that any odd number greater than one can be written as the sum of a prime number and twice a perfect square.  (ie 3=3+2*0², 15 = 7+2*2²)  I would like to know if this holds for all the numbers from 3 to 9,000 and if not, which odd numbers fail this criteria."

I have very little coding experience so I am looking to see if anyone has an idea where I can go with this.  I tried the following:

X=P+2*N where

X:=seq(3+2i, i=0..4499);

P:={seq(ithprime(i), i=1..1250)};

N:=seq(i^2, i=0..100);

Once I get this far, however, I don't know where to go.  I would like the left side of the equation to run from 0 through 4499 (giving me a value of over 9000) while the right side does the same thing.  If there are any situations where the left side and right side do not equal each other, I want it to tell me.

My apologies if this is too wordy or confusing.  Thank you.

Dear all,

I cannot quite seem to figure out what it is that I am doing wrong. In the plots that I am trying to plot, I keep getting an error "Error, bad index into Matrix". I understand what the error means (having read through the webpage that it takes me to), but I cannot figure out what my mistake is. Can someone please help me understand what I am doing wrong?

Simple_Plot.mw

In order to use a physical constant I have to load the package "ScientificConstants". For instance, if I need the speed of light, c, I have to assign c:= Constant(c). In other words just by loading ScientificConstants the constants are not automatically assigned. Is that correct?

Hi folks

I'm trying to integrate products of Bessel functions and reproduce known results in the literature. For example, the integral

should evaluate to 1/q if q>k and 0 if q<k, but maple only gives the first solution, even when 'assuming q<k' is specified. It should give 1/2q when k=q -- which it does but only when I put k=q in the integral. For some reason, it always assumes q>k whatever I do. 

There's a lot of examples like this. Am I missing something? Is there a way to return an answer with all the (correct!) solutions?

Thanks in advance.

I am trying to solve (unsuccessfully) the title equation with MAPLE 2020.1

If someone could help me...

Here is my code, inspired by the above code :)

restart;
with(plots);
with(PDEtools);
a := 1;
b := 1;
rho := 8166;
h := 10^(-3);
Dp := 19.2615;
f := sin(7*x)*cos(4*y);
WavEq := [rho*h*diff(w(x, y, t), t $ 2)/Dp + diff(w(x, y, t), x $ 4) + 2*diff(w(x, y, t), x $ 2, y $ 2) + diff(w(x, y, t), y $ 4) = 0, u(0, y, t) = 0, u(a, y, t) = 0, u(x, 0, t) = 0, u(x, b, t) = 0, u(x, y, 0) = f];
sol := pdsolve(WavEq);

This get the following error message : "Error, (in pdsolve/sys) the input system cannot contain equations in the arbitrary parameters alone; found equation depending only on _F1(y,t): _F1(y,t)"

Hello, 

is there any possibility how to write just a part of axis label e.g. in Italic? For example if I want want a label to be "ln(h)" and i want just the "h" to be written in Italic font.

Thanks a lot.

Today I got some useful  graph datas  in the form of adjacency list like following:

I originally planned to enter them in maple one by one manually,

but I will consider a lot of graph datas, I don’t know if there is a program that can help, Thanks advance!

 

Dear all

I run a code But I obtained the following error 

Error, (in PW) cannot determine if this expression is true or false: (1/4)*10^(1/2) < 0.1000000000e-3


and  how can I display all iterations obtained from  PW and WDEF2 in a table.

PWW.mw

Thank you for any help 

Ennep := <u - u^3/3 + u*v^2, v - v^3/3 + v*u^2, u^2 - v^2>

plot3d(Ennep, u = -2 .. 2, v = -2 .. 2, shading = zhue, lightmodel = light1, orientation = [89, 54])

Do you know how I can make contours for this plot? Thanks in advance