## why evalf dosnt work????...

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

## How to write that matrix in Maple?...

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...

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

## problem with solve 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
/

## Get Error in "DGsetup([t, r, theta, phi], M1, verb...

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

## Solving system of equations...

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

## Differential Transforms Method...

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.

## Error, (in dsolve/numeric/bvp/convertsys) unable t...

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

## How to determine unknown value in boundary conditi...

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.

## plot3d problems...

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

## How to preview input like wolfram ...

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.

Cheers and thanks!

## expand differential operators...

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)

## QPochhammer wrong type of arguments...

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???

## ODE- convergence problem...

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