I'm receiving a very interesting error when attempting to solve a system of differntial equations:

{sqrt(2)*sqrt(F1(t))*diff(F2(t),t) = 0}:
% union {F3(t) = 0}:
pdsolve(%, {F1(t),F2(t),F3(t)});

When attempting to run this code, I get the following error:

Error, (in pdsolve/sys) invalid input: DifferentialAlgebra:-DifferentialRing expects value for keyword parameter blocks to be of type {undefined, Or(list({name, unknown, list({name,unknown})}),name,unknown)}, but received [F3, F2, _F1 = (NULL), F1 = [_F1, F1], _F2] 

Interestingly, removing the second line, or even the sqrt(2), allows the calculation to proceed correctly. It this a bug, or am I doing something wrong?

Thank you very much!

(I'm aware this would technically be a system of ODEs, not PDEs; my actual code includes additional partial derivatives. Additionally, dsolve appears to only generate a trivial solution in this case).

Hey folks.

I need to do a kind of loop calc. in maple and im not shure how to. 

lets say i guess a temperature T.

i then put it through all soughts of eqations where the variabel T is necessary for the outcome of my result, which i then use again as T. the loop should finish when T as a result is the same as my guess in the beginning. like:

T:=45 

bla bla calculation

T_new:= 35

again 

T:=35

bla bla calculation

T_new= 35 

Finish!!

Any ideas?

Dear sir I want to plot the graph for different values of F=0.1, F=0.2,F=0.5 . My codes are given and i am attaching the sample graph. For F=0.1, Nb=0.1,0.2,0.3 three curves

Similarly For F=0.2,Nb=0.1,0.2,0.3 three curves and for F=0.5, Nb=0.1,0.2,0.3 three curves. Total nine curves 

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):

L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);

please how can I simplify the results of assigned K cum each of the coefficients further in the following code? the results are too large.

restart:
#k=2:
omega:=u/h:
psi:=v/h:
t:=(sum(a[j]*x^j,j=0..2)+a[3]*sin(omega*x)+a[4]*cos(omega*x)+a[5]*sin(psi*x)+a[6]*cos(psi*x)):
F:=diff(t,x):
G:=diff(t,x,x):
p1:=simplify(eval(t,x=q+h))=y[n+1]:
p2:=simplify(eval(F,x=q))=f[n]:
p3:=simplify(eval(F,x=q+h))=f[n+1]:
p4:=simplify(eval(F,x=q+2*h))=f[n+2]:
p5:=simplify(eval(G,x=q))=g[n]:
p6:=simplify(eval(G,x=q+h))=g[n+1]:
p7:=simplify(eval(G,x=q+2*h))=g[n+2]:
vars:= seq(a[i],i=0..6):
Cc:=eval(<vars>, solve({p||(1..7)}, {vars})):
for i from 1 to 7 do
	a[i-1]:=Cc[i]:
end do:
Cf:=t:

K:= collect(combine(simplify(eval(Cf,x=q+2*h),size),trig),{y[n+1],f[n],f[n+1],f[n+2],g[n],g[n+1],g[n+2]},factor):


alpha[1]=simplify(coeff(K,y[n+1]));
beta[0]=simplify(coeff(K,f[n]),size);
beta[1]=simplify(coeff(K,f[n+1]),size):
beta[2]=simplify(coeff(K,f[n+2]),size):
gamma[0]=simplify(coeff(K,g[n]),size):
gamma[1]=simplify(coeff(K,g[n+1]),size):
gamma[2]=simplify(coeff(K,g[n+2]),size):

I need to plot a trajectory where the horizontal equation is x(t) = 4t and the vertical is y(t)=4t-0.5t^2

How would I plot the two equations at once to make a trajectory?

If someone could please explain how to do it (also doing it would be nice) because I'm sure I'll have to do it again.

How do I apply select to a set element-wise?

I have defined a procedure which computes the order of an element g mod p.

findOrderOf:=proc(g,p)
  local i:=1;
  local gpwr:=1;
  for i from 1 to p-1 do
    gpwr:=gpwr*g mod p;
    if(gpwr = 1) then:
      return i;
    end if:
  end do:
end proc:

I then want to select the primes which are generated by 2 and 3.

primes:=select(isprime,[`$`(2000 .. 3000)]);
select(findOrderOf(2,x=p)=p-1 and findOrderOf(3,x=p)=p-1,primes);

I am stuck with how to make the last line work.

hi

please help me for simplify (factor) this equations.

thanks

vel.mw
 

Download vel.mw

can maple code saved into database for self modifying easily?

self modifying , metaprogramming , database, snapshot, blockchain for code integrity etc

LetterA_Axes_query.mws

  The attached shows animation of the letter A. The animation works fine, and in the plot diagram of the large letter A, the axes are shown.  However, in the second plot (animated), the axes are not being shown, despite clicking on Axes in the animation.  The choices here are Boxed, Framed, Normal or None.

  In the Maple 7 documentation on Axes it states "From the Axes menu choose Ranges...."   ...but I can't see Ranges. 

Any reasons for this?

Hi all, since I cannot use Latex over here I try to formulate as good as possible. Suppose we have the standard integral from a to b of a function f(x), which is equal to:

h * sum (k=1 to N) [c_k f(a+kh)] +O(h^p), where h=(b-a)/N+1, an p(N)>= N+1;

Would it be possible to write a procedure to evaluate the coefficients c_k , k=1,2,3..N and P(N)

Perhaps a stupid question, because I don't know if it is possible.

Best wishes, Math

Hey, 

I run into the same issue quite often. Let's say I have some polynomial over variable like this

I do a substitution on s and get something quite messy as a result. Example:

Now I would like to get the polynomial in the same pole / zero form. I'm sure maple is capable of doing this but I have no clue where to click. Is there an easy way to bring a polynomial over some variable in a standarized form?

Cheers!

Download simplln.mw

with(CodeGeneration)

ua=[ 1 2 3 4];

Matlab(ua);

%How can i write the Output of "Matlab(ua)"  on file so I can use the Matlab "load" command
later on to get the result in Matlab?  

For my own use, I am attempting to port Joe Riel’s glyph package for geometric algebra into a module more compatible with recent versions of Maple. To this end, I have been testing individual procedures extracted from the package into Maple 2016, both to understand the algorithms and to check for glitches caused by the code running in more current Maple 2016. The procedure for carrying out the exterior multiplication of blades does not seem to work reliably, and I haven’t the necessary knowledge of Maple language to determine whether this is due to an error on my part or a feature of Maple V that no longer works.  I have attached a worksheet,tablemultiplyexample.mw,  that includes the procedures necessary for the multiplication routine to work, but I can’t get any consistency in the results.  Can anyone advise me what is the problem?  

As I understand the routine, setup defines a anti-symmetric root blade table with an indexing function that precludes assignment to the table. Clifford blades are then represented as indexed variables using the root table. The process is as follows see worksheet for actual code):

initialize := proc ()
 global _e, tableroot;
tableroot := table(antisymmetric, blade);
tableroot[] := 1;
_e := tableroot;
end proc:
#The index function `blade` is as follows:

`index/blade` := proc (Indices, tableau)
 if nargs = 2 then if Indices = [] then 1
          else tableau[op(checkindices(Indices))] end if
elif Indices = [] then tableau[Indices[]] := 1
 else ERROR("cannot assign to a blade", Indices) end if
end proc;

#Exterior multiplication is performed by the following routine.
b_exteriorp := proc (u, v)
option remember;
 if  u = 1 or v = 1 then u*v
else _e[op(u), op(v)] end if
end proc:

As near as I understand, the procedure joins the lists representing the two input blade into a single list that is processed by the antisymmetric indexing function and outputs the indexes as the product blade. I don’t understand how the case of duplicate indexes (which should return 0) is supposed to be handled by the procedure.  What the procedure usually returns is simply the appended list of the two blades without modification by the indexing function.

Can anyone give me a hint about how to fix this procedure?

tablemultiplyexample.mw

I am trying to input data via a data table. Have several problems here.

I used an array because I want the row and column numbers to start at 0.

1st When the table appers after that the document runs hediously slow as in a second or two to enter a digit or letters appear after typing.  Like something is absorbing the computer resources. But I have a fast machine.

2nd Any data I  enter to the table vanishes but does get stored.

3rd I tried to turn it all into a procedure but cant get that to work.
 

Download DataTable_Experiment.mw