I have a proc f(a,b)

For example

m := n->piecewise(type(n,even),0, 1):
f := proc(a,b)
    m(a*b);
end;

when I execute f(a,b) I get 1!!

for some reason the type is returning that it is not even, even though the expansion needs to be delayed.

I did this because mod was also doing this. I replaced it so I could undefine m and actually get a general expansion in terms of m.

Is there any way to get it to delay the expansion? I guess I need to check if the type is a "variable" or not an itneger and just return it.

I've created

m := n->piecewise(type(n,even),0, type(n, integer), 1, M(n)):

which seems to work but this is not ideal because sometimes M is returned even if an integer is used when used in other expressions. I suppose it is how ever maple figures out how evaluate things.

Hi,

Sorry to ask such a stupid question but I can't find out where my error is. Probably it's so huge it blinds me!

The double loop and the matrix product F^+ . F should give the same result, no? (it seems that F^+ . F has its rows reordered ?)


 

Download Too_Blind_To_Find_My_Mistake.mw

I have several plots where I use display(a,b,c). Maple displays them with display(actual plot of a, actual plot of b, ...).
I've tried adding insequence but nothing and also tried passing them as sets or lists. In earlier versions I could use display(plots) and it would display all the plots on one plot together but now maple always separates them. I'd rather have the old behavior.

How do I get them to plot on the same plot and is there any way to set it as a global ooption?

what is the wrong in this operations for isentropic process ( isentropic-process) 

Restart: with(RealDomain) : interface(dispalyprecision=4) : ; Isentropical := proc(N,v) local K,PTR,KM1,KP1,GO, M,MS,C,PPT,TTT,RRT,AAS,eq,o ; K : = G() : KM1 : = k - 1 : KP1 : = K + 1 : GO := 1 : If N=1 then M := v Elif N=2 then PPT := v : M := sqrt ( 2 * (PPT^(-KM/K) -1 ) : Elif N=3 then MS := v : M := sqrt ( 1 / (KP1/ (2*MS^2) –KM1/2 ) : Elif N=4 then C := v : M := sqrt ( 2 / (KM1* (1 /(C^2) -1 )) ) : Elif N=5 then TTT := v : M := sqrt ( 2 * (1/TTT-1 ) / KM1 : Elif N=6 then RRT := v : M := sqrt ( 2 * (1/RRT^KM1 -1 ) / KM1 : Elif N=7 then AAS := v : eq :=AAS- (2*(1+KM1*M^2/2)/ KP1) ^ (KP1/ (2*KM1)) /M : M := fsolve(eq,M,0..1) : GO :=2 : Elif N=8 then AAS := v : eq :=AAS- (2*(1+KM1*M^2/2)/ KP1) ^ (KP1/ (2*KM1)) /M : M := fsolve(eq,M,1..infinity) : GO :=2 : Fi:fi GO<>2 then ASS := (2 *(1+ KM1*M^2/2/KP1) ^ (KP1/ (2*KM1)) /M fi: PTR := (1+ KM1*M^2/2/KP1) : PPT := PTR^(-K/KM1) : TTT:=PTR^ (-1) : RRT:+PTR^(-1/KM1) : If N <> 3 then MS := sqrt ( KP1 / ( 2 / ( M^2)+KM1) ) If N <> 3 then C := sqrt( 1 / (2 / (KM1*M^2) + 1 ) ) O := evalf([m,1./PPT,MS,C,1./TTT,1./RRT,AAS]) end : irp := proc (N, v) local f; f:=Isentropic1(N,v) : printf(cat(`% 9.4f `$7, ` \ n`) , f[ii$ii=1..7]) end : G := ( ) -> 1.4 : For mi form .1 to .5 do irp(1.mi) od : > > > Eq := table( [ (M2) = M2 , (MS2) = (K+1) /2 (2 / ( M2)+ (K-1) ) (C2) = 1 / (2 / (K-1)*M2+1 ) (PP0) = (1+(K-1)*M2/2)^(-K/(K-1) ) (RR0) = (1+(K-1)*M2/2)^(-1/(K-1) ) (TT0) = (1+(K-1)*M2/2)^(-1) (AAS) = (2 *(1 +(K-1)*M2/2 ) / (K+1) )^( (k-1) /(2*(k-1)))/ squt(M2) ]) : if vL=M then M2v := solve ( M2 -Eq[M2 ] , M2 ) : M2v := eval(M2v ,M2 = VR^2 ) : elif vL=MS then M2v := solve( M2 -Eq[MS2 ] , M2 ) : M2v := eval(M2v ,MS2 = VR^2 ) : elif vL=C then M2v := solve( M2 -Eq[C2 ] , M2 ) : M2v := eval(M2v ,C2 = VR^2 ) : elif vL=PP0 then M2v := solve( M2 -Eq[PP0] , M2 ) : M2v := eval(M2v ,PP0 = VR ) : elif vL=RR0 then M2v := solve( M2 -Eq[RR0] , M2 ) : M2v := eval(M2v ,RR0 = VR ) : elif vL=TT0 then M2v := solve( M2 -Eq[TT0] , M2 ) : M2v := eval(M2v ,TT0 = VR ) : elif vL=AASsub then M2v := fsolve( M2 -Eq[AAS] , M2 ) : M2 :=0..1 : elif vL=AASup then M2v := fsolve( M2 -Eq[AAS] , M2 ) : M2 :=1..infinity : fi: Eqsva1 := eva1(Eq, M2 = M2v) : Eqava1 [M] := sqrt(EqsVa1[M2] ) : Eqsva1[MS] := sqrt(EqsVa1[MS2] ) : Eqava1 [C] := sqrt(EqsVa1[C2] ) : 1var := [M, MS, C, PP0, RR0, TT0, AAS] : [seq(1var[i]=EqsVA1[1var[i] ] , I=1…nope(1var))] end: Egs :=[M=0.5, MS=0.53452, C=0.2182, P0P=0.84302, RR0=0.88517, TT0=0.95238, AASsub=1.33985, AASsup=0.95238] : For i from 1 to nops(Egs) do print(I,Isentropic2(Egs[i])) od : (process)

I have functions with sin(x) and  cos(x) terms with x values are degrees.

How to calculate sin(30 degree)?

How to insert degree symbol in maple?

Suppose evaluate  sin(Pi/2), my out put becomes sin(Pi/2)...How to solve this problem?

Hi,

Idea for trigonometric animation like this? Thanks

Please I need help with my proc.

> RK2skritt:=proc(FR::procedure,xo,yo,vxo,vyo,h) local x, y, vx, vy, r, kx1, kx2, kx3, kx4, ky1, ky2, ky3, ky4, l x1, l x2, l x3, l x4, ly1, ly2, ly3, ly4, tmp; r:=sqrt(xo*xo+yo*yo); tmp:=-(h*FR(r))/(r); lx1:=h*vxo; ly1:=h*vyo; kx1:=tmp*xo; ky1:=tmp*yo; lx2:=h*(vxo+0.5*kx1); ly2:=h*(vyo+0.5 ky1); r:=sqrt((xo+0.5*lx1)^(2)+(yo+0.5*ly1)^(2)); tmp:=-(h*FR(r))/(r); kx2:=tmp*(xo+0.5*lx1); ky2:=tmp*(yo+0.5*ly1); lx3:=h*(vxo+0.5*kx2); ly3:=h*(vyo+0.5*ky2); r:=sqrt((xo+0.5*lx2)^(2)+(yo+0.5*ly2)^(2)); tmp:=-(h*FR(r))/(r); kx3:=tmp*(xo+0.5*lx2); ky3:=tmp*(yo+0.5*ly2); lx4:=h*(vxo+kx3); ly4:=h*(vyo+ky3); r:=sqrt((xo+lx3)^(2)+(yo+ly3)^(2)); tmp:=-(h*FR(r))/(r); kx4:=tmp*(xo+lx3); ky4:=tmp(yo+ly3); x:=xo+(lx1+2*lx2+2*lx3+lx4)/(6); y:=yo+(ly1+2*ly2+2*ly3+ly4)/(6); vx:=vxo+(kx1+2*kx2+2*kx3+kx4)/(6); vy:=vyo+(ky1+2*ky2+2*ky3+ky4)/(6); [x,y,vx,vy]; end proc: Error, invalid local declaration Please someone help me!!!!!!!!

in maple after fitting an equation (A + B.x) in the experimental data, how to assign to parameters A and B the adjusted values will be assigned in the equation?

Thanks

I have a surface described by the equation.

f := (x, y) -> 3*(1 - x)^2*exp(-x^2 - (y + 1)^2) - (2*x - 10*x^3 - 10*y^5)*exp(-x^2 - y^2) - 1/3*exp(-(x + 1)^2 - y^2)

with(plots):
plot3d(f(x,y),x=-3..3,y=-3..3,orientation=[-40,60,0])

If a drop of water were put at point f(0.2,1.4), how would we draw a line to follow it's path?  (is that the gradient curve?)

I expect the output for the following code to be a Matrix containing 720 and a Matrix containing 21. What I get is a Matrix containing x^6 and one containing 21, which shows an unexpected difference between the parameter evaluation rule for Matrix products and that for Matrix sums. Is there a good reason for this, or is it a bug? Do any other binary operators between rtables also show this weirdness? 

proc(n::nonnegint) 
local x, q, v:= <<x>>, s:= <<0>>, p:= <<1>>; 
   for x to n do p:= p.v; s:= s+v od; 
   p, s 
end proc(6);

Why does the following code return a vector containing x rather than containing 1 and how can I correct it?

proc() local x, v:= <x>; x:= 1; rtable_eval(v) end proc();

My actual goal is to have a Matrix with a parameter inside a procedure and to create different instantiations of the Matrix in a loop by changing the parameter. 

What is the problem with the initialconsiditions that caused this error?

"Error, (in dsolve/numeric/DAE/initial) missing initial conditions for the following: {Integer}"

Thanks.

Hi

I use Maple 2016.

The following command calculates semicircle perimeter, but it returns infinity.

`assuming`([int(sqrt(1+(diff(sqrt(R^2-(x-R)^2), x))^2), x = 0 .. 2*R)], [R > 0])

I am running a for loop which for some loop index may give error message but it is fine for me to go with next loop indexes. but maple stop after the error message. How can I make maple to run even after the error message. Thank you

Can I do something like it?