Maple 17 Questions and Posts

These are Posts and Questions associated with the product, Maple 17

hi.please help me .why rule ''''evalf'''' dose not work properly?and I should copy the former answer in a another 'evalf' rule a gain??can I use only one rule 'evalf' instad twice useage??

this program attached blow.thanks alot

evalf.mw

 

I have the following transition rate matrix. Upto m the matrix follows a patern and after that the patern changes but repetative structure is there upto n. In this problem m=4, n=9. How to write these entries using for loop in maple?

Matrix.pdf

why the the software can't solve the integral like ∫xdlnx?

Thanks in advance for your help.

hi.please help me for solve this integral

thanks aloterror.mw

(Phi)(x,theta,z,t):=-Phi[0].sin((m*Pi*(x-x[0]))/(L)).cos(n*theta)*cos((Pi.z)/(h)).sin(omega*t):n:=1:L:=2:m:=3:x[0]:=6:h:=1/(2);
print(`output redirected...`); # input placeholder
1
-
2
evalf(Int(diff(Phi(x, theta, z, t), theta, theta), z = -(1/2)*h .. (1/2)*h));
print(`output redirected...`); # input placeholder
/0.2500000000
| / / /3 \\
| |(Phi[0]) . |sin|- Pi (x - 6)|| . (cos(theta))
| \ \ \2 //
/-0.2500000000

\
cos(2 Pi z)| . (sin(omega t)) dz
/

 

Hi, I am learning GR and new to maple. I played with maple's DifferentialGeometry package.

While duplicated the maple GR example, in DGsetup([t, r, theta, phi], M1, verbose),

I got the error,

DGsetup([t, r, theta, phi], M1, verbose)

The following coordinates have been protected:
[t, r, theta, phi]


The following vector fields have been defined and protected:


Error, (in Typeset:-Tdisplay[true]) `alias` does not evaluate to a module

 

 I don't know how to fix it. Can anyone kindly help?

Thanks

 

 

 

 

 

 

 

Hi, I have 10 equation system and 10 unknown variables. I however, want to reduce the equations to 2 with two unknowns. I'm wondering how this could be done in maple. The variables are Y, q, yd, y*, yx, H,  pd, w, P and Pv. I intend to solve the equations for Y and w. 

Thanks in advance for your help. The maple file has also been attached.

 

Regards

Small_Open_Economy.mw

 

hello everyone. I have an undergradute project i'm currently working on and I'm stuck where I have to use the Differential Transforms Method to solve a problem with boundary conditions at infinity


restart;

Digits := 5;

F[0] := 0; F[1] := 0; F[2] := (1/2)*A; T[0] := 1; T[1] := B; M := 2; S := 1;

for k from 0 to 10 do F[k+3] := (2*(sum((r+1)*F[r+1]*(k+1-r)*F[k+1-r], r = 0 .. k))-T[k]-3*(sum((k+1-r)*(k+2-r)*F[r]*F[k+2-r], r = 0 .. k))-M*(k+1)*F[k+1])*factorial(k)/factorial(k+3);

T[k+2] := (-3*(sum((k+1-r)*F[r]*T[k+1-r], r = 0 .. k))-S*T[k])*factorial(k)/factorial(k+2)

end do; f := 0; t := 0;

for k from 0 to 10 do

f := f+F[k]*x^k;

t := t+T[k]*x^k end do;

print(f);
print(t);

but the problem is that i cant seem to evaluate

or higer diagonal pade-approximant. any help will be greatly appreciated. thank you.

Hi .please help me for solve this nonlinear equations , that attch below

thanks a lots.....

 

 

dsys3 := {8*(diff(f2(x), x, x, x, x))+9*(diff(f2(x), x, x))+10*f2(x)+11*(diff(f1(x), x, x, x))+12*(diff(f1(x), x))+13*(diff(f3(x), x, x))+14*f3(x)+f3(x)*f3(x)+(diff(f3(x), x))*(diff(f3(x), x))+(diff(f3(x), x, x))*f3(x) = 0, 16*(diff(f3(x), x, x, x, x))+18*(diff(f3(x), x, x))+19*(diff(f3(x), x, x))+22*(diff(f1(x), x))+23*(diff(f1(x), x))+24*(diff(f2(x), x, x))+25*f2(x)+26*f2(x)+27*f3(x)+29*f3(x) = 0, 2*(diff(f1(x), x, x))+3*(diff(f2(x), x, x, x))+4*(diff(f2(x), x))+6*(diff(f3(x), x))+7*f1(x)+(diff(f3(x), x, x))*(diff(f3(x), x))+(diff(f3(x), x))*f3(x)+(diff(f3(x), x))*f3(x) = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)12.mw

Hello everyone,

How can unknown values in boundary conditions of system ODEs be determine by using shooting technique. See the work sheet here Sht.mw

Thanks.

Hi,

I am experiencing some problems with the plot3d command. It does not return a nice 3d-graph, but just a bunch of dashed lines. I have entered something like this:

with(plots):

plot3d(x*y,x=-1..1,y=-1..1);

I am using maple 17 on windows 7

Thanks in advance.

Hi,

Might be a silly question but how do I preview/display the expression I type instead of maple automatically solving it.

For example, If I type in '1+1' in wolfram, it displays 'input: '1+1'' right below it.

I would like to do this in maple but cant seem to figure it out.

Hopefully that made some sense.

 

Cheers and thanks! 

I just wanted to ask whether one can get out of the expansion

(D-f(x))@@2(g)(x)=(D@@2)(g)(x)-(D((f(x))(g)))(x)-((f(x))(D(g)-(f(x))(g)))(x)

the intended result

(D@@2)(g)(x)-D(f*g)(x)-f(x)*D(g)(x)+f(x)^2*g(x)

I'm just using it this way:

with(QDifferenceEquations):
QPochhammer(-1,5,10)

Error, (in QDifferenceEquations:-QPochhammer) wrong type of arguments

Am I doing something terribly wrong???

Dear Colleges

I have a problem with the following code. As you can see, procedure Q1 converges but I couldn't get the resutls from Q2.

I would be most grateful if you could help me on this problem.

 

Sincerely yours

Amir

 

restart;

Eq1:=diff(f(x),x$3)+diff(f(x),x$2)*f(x)+b^2*sqrt(2*reynolds)*diff(diff(f(x),x$2)^2*x^2,x$1);
Eq2:=diff(g(x),x$3)+diff(g(x),x$2)*g(x)+c*a^2*sqrt(2*reynolds)*diff(diff(g(x),x$2)^2*x,x$1);
eq1:=isolate(Eq1,diff(f(x),x,x,x));
eq2:=subs(g=f,isolate(Eq2,diff(g(x),x,x,x)));
EQ:=diff(f(x),x,x,x)=piecewise(x<c*0.1,rhs(eq1),rhs(eq2));
Eq11:=diff(theta(x),x$2)+pr*diff(theta(x),x$1)*f(x)+pr/prt*b^2*sqrt(2*reynolds)*diff(diff(f(x),x$2)*diff(theta(x),x$1)*x^2,x$1);
Eq22:=diff(g(x),x$2)+pr*diff(g(x),x$1)*f(x)+pr/prt*a^2*c*sqrt(2*reynolds)*diff(diff(f(x),x$2)*diff(g(x),x$1)*x^1,x$1);
eq11:=isolate(Eq11,diff(theta(x),x,x));
eq22:=subs(g=theta,isolate(Eq22,diff(g(x),x,x)));
EQT:=diff(theta(x),x,x)=piecewise(x<c*0.1,rhs(eq11),rhs(eq22));
EQT1a:=eval(EQT,EQ):
EQT2:=eval(EQT1a,{f(x)=G0(x),diff(f(x),x)=G1(x),diff(f(x),x,x)=G2(x)}):
bd:=c;
a:=0.13:
b:=0.41:
pr:=1;
prt:=0.86;
reynolds:=12734151.135786774055543653356602;     #10^6;   #1.125*10^8:

c:=88.419896050808975395120916434619:
;
Q:=proc(pp2) local res,F0,F1,F2;
print(pp2);
if not type(pp2,numeric) then return 'procname(_passed)' end if:
res:=dsolve({EQ,f(0)=0,D(f)(0)=0,(D@@2)(f)(0)=pp2},numeric,output=listprocedure);
F0,F1,F2:=op(subs(subs(res),[f(x),diff(f(x),x),diff(f(x),x,x)])):
F1(bd)-1;
end proc;
fsolve(Q(pp2)=0,pp2=(0..1002));
se:=%;
res2:=dsolve({EQ,f(0)=0,D(f)(0)=0,(D@@2)(f)(0)=se},numeric,output=listprocedure):
G0,G1,G2:=op(subs(subs(res2),[f(x),diff(f(x),x),diff(f(x),x,x)])):
plots:-odeplot(res2,[seq([x,diff(f(x),[x$i])],i=1..1)],0..c);



Q2:=proc(rr2) local solT,T0,T1;
print(rr2);
if not type(rr2,numeric) then return 'procname(_passed)' end if:
solT:=dsolve({EQT2,theta(0)=1,D(theta)(0)=-rr2},numeric,known=[G0,G1,G2],output=listprocedure):
T0,T1:=op(subs(subs(res),[theta(x),diff(theta(x),x)])):
T0(bd);
end proc;
fsolve(Q2(rr2)=0,rr2=(0..100));


shib:=%;
sol:=dsolve({EQT2,theta(0)=1,D(theta)(0)=-shib},numeric,known=[G0,G1,G2],output=listprocedure):
plots:-odeplot(sol,[x,theta(x)],0..c);
#fsolve(Q2(pp3)=0,pp3=-2..2):

Amir

I am facing a problem in maple 17. When I am launching maple 17, each time  the evaluating process is automatic started in server 1 without any input. I could not interrupt this evaluating process. For removing this problem I reinstall maple 17, but the problem is exist.  Please inform me soon how remove this problem ?

First 22 23 24 25 26 27 28 Last Page 24 of 58