Sort of revising fourier. I have a graph  however it gives result from -2pi to + 2pi. there must be a comand to modify plot to go from 0 to 2pi can someone give me a clue.

Rick

Got a question regarding numeric indices in tables. Apparently every index in a table is considered alfanumeric, which means that 90 is not the same as 90.0. 

Or test["text", 90] is not the same as test["text", 90.0].

Here's the example.

Table.mw

Apparently I will need to do some numeric formatting before using the numbers as indices for the table, right?

I have a graph with 267 nodes and 727 edges. I am trying to produce an eligible visualization of my graph but I'm getting a very crowded network. 

Code:

with(GraphTheory);
DrawGraph(Graph(M1), style = spring, symbolsize = 10, scaling = unconstrained); 

I have attatched the a JPEG file to my question to show what it looks like. I would like to space out the nodes with a more zoomed out view. How would I go about this?
 

I am trying to graph x^(1/2) and (x^2/8) revolving around the y-axis using the ring method but I have not been able to do it. I only get the shell method. I have adjusted the functions to read x=y^2 and x=sqrt(8y). The bounds for this method is 0 to 2. No matter the commands I have tried I have not been able to get the washer method to show. I do not want to use the tutor program.

Dear all

I applied the definition of a permutation matrix to a given matrix A to obtain a Tridiagonal Toeplitz matrix, but I haven't any idea why the final result looks strange.

Permutation_matrix.mw

Many thanks

d/dt(n_2)+d/dx (n0u2)

Thanks

Noticed an issue with Vector and excel.  Not even really excel but I found it through that.  

If I select a row in excel and do a copy paste into Maple it sometimes appears as if it's a list.  ie it shows up as [1,2,3,4]  other times it shows up as [ 1 2 3 4 ] signifying an array (for this example I've saved as the data under the variable a).  Running whattype(a) on the data shows they're both arrays.  Of course it's probably wise to instead load data through Excel into Maple by using the ExcelTools and Import option, however I didn't think there was no reason to use the copy paste option.  It does, to some extent, or maybe, well, it shows a limitation.

In the first instance when it shows up as a list, I thought it might be simple enough just to add other numbers in the list ie/ just tack on ,5,6,7,8 after the 4.  But after pressing enter we get the error

Error, invalid sequence

I don't understand why I couldn't just add to the copied data but maybe it's the formatting of the data that was pasted into maple? 

With the original copied data, I can use with(plots):  listplot(a) no problem however, we want to add data.  ArrayTools will let us Append data to the vector

with(plots):
with(ArrayTools):
Append(a,5,6,7,8)

listplot(a) # will throw another error

                  Error, (in plots:-pointplot) points cannot be converted to floating-point values (this doesn't show in Maple 2019 but does come up in Maple 18)

For Maple 2019 it with responds with a returned array [ 1 2 3 4 5,6,7,8]

Running listplot on that Appended Array.  Maple only recognizes the next 1 added value.

So this means appending values using ArrayTools and Append only allows us to append one value at a time?  Is that a bug?  Are there better ways to append values to vectors?  The only solution for that is to convert to a list and use op to append more values.  But perhaps Maple should have understood the copied values from Excel should have been a list or listlist of values?

Can someone help me with the following question?

I would like to plot the function f(x)=a*x*(1-x)+x*ln(x)+(1-x)ln(1-x) in the interval [x1,x2], where x1 and x2 are points that satisfy the following: df/dx(x=x1)=df/dx(x=x2)=0 and f(x1)=f(x2), where a is some parameter that is greater than 0.

I guess I need to use here fsolve, but I am not sure how.

Can you please tell me the exact weighting scheme for LinearFit() and NonlinearFit() when using the weights= option?

Is it wi*(yi_exp-yi_calc)^2 or wi^2*(yi_exp-yi_calc)^2 ?

(wi=weight of point i, yi_exp=value i measured, yi_calc= value i calculated from the model function)

Best regards,

 Anthrazius

Any integer in [129,129+13^2[ 1s written as the sum of squares of separate integers equal to or less than 12.
examples:: 129=10^2+5^2+2^2, 130=11^2+3^2... 132=9^2+5^2+4^2+3^2+1^2...
179=12^2+5^2+3^2+1^2...297=12^2+10^2+7^2+2^2. Thank you.

Hompotopy perturbation method

restart; with(LinearAlgebra);
PDEtools[declare](f(x), prime = x);
PDEtools[declare](g(x), prime = x);
N := 3;
F := sum(p^i*f[i](x), i = 0 .. N);
G := sum(p^i*g[i](x), i = 0 .. N);
FEq := (1-p)*(diff(F, x$3))+p*(diff(F, x$3)-(2*(diff(F, x$1))*(diff(F, x$1))-(5/2*F)*(diff(F, x$2))+M*(diff(F, x$1))) . A);
GEq := (1-p)*(diff(G, x$2))+p*(diff(G, x$2)+(5/2*(1/(1+R)) . Pr)*F*(diff(G, x$1)))*K(L);

M := 1;
A := (((1-W)^2.5*(1-W))*`ρf`+W*`ρt`)/`ρf`;
Pr := 6.2;
K := ((1-W)*`ρf`*Cf+W*`ρt`*Cs)/`ρf`;
L := .5;
`ρf` := 998.3;
`ρt` := 3970;
W := .2;
R := 1.0;
Cs := 765;
Cf := 4182;
coeff(FEq, p, 0);
coeff(GEq, p, 0);
for i from 0 to N do Fequ[i] := coeff(FEq, p, i) = 0 end do;
for i from 0 to N do Gequ[i] := coeff(GEq, p, i) = 0 end do;
Fcond[1][0] := f[0](0) = 0, (D(f[0]))(0) = 1, (D(f[0]))(5) = 0; for j to N do Fcond[1][j] := f[j](0) = 0, (D(f[j]))(0) = 0, (D(f[j]))(5) = 0 end do;
Gcond[0] := g[0](0) = 1, g[0](5) = 0; for j to N do Gcond[j] := g[j](0) = 0, g[j](5) = 0 end do;

for i from 0 to N do dsolve({Fequ[i], Fcond[1][i]}, f[i](x)); f[i](x) := rhs(%) end do;
for i from 0 to N do dsolve({Gcond[i], Gequ[i]}, g[i](x)); g[i](x) := rhs(%) end do;
Fa := simplify(sum(f[n](x), n = 0 .. N)); dFa := diff(Fa, x); subs(x = 2.4, dFa);
Ga := simplify(sum(g[n](x), n = 0 .. N)); dGa := diff(Ga, x); subs(x = 2.4, dGa);
plot(Ga, x = 0 .. 5);
plot(dFa, x = 0 .. 5);

ND sove solution 

Eq1 := diff(F(x), x$3)-(2*(diff(F(x), x$1))*(diff(F(x), x$1))-(5/2*F(x))*(diff(F(x), x$2))+M*(diff(F(x), x$1))) . A = 0;
Eq2 := (diff(G(x), x$2)+(5/2*(1/(1+R)))*Pr*F(x)*(diff(G(x), x$1)))*K(L) = 0;
M := 1;
Pr := 6.2;
A := (((1-W)^2.5*(1-W))*`ρf`+W*`ρt`)/`ρf`;
K := ((1-W)*`ρf`*Cf+W*`ρt`*Cs)/`ρf`;
L := 1;
`ρf` := 998.3;
`ρt` := 3970;
W := .2;
R := 2.0;
Cs := 765;
Cf := 4182;

Cd1 := F(0) = 0, (D(F))(0) = 1, (D(F))(5) = 0;
dsys := {Cd1, Eq1};
dsol := dsolve(dsys, numeric, output = operator);
dsol(.1);
plots[odeplot](dsol, [x, diff(F(x), x$1)], 0 .. 5, color = green);
Cd2 := G(0) = 1, G(5) = 0;
dsys := {Cd1, Cd2, Eq1, Eq2};
dsol := dsolve(dsys, numeric, output = operator);
plots[odeplot](dsol, [x, G(x)], 0 .. 5, color = green)

I find some book teach motion planning in topology 

but do it need to formulate the equations for the environments such as a map with obstacle in 2d or 3D?  the environment is quite complex, how can these equations be formulated?

Dear all

Using maple I can't obtain an exact solution, my goal for which initial condition can the solution diverge to infinity if maple can not give us a solution.

special_condition_diverge.mw

many thanks

I have two functions:

f(x)= (x)^(1/2)

g(x)=(x^2)/8

I want to get each function to rotate the y-axis and build the shape and I want the shaded region to build the 3D shape as it goes around the y-axis.

The domains for all would be 0 to 4.

The url link to what I am looking for it below.

https://en.wikipedia.org/wiki/Solid_of_revolution#/media/File:Rotationskoerper_animation.gif

Thank you in advanced!

Maple_question.mw

Hello everyone .

I asked Maple to compute a definite integral of a function, numerically. Unfortunately it does not give me any values . I briefly explain what I did.

As you can see in the picture first I define a real valued function BW_dilog(x) (with complex domain), where the function dilog  already exists in maple. Then I defined another function, V, (with real domain and codomain ). Now I need the approximate value of the bi integral of V, on a triangle. The Maple file is attached.

I would be grateful if any one could help me.