So I have this system of equations with which I am not sure if the result is the same or not using "series" and "limit" or what is going on here.

I hope it is clear what I mean.

Download Mapleprimes_-_Ionproduct.mw

restart;
N:=4;alpha:=5*3.14/180;r:=10;Ha:=5;H:=1;
dsolve(diff(f(x),x,x,x));
Rf:=diff(f[m-1](x),x,x,x)+2*alpha*r*sum*(f[m-1-n](x)*diff(f[n](x),x),n=0..m-1)
+(4-Ha)*(alpha)^2*diff(f[m-1](x),x);
dsolve(diff(f[m](x),x,x,x)-CHI[m]*(diff(f[m-1](x),x,x,x))=h*H*Rf,f[m](x));
f[0](x):=1-x^2;
for m from 1 by 1 to N do
CHI[m]:='if'(m>1,1,0);
f[m](x):=int(int(int(CHI[m]*(diff(f[m-1](x),x,x,x))+h*H(diff(f[m-1](x),x,x,x))
+2*h*H*alpha*r*(sum(f[m-1-n](x)*(diff(f[n](x),x)),n=0..m-1))+4*h*H*alpha^2*
(diff(f[m-1](x),x))-h*H*alpha^2*(diff(f[m-1](x),x))*Ha,x),x)+_C1*x,x)+_C2*x+_C3;
s1:=evalf(subs(x=0,f[m](x)))=0;
s2:=evalf(subs(x=0,diff(f[m](x),x)))=0;
s1:=evalf(subs(x=1,f[m](x)))=0;
s:={s1,s2,s3}:
f[m](x):=simplify(subs(solve(s,{_C1,_C2,_C3}),f[m](x)));
end do;
f(x):=sum(f[1](x),1=0..N);
hh:=evalf(subs(x=1,diff(f(x),x))):
plot(hh,h=-1.5..-0.2);
A(x):=subs(h=-0.9,f(x));
plot(A(x),x=0..1);

A parameterization of the function i am studying this morning produced what seems to be not single valued on the domain i choose,

plot(floor((x+1)^2/x^2)-(x+1)^2, x = -10 .. 10, coords = logarithmic);

plot(floor((x+1)^2/x^2)-(x+1)^2, x = -10 .. 10, coords = logcosh);

but the cartesian is fine:

Dear Users!

I am facing problem to compare the coefficient of x^i*y^j for i, j =1..,Equations. Please my effort and fix the problem.

H1 := 3*y^4*a[1]^5*b[1]+6*y^4*a[1]^3*b[1]^3+3*y^4*a[1]*b[1]^5+6*x*y^3*a[1]^5+6*x*y^3*a[1]^4*b[1]+12*x*y^3*a[1]^3*b[1]^2+12*x*y^3*a[1]^2*b[1]^3+6*x*y^3*a[1]*b[1]^4+6*x*y^3*b[1]^5+6*y^3*a[1]^5*b[2]+6*y^3*a[1]^4*a[2]*b[1]+12*y^3*a[1]^3*b[1]^2*b[2]+12*y^3*a[1]^2*a[2]*b[1]^3+6*y^3*a[1]*b[1]^4*b[2]+6*y^3*a[2]*b[1]^5+18*x^2*y^2*a[1]^4+36*x^2*y^2*a[1]^2*b[1]^2+18*x^2*y^2*b[1]^4+18*x*y^2*a[1]^4*a[2]+18*x*y^2*a[1]^4*b[2]+36*x*y^2*a[1]^2*a[2]*b[1]^2+36*x*y^2*a[1]^2*b[1]^2*b[2]+18*x*y^2*a[2]*b[1]^4+18*x*y^2*b[1]^4*b[2]+18*y^2*a[1]^4*a[2]*b[2]+36*y^2*a[1]^2*a[2]*b[1]^2*b[2]+18*y^2*a[2]*b[1]^4*b[2]-5*delta*y^2*a[1]^4-8*delta*y^2*a[1]^3*b[1]-10*delta*y^2*a[1]^2*b[1]^2-8*delta*y^2*a[1]*b[1]^3-5*delta*y^2*b[1]^4+12*x^3*y*a[1]^3+12*x^3*y*a[1]^2*b[1]+12*x^3*y*a[1]*b[1]^2+12*x^3*y*b[1]^3+36*x^2*y*a[1]^3*a[2]+36*x^2*y*a[1]^2*b[1]*b[2]+36*x^2*y*a[1]*a[2]*b[1]^2+36*x^2*y*b[1]^3*b[2]+18*x*y*a[1]^3*a[2]^2+36*x*y*a[1]^3*a[2]*b[2]-18*x*y*a[1]^3*b[2]^2-18*x*y*a[1]^2*a[2]^2*b[1]+36*x*y*a[1]^2*a[2]*b[1]*b[2]+18*x*y*a[1]^2*b[1]*b[2]^2+18*x*y*a[1]*a[2]^2*b[1]^2+36*x*y*a[1]*a[2]*b[1]^2*b[2]-18*x*y*a[1]*b[1]^2*b[2]^2-18*x*y*a[2]^2*b[1]^3+36*x*y*a[2]*b[1]^3*b[2]+18*x*y*b[1]^3*b[2]^2+18*y*a[1]^3*a[2]^2*b[2]-6*y*a[1]^3*b[2]^3-6*y*a[1]^2*a[2]^3*b[1]+18*y*a[1]^2*a[2]*b[1]*b[2]^2+18*y*a[1]*a[2]^2*b[1]^2*b[2]-6*y*a[1]*b[1]^2*b[2]^3-6*y*a[2]^3*b[1]^3+18*y*a[2]*b[1]^3*b[2]^2-16*delta*x*y*a[1]^3-20*delta*x*y*a[1]^2*b[1]-20*delta*x*y*a[1]*b[1]^2-16*delta*x*y*b[1]^3-10*delta*y*a[1]^3*a[2]-6*delta*y*a[1]^3*b[2]-10*delta*y*a[1]^2*a[2]*b[1]-10*delta*y*a[1]^2*b[1]*b[2]-10*delta*y*a[1]*a[2]*b[1]^2-10*delta*y*a[1]*b[1]^2*b[2]-6*delta*y*a[2]*b[1]^3-10*delta*y*b[1]^3*b[2]+12*x^4*a[1]*b[1]+12*x^3*a[1]^2*a[2]-12*x^3*a[1]^2*b[2]+24*x^3*a[1]*a[2]*b[1]+24*x^3*a[1]*b[1]*b[2]-12*x^3*a[2]*b[1]^2+12*x^3*b[1]^2*b[2]+18*x^2*a[1]^2*a[2]^2-18*x^2*a[1]^2*b[2]^2+72*x^2*a[1]*a[2]*b[1]*b[2]-18*x^2*a[2]^2*b[1]^2+18*x^2*b[1]^2*b[2]^2+6*x*a[1]^2*a[2]^3+18*x*a[1]^2*a[2]^2*b[2]-18*x*a[1]^2*a[2]*b[2]^2-6*x*a[1]^2*b[2]^3-12*x*a[1]*a[2]^3*b[1]+36*x*a[1]*a[2]^2*b[1]*b[2]+36*x*a[1]*a[2]*b[1]*b[2]^2-12*x*a[1]*b[1]*b[2]^3-6*x*a[2]^3*b[1]^2-18*x*a[2]^2*b[1]^2*b[2]+18*x*a[2]*b[1]^2*b[2]^2+6*x*b[1]^2*b[2]^3+6*a[1]^2*a[2]^3*b[2]-6*a[1]^2*a[2]*b[2]^3-3*a[1]*a[2]^4*b[1]+18*a[1]*a[2]^2*b[1]*b[2]^2-3*a[1]*b[1]*b[2]^4-6*a[2]^3*b[1]^2*b[2]+6*a[2]*b[1]^2*b[2]^3-10*delta*x^2*a[1]^2-16*delta*x^2*a[1]*b[1]-10*delta*x^2*b[1]^2-16*delta*x*a[1]^2*a[2]-4*delta*x*a[1]^2*b[2]-16*delta*x*a[1]*a[2]*b[1]-16*delta*x*a[1]*b[1]*b[2]-4*delta*x*a[2]*b[1]^2-16*delta*x*b[1]^2*b[2]-5*delta*a[1]^2*a[2]^2-6*delta*a[1]^2*a[2]*b[2]+delta*a[1]^2*b[2]^2-2*delta*a[1]*a[2]^2*b[1]-12*delta*a[1]*a[2]*b[1]*b[2]-2*delta*a[1]*b[1]*b[2]^2+delta*a[2]^2*b[1]^2-6*delta*a[2]*b[1]^2*b[2]-5*delta*b[1]^2*b[2]^2+delta^2*a[1]^2+delta^2*a[1]*b[1]+delta^2*b[1]^2+16*y^2*a[1]^2+48*y^2*a[1]*b[1]+16*y^2*b[1]^2+80*x*y*a[1]+80*x*y*b[1]+32*y*a[1]*a[2]+48*y*a[1]*b[2]+48*y*a[2]*b[1]+32*y*b[1]*b[2]+80*x^2+80*x*a[2]+80*x*b[2]+16*a[2]^2+48*a[2]*b[2]+16*b[2]^2-8*delta;

Equation := 12;

for i from 0 to Equation do;

for j from 0 to Equation do;

C[i, j] := coeff(H1, x^i*y^j) = 0;

end do;

end do;

I got this error
Error, invalid input: coeff received 1, which is not valid for its 2nd argument, x
 

Dear Users!

Hope you would be fine. I want to write an expression in sigma notation which control ny n (any constant >0);
for n =1 expression expand as

E[1]+1

for n =2 expression expand as
E[1]*E[2]*a[12]+E[1]+E[2]+1;

for n =3 expression expand as

E[1]*E[2]*E[3]*a[123]+E[1]*E[2]*a[12]+E[1]*E[3]*a[13]+E[2]*E[3]*a[23]+E[1]+E[2]+E[3]+1;

for n =4 expression expand as

E[1]*E[2]*E[3]*E[4]*c[1234]+E[1]*E[2]*E[3]*a[123]+E[1]*E[2]*E[4]*a[124]+E[1]*E[3]*E[4]*a[134]+E[2]*E[3]*E[4]*a[234]+E[1]*E[2]*a[12]+E[1]*E[3]*a[13]+E[1]*E[4]*a[14]+E[2]*E[3]*a[23]+E[2]*E[4]*a[24]+E[3]*E[4]*a[34]+E[1]+E[2]+E[3]+E[4]+1;

and so on.

I am waiting your kind respons. Thanks

Hello,

how to calculate the laplace transform for the following equations?

L1:=laplace(psi1(t)*(diff(z1(t), t)), t, s):
L2:=laplace((diff(psi1(t), t))^2, t, s):

I can't final an equivalent to Mathematica's Flatten for sets. I know Maple has ListTools:-Flatten for lists.   

For example, given set r:={a,{b,c},d,{e,f,{g,h}}}; How to convert it to  {a,b,c,d,e,f,g,h}; 

does one have to convert each set and all the inner sets to lists, then apply ListTools:-Flatten to the result? How to map convert(z,list) for all levels?

     map(z->convert(z,list),r);

does not work, since it only maps at top level, giving {[a], [d], [b, c], [e, f, {g, h}]}

So doing

   ListTools:-Flatten(convert(map(z->convert(z,list),r),list));

Gives [a, d, b, c, e, f, {g, h}] 

Hey there, I'm trying to count how many letter arrangements are possible by using the algorithm below (for Maple TA). It's a bit crude but the console tells me countcharacteroccurences second argument must be a string, but it works for the earlier letters. Can someone please give me a bit of guidance?

$temp = maple("
randomize();
with(MathML):
with(StringTools):
with(combinat):

rintS := rand(1..27);
word := ARITHMETIC,ALGORITHM,ASYMPTOTE,AVERAGE,CARTESIAN,CALCULUS,COEFFICIENT,COORDINATE,NUMERATOR,
DENOMINATOR,DIFFERENTIATE,DERIVATIVE,DIAMETER,DYNAMICS,EXTRAPOLATION,FACTORIALS,GEOMETRIC,
HYPOTENUSE,INTEGRATION,IRRATIONAL,INVERSE,ITERATION,POLYNOMIAL,COMBINATIONS,PERMUTATIONS,POLYGON;

disp := word[rintS()];
n := length(disp);

n1 := CountCharacterOccurrences(disp,A);
n2 := CountCharacterOccurrences(disp,B);
n3 := CountCharacterOccurrences(disp,C);
n4 := CountCharacterOccurrences(disp,D);
n5 := CountCharacterOccurrences(disp,E);
n6 := CountCharacterOccurrences(disp,F);
n7 := CountCharacterOccurrences(disp,G);
n8 := CountCharacterOccurrences(disp,H);
n9 := CountCharacterOccurrences(disp,I);
a1 := CountCharacterOccurrences(disp,J);
a2 := CountCharacterOccurrences(disp,K);
a3 := CountCharacterOccurrences(disp,L);
a4 := CountCharacterOccurrences(disp,M);
a5 := CountCharacterOccurrences(disp,N);
a6 := CountCharacterOccurrences(disp,O);
a7 := CountCharacterOccurrences(disp,P);
a8 := CountCharacterOccurrences(disp,Q);
a9 := CountCharacterOccurrences(disp,R);
b1 := CountCharacterOccurrences(disp,S);
b2 := CountCharacterOccurrences(disp,T);
b3 := CountCharacterOccurrences(disp,U);
b4 := CountCharacterOccurrences(disp,V);
b5 := CountCharacterOccurrences(disp,W);
b6 := CountCharacterOccurrences(disp,X);
b7 := CountCharacterOccurrences(disp,Y);
b8 := CountCharacterOccurrences(disp,Z);

z1 := n1!*n2!*n3!*n4!*n5!*n6!*n7!*n8!*n9!;
z2 := a1!*a2!*a3!*a4!*a5!*a6!*a7!*a8!*a9!;
z3 := b1!*b2!*b3!*b4!*b5!*b6!*b7!*b8!;
ans := n!/(z1*z2*z3);

Export(disp),convert(ans,string),convert(z1,string),convert(z2,string),convert(z3,string);
");

$ans = switch(1,$temp);
$disp = switch(0,$temp);

In attempting to store procedures in a table to create an extensible module, I used the following procedure to test whether a procedure could be queried from the table and executed.

proc() print("Metric is default") end proc;

The result was proc () print("Metric is default") end proc

I have checked the statement with mint in the code editor, but it reports no errors. I get the same result when I use an eval statement around the procedure.  I would be very interested in understanding what is wrong with the procedure as defined.

Do implicit plots need to be strictly real valued because of the multivalued problem with complex functions? because i get this really funny output for this one, and if I increase Digits, it produces a blank axis for some reason.

Hi everyone, now I try to plot a graph by varying more than one variable. Is it possible to vary for more than one variable at a time (vary the two or more variable at one once) Please anybody can help in this regard?
 

Download MHD_cchf.mw
 

Download MHD_cchf.mw

This is a picture of what maple looks like on  my laptop with a 14 inch screen (which ahas a 3200*1800 resolution).

As you can see the icons at the top are small, I'd love to make them permanently bigger; however the larger (smaller) probelm  is the default text size on the plot. Its microscopic. Is there anyway to change it so by default the font is much bigger every time i make a plot?

 

Does anyone know why the following gives an error of invalid loop statement termination?

if true then
  for i from 1 to 2 do
  end do
else
end if

while putting a semicolon after `end do' fixes it.

This seems to only occur in if else statements, when end do appears directly before the else

I can't figure out why this pattern is failing

restart;
r:= y=3;
patmatch(r,y=z::integer,'la');
patmatch(r,identical(y)=z::integer,'la');
patmatch(r,y::anything=z::anything,'la');

What Am I doing wrong?

I try to calculate the alpha1,alpha2…which were showed in the reference paper.  Actually, the parameter of A and B is identical to the reference result, but the parameter alpha1,alpha2…… do not as same as shown in the reference paper. The attachment is my Maple code about this problem. Please tell me why this does not work.  By the way, I can not get the integration result in which the integral function contains the differential, such as  alpha3 and alpha4.
 

Download feijisuan.mwfeijisuan.mw