having its diameter equal to three and its number of vertices equal to 32? Is it possible to find the one with Maple? The simpleminded approach
>with(GraphTheory):
>with(RandomGraphs):
> G := RandomRegularGraph(32, 4);
> Diameter(G);
produces 4 or 5.
What I know is http://www-ma4.upc.es/~comellas/delta-d/desc_g/desc_g3.html ( the second graph).

I mean not open these files, and ctrl+A select whole file and copy to one file.

I have been trying to find a solution to the following equation for ages! 

Vt:=(a*t)+((1/(2*Pi))*((b*cos(2*Pi*t))-(c*sin(2*Pi*t))-b)-295);

I need to find what t converges to, I have solved the equation using a different method and I know the answer is around 0.787 but I need to find the result using the Newton Raphson method, or another iterative method. I am not allowed to use dsolve. or the Newton command.

Dear All,

I am trying to find the exact amplitude and phase angle for the signal, with values attached using MATLAB.  BUT the phase angles doesn't make sense.  Can you guys help me in doing the FFT of the signal and find the amplitude and phase using MAPLE.   I am interested in,  fundamental component, 5th , 7th , 11th, 13th 15th, 17th and 19th harmonics with their phases.

I shall be very thankful if anyone of you can actually solve the problem.

below code is calculate basis of kernel and kernel

i guess basis of image is 

remove(has, Ga, [r,u,v,w]); if this correct, i eliminate this, i can get the image

however it include variable 'a'

is it correct? if not, how to calculate? 

my final goal is to make unexact sequence into exact sequence

after union ideal I and J, how to remove the intersection of I and J in order to do primedecomposition?

I have a Matrix 'A' in maple with symbolic entries and I need to use that Matrix in matlab M file. The symbols with in 'A' are already initialized in M file.

How is this possible?

For Matrices 'A' of small dimension, I can use the codegeneration package to convert 'A' to Matlab Matrix and then copy paste.

But for matrices of 1000x1000, lets say, does not seem as a option.

r-x^4
u-(x^3)*y
v-x*y^3
w-y^4

Step 1
with(Groebner):
K := {r-x^4,u-(x^3)*y,v-x*y^3,w-y^4};
G := Basis(K, 'tord', degrevlex(r,u,v,w));

Step 2
G intersect k[r,u,v,w]

G2 := Basis({w*G[1],w*G[2],w*G[3],w*G[4],w*G[5],w*G[6],w*G[7],w*G[8],w*G[9],w*G[10],w*G[11],w*G[12],w*G[13],w*G[14], (1-w)*K[1], (1-w)*K[2], (1-w)*K[3], (1-w)*K[4]}, 'tord', degrevlex(r,u,v,w));

G2 has too much though include the correct solution

correct solution...

3*rho1 - 2*rho2 + rho3 - rho4 = -1

4*rho1 +   rho2 - rho3        = 5

original without cost function:

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,(y1^1000)*(y2^1)*(y3^1)*(y4^100)- x1*x2*w + 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));

after have cost function 1000*rho1 + rho2...

Here is my sheet. non_linear_P2_last_q.mw

In this problem I have achieved every thing what I aimed, thnaks to the wonderful people here!!!

I am looking for a way to do it in a more structured way meaning the calculation of Jacobian in the attached sheet is somewhat manuel, though I wonder does maple has some other way?

Hi all I have the following ODE

ODE_1:=diff(w(r),r)+R*(diff(w(r), r))^3=K/r;

ODE_T_1:=collect(algsubs(w(r)=u(r)+(r-1)*(U-0)/(delta-1), ODE_1),diff(u(r), r)) ;

 and when I try to do the following integration,

eq1:=int(phi[i](r)*ODE_T_1,r=1..delta) assuming delta > 1;

it gives the following error

"

Error, (in assuming) when calling 'int'. Received: 'wrong number (or type) of arguments:

i am doing a Q Laplace and finally summation it

however, i find an AccurateQSummation which is not like Summation

Should i use it, if not, use SumTools' summation how to simplify it?

restart;
with(QDifferenceEquations):
qexp := QPochhammer(-(-p*t), q, infinity);
f:= qexp*t;
b := subs(t=q^n,f);
QSimplify((1-q)*AccurateQSummation(b*q^n,n=-infinity..infinity));

with(SumTools):
QSimplify((1-q)*Summation(b*q^n,n=-infinity..infinity));

Hi

i know to make an asymptotic expansion one uses asympt(1+x,x); how do you make it go to -infity rather than plus infinity?

I have a huge PDE system included reaction-diffusion system with 4 independent and 12 independent variables. I read on the maple help for pdsolve/numeric/method. I found it is only for single partial differential equations in two independent variables. Is there anyone knows the ability of Maple for such a system?

1.

qRacah := qhyperterm([q^(-x), alpha*q^(-x)/(gamma*delta)], [alpha*q], q, gamma*delta*q^(x+1)*t, k)*qhyperterm([beta*delta*q^(x+1), gamma*q^(x+1)], [beta*q], q, q^(-x)*t, k);

i use inverse laplace to transform ordinary generating function into exponential generating function

can not further evaluated

2.

i do not know whether x is main variable or q is main variable or both are.

ff := eval(subs(x=0, exponential_genfun));