Hi every one,

I am facing a situation here. Whenever, I sub n = 0 or 1, the following ode has a solution. 

But when I adopt other values except the above ones then the maple gives me the error

"division by zero". 

eq11:=(diff(((abs(diff(f(eta),eta,eta)))^(n-1)*diff(f(eta),eta,eta)),eta)

+n/(n+1)*f(eta)*diff(f(eta),eta,eta))=0;

resx:=dsolve({eq11,f(0)=1,D(f)(0)=0,D(f)(12)=0},numeric,method = bvp[midrich],

> f := a+D(x); 

> convert(f(t), diff);

Is there a way to tell Maple that "a" is a constant, so (t) won't attach to a? thanks, kyle

Hi there, I'm trying to pass some calculation to Matlab then retrive the result by getvar(), but seems something is getting wrong.

my matlab function is very simple:

function out =maple_fun(a, b)

out = a + b;

return

this function is saved as maple_fun.m.

in Maple,

restart:
with(Matlab):
openlink():
DirectoryLocation := " cd  '/Users/Kyle/Documents/MATLAB/'  ":

sol := dsolve(t*diff(phi(x),x$2)-x*diff(phi(x),x)+n*phi(x),phi(x));
phi := unapply(rhs(sol), x);
BC := [phi(0)=0,phi(1)=0];
with(linalg):
Ccoef := genmatrix(BC, [_C1,_C2]);
CharEqn := det(Ccoef) = 0;

restart;
sol := dsolve(t*diff(phi(x,t,n),x$2)-x*diff(phi(x,t,n),x)+n*phi(x,t,n),phi(x,t,n));
phi := unapply(rhs(sol), x);
BC := [phi(0,0,0)=0,phi(1,1,1)=0];
with(linalg):
Ccoef := genmatrix(BC, [_C1,_C2]);
CharEqn := det(Ccoef) = 0;

how can i specify (range and domain) of a function in the softwere,please answere if you are able to.

Division by zero, overflow error

A recent question brought to me some old applications in the application center, loading those applications in newer Maple versions shows critical information missing.

This prompted me to install MapleV on my computer to properly view these applications however after installing, sadly we get a division by zero, overflow error message. 

Since probably no has done this in this day and age, this is a real stab in the...

Hi,

I'm currently working on chemical process thermal exchange and particularly on the solving of the heat equation using a time dependant boundary condition.

Briefly, the process consists in two layers of different materials (M1 and M2, thickness L1 and L2). The bottom part of the material M1 (z=0) is cooled down from Ti to Tf with the function T(0,t)=Ti-R*t (R is the cooling rate in °C.min-1) until T(0,t)=Tf. Here the equilibrium is reached in t=(Tf-Ti...

Useing the following procedure I'd like to collect a set of  roots in a list of lists, so they can be manipulated and presented in various plot options. Some of which could be very interesting (to me only perhaps).

One option is to create a loop which produces values of newton roots x0 for a given function. For example:              f:=x->(x^7)-5

for x from -1 to 1 by 0.05 do

Dear Maple Users

I am pretty new to programming in Maple, so maybe this is not a very difficult question to answer. I want to write some code, which generates a table to be displayed. In every cell I want some content, which is generated by some code as well. How can it be done, and if there are several ways, which one is the best?

Regards

Erik

I guess there is a limit bug here somewhere.  With MultiSeries things are OK...

L:=(x*BesselI(k-1, x)+x*BesselI(k+1, x))*BesselI(-k, x);
L := (x BesselI(k - 1, x) + x BesselI(k + 1, x)) BesselI(-k, x)
simplify(limit(L,x=0)) assuming k>0;
                               0
with(MultiSeries):
simplify(limit(L,x=0)) assuming k>0;
                          2 sin(Pi k)
                          -----------
                              Pi     

BesselLimit.mw

Hello,

I'm solving a system of 3 second order ODE's for 3 unknown functions, using dsolve. But in the output, I receive various sets of solutions which contain only 3 constants of integration. A friend solved the same system but in Mathematica, and received more constants of integration. 

I guess this has to do with the issue that Maple omits additive constants of integration almost everywhere. 

Is there some solution? i.e. some way of using...

How do I uninstall everything associated with Maple Prime?  I don't know how I got it and don't want it!  I don't know what it does and don't need it.  Operating with new iMac running OS 10.7.4.  It's 3:36 AM and I've been trying to eliminate it fot the last several hours!!  Help now please before I go nuts!!

Charles Stiles

We need to compute some integral, containing product of some "regular" function (e.g. gaussian)
and Dirac delta-function. The problem is that argument of delta-function is itself a function.

A simplified example is

d1:=s->int( g(x,5,0.025)*Dirac( s-f1(x) ) , x=0..infinity);

By trials and errors, we figured out that Maple runs such calculation successfully
when this function is very simple and fails in more complex cases, i.e., integration  returns zero.

I am looking for a numerical scheme which can easly handle such type of problem.

I already had the analytical results of this problem, now I want to know how to treat this problem numerically.

Any idea, then please let me know.

Thanks

I am trying to solve a transcendetal equation that has no analytical solution (AFAIK).
My question is how to find different branches of solutions with Maple. So far, Maple
only returns one particular solution.

To fix ideas, my equation is of the type,
x-a*(1-x)* exp(bx)+c=0
For suitable values of a and b this equation seems to have, at least, two solutions in
the real axis and I would expect it to have a numerable infinite of solutions in the
complex...