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 ?

Hallo:

I am trying to find the inverse fourier transform given a vector (of real) as input. I get the following error message:

Error, (in DiscreteTransforms:-InverseFourierTransform) entries of input data must all evaluate to complex floats, got '[.7644453211]'

I am unable to understand this message. I would much appreciate help to resolve this issue.

Thanks

Good day everyone,

please how can one solve this pde in terms of Bessel function or any other analytic solution with the plot.

See the file ID.mw

Thanks.

While I was using DEplot3d, I want to draw a space curve, the equation of which is a list of functions of the vars in DEs. I put them in "scene", but Maple told me this:

Error, (in DEtools/DEplot) Invalid scene; must be list of vars: scene = ...

Now I know that I cannot put functions in "scene", but I am wandering if there are some other way so I can draw this kind of space curve.

Thank you very much!

My task is to develop a mathematical model of time-variant temperature
distribution in a bar with Maple. 
The bar is made of aluminum. The length is 202 mm and the diameter
is 8mm.
Heat is supplied from the kitchen lighter where the flame burns on butane.
The height of the flame is about 4 mm.
Your model should be able to answer how long the time it takes to reach
a certain value of temperature at the distance of 10 cm from the heat source.

anyone can help to answer this? im totally new to maple.. hopefully some1 can help me to answer this..

Good day everyone,

Please I do get numerical output/values in this solution U1.mw

Best regards

Does anyone know how to using some softward convert ?  for example  I have maple code , but I want to using mathematica code .I need fast way.