I want to compute the following solve:

C:=A, B, F, H(It has 4 unknowns):

V:=R,Y,E,I (It has 4 equations):

solve({V},{C});

By hand calculations I just found that the solve result is real ,but in maple there is no solution and no error (the maple ignores to solve this equation !!!) 

can anybody help me to find the solve soltuion ?

Thank you 

H. Kabir

The font size in the output pane of the standard interface, interactive debugger is really, really small and I can't figure out how to increase it. (Other than drop my screen resolution from its native 2560x1440, which makes it bigger but "fuzzier" - not a huge improvement for tired eyes!!)

Hi dear friends

How do I improve my ability to manipulate code?

Thanks

_________________________________________________________________________________________

restart:

Digits := 30: m := 20: p0 := 8.: g := -0.3: nu := 0.3: a := 1:

w := sum(b[n]*r^(1.02*n), n = 1 .. m):

ODE := r^3*(g*r+a)^3*(diff(w, r, r, r, r))+2*(diff(w, r, r, r))*r^2*(4*g*r+a)*(g*r+a)^2+(diff(w, r, r))*r*(3*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2+11*g^3*r^3+3*a^2*g*nu*r+15*a*g^2*r^2+3*a^2*g*r-a^3)+(diff(w, r))*(6*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2-2*g^3*r^3-3*a*g^2*r^2+a^3):

for P from p0 by 0.1e-1 to p0+1 do
W := simplify(subs(solve({seq(evalf(subs(r = j/(m+2), ODE)), j = 2 .. m-1), subs(r = 1, diff(w, r))}, {seq(b[j], j = 1 .. m)}), w/b[m])):
F := (int(((diff(W, `$`(r, 2))+(diff(W, r))/r)^2-(2*(1-nu))*(diff(W, `$`(r, 2)))*(diff(W, r))/r)*r*(1+g*r/a)^3, r = 0 .. a))/(int((diff(W, r))^2*r, r = 0 .. a)):
delta[trunc(100*(P-p0))] := abs(F-P) end do:

plot(W, r = 0 .. a);

Min := min(seq(delta[j], j = 0 .. 100)):

for j from 0 to 100 do if delta[j] = Min then print(p0+j/(100.)) else  end if end do;

Hi all,

I have this matrix for example:

A := <<1,2>|<3,4>>;

How to determine: 2^A

Thank you

plz help me, how do i solve singular ODEs of lane Emden type equation for homotopy analysis method in maple? there is arising an arror, invalid fraction

To begin with my code is:

with(Statistics);

N:=RandomVariable(Normal(5,2))

for i from 1 by 1 to 100 do

OS[i]:=Sample(N,50)

end do

for i from 1 by 1 to 100 do

 print(i)

  from j from 1 by 1 to 1000 do

    BS[i,j]:=Bootstrap('Mean',OS[i],replications=1)

  end do

end do

My goal is to produce as an output in the software's screen the values of 1000 bootstrap means for each of 100 original samples in order to export the results to Excel. Based on my basic knowledge of programming I thought that I can do that through the use of the vectors. Unfortunately though when I trigger the operation the software seems to be evaluating (with the homonymous indication in the bottom left corner of the screen) it for an indefinite period of time, obliging me to stop it at some point. However, let me point out that when I do the same process for i and j 4 and 10 for example then I have results, therefore presumably my code is correct. I know that I can do the same operation in Excel and I have already created a respective macro, however I wanted to use Maple in order to take advantage of its fast operation of calculations, which is pointless in case that I have to wait more time for the results than I have already waited to be eventually produced.

Any thought and idea would be appreciated, thank you very much in advance!

As i  want to build a model of engine an want to supply compressed air for same ii prepared a engine model but now i am not geeting how to prepare model flow of compressed air and its expansion in enige and force crated on engine head.

How do I solve canonical diffusion mathematical equation in Maple?

http://www.maplesoft.com/applications/view.aspx?SID=5083&view=html&L=F

Error, (in hanoi) cannot determine if this expression is true or false: 0 < n

N:=n:moves:=[];
hanoi:=proc(N,source,dest,aux) global moves;
if N>0 then
hanoi(N-1,source,aux,dest):
moves:=[moves[],[N,source,dest]]:
hanoi(N-1,aux,dest,source);
end if:
end proc:
sol:=hanoi(n,1,3,2);

There is million items in list, 

what is the difference between seq and for loop ?

is seq faster than for loop?

it is very slow when running code below, how to speed up this part of code?

https://gist.github.com/LovelyYanki/c0b61fbb9d5954b34e03#file-gistfile1-txt

HSKeyIn := Table();
for ij from 1 to 1685159 do
s := solve([mm[ij][1,1]=1,mm[ij][1,2]=1,mm[ij][1,3]=1],[a,b,c]):
if nops(s) > 0 then
if (rhs(s[1][1]) = 1 or rhs(s[1][1]) = 0) and (rhs(s[1][2]) = 1 or rhs(s[1][2]) = 0) and (rhs(s[1][3]) = 1 or rhs(s[1][3]) = 0) then
#print(lhs(indices(T3, pairs)[ij])):
#print("***"):
h := HilbertSeries([mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]], {a,b,c}, z):
#print(h):
if not assigned(HSkeyIn[h]) then
if mod(ij, 100) = 0 then
print(ij):
end if:
HSkeyIn[h] := [[mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]]]:
else
if mod(ij, 100) = 0 then
print(ij):
end if:
HSkeyIn[h] := [op(HSkeyIn[h]), [mm[ij][1,1],mm[ij][1,2],mm[ij][1,3]]]:
end if:

end if:
end if:
od:

mm := [[a, 1], [a, b]];
mmseq := {seq(mm[ij],ij=1..2)};

how to filter above in seq when any items in the [] has 1

criteria is below

if mm[ij][1] = 1 or mm[ij][2] = 1 or mm[ij][3] = 1 then
Dummy := 0:
else
mm2 := [op(mm2), mm[ij]]:
end if

People im with error to show a plot, this is a wor for universty about mechanic materials, and i'm with this error... do not show the plot... before, i've a error because the variable x has values, and i've to unassign to have this.. without unassign, has error in domin of the plot... i dont know how to do this ! And ive to send the work today. 

Someone can help me about this?


Value of variables:

plot algoritm...

Thanks !

in page 137 of an introduction to groebner bases

how to eliminate the redundant solution (y^2-x*z, 0, -x^2+y*w)

from 3 of them?
eliminate({y,y^2-x*z,-z}, {x, y, z, w});
eliminate({-x,0,y}, {x, y, z});
eliminate({w,-x^2+y*w,-x}, {x, y, z, w});

Please help me to solve this integration

restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity);

restart; with(LinearAlgebra); int(exp(-(ln(z/(snr*B^2))+4*sigma^2)^2/(32*sigma^2))*eta^2*(y/z)^((1/2)*eta^2-1)/(z*sqrt(32*Pi*sigma^2)*(2*sqrt(y*z))*(2*A[o]^(eta^2))), z, z = y/A[o]^2 .. infinity)

I have created and saved a MAPLE module in an .mla archive. The module contains three procedures A, B, C, where

A calls, B and C.  

Once the module library has been loaded, A acccepts inputs and generates outputs.

Is it possible to create a MAPLE player worksheet which calls the module and share it with a Maple Player (only) user, so that they can then supply the inputs and observe the outputs from A using the Maple Player programme components?

Can anyone help?

MRB