This is the procedure for the question I posted on http://www.mapleprimes.com/questions/128600-Romberg-Approximation-Procedure-In-Maple?sq=128600 

Can anyone see where the problem is?

RombInt:=proc(f,a,b,m,n)
local j,k,R;
if type(m, numeric) then
   if not type(m,nonnegint) then ERROR("m must be a non-negative integer"); fi;fi;

Hi, I want my output to be column vectors, how do I use the printf statement?

Hello there,

When I was trying to solve a differential equation with dsolve, I've got such a error report:

> dsolve({Q(0) = F[Az], diff(Q(x), x) = q[2](x)});

Error, (in dsolve) found the following equations not depending on the unknowns
of the input system: {Q(0) = -(1/4)*l^2*q[0]-(1/4)*l^2*q[0]*exp(2)+F}

But if I only write:

> dsolve(diff(Q(x),x)=q[2](x));

then there is no problem. 

How to find the discontinuities of the function
g:=unapply(int((1-z*(cos(t)+cos(2*t)))/(1-2*z*(cos(t)+cos(2*t))+2*z^2*(1+cos(3*t))),t=0..Pi),z);

over the reals with Maple?

Hello,

I'm using Maple 7 and got a question:

I've got the function y

y:=x->(x)^0.5; ug:=0; og:=10;

and

xs:=2;

ys:=3;

How do I create a point with the coordinates [xs,ys] and plot it with the function y in one graphic?

I'm using plot(y(x),x=ug..og); to plot the function y.

Thanks for your help!

flash20001

hi,i want my output be a complex number,how do i use the printf command?for example  s=3+4*i be a output.

Hey, I need to write a procedure to calculate the romberg approximation by combining the trapezium rule approximations. The command will be RombInt(f(x),a,b,m,n), where a and b are the 2 limits and m is a number more or equal to 0 which indicates the number of iterations and n is a integer more or equal to m. Also given is that when m=0, this approximation can be called TrapRule(f(x),2^n,a,b). (I have got the procedure for this traprule, if that helps)

So basically...

Hey, I have made this procedure for the trapezium rule but I get an error when I perform this calculation.

Procedure is 

TrapRule:=proc(f,N,a,b)
local k,sum1,h,i:
h:=(b-a)/N:
sum1:=0:
for i from 1 by 1 to n-1 do
sum1:=sum1+f(a+i*h):
end do:
k:=(h/2)*(f(a)+2*sum+f(b)):
return (k):
end proc:

When I try calculating (TrapRule(exp(x),4,-3,1)); I get 

Error, (in TrapRule...

Simple question:

Solve y''(x)-y'(x)+2y(x)=10e^(-x)sin(x)

Mathematica answer:

which is correct.

Maple 15 can get the homogeneous solution (the first two terms) but the particular solution? Not as we know it.

y(x) = exp((1/2)*x)*sin((1/2)*sqrt(7)*x)*_C2+exp((1/2...

I am evaluating a difficult integral of an oscillating function.

f:=k->evalhf(piecewise(k=0,1,k>0,cos(k*(dz))*BK0(k*rho1)/BK0(k*a)));

evalf(Int(f, kmin..kmax, method = _d01akc,epsilon=eps,maxintervals=1000000))

with dz>0 and possibly very large
kmin is Pi/(2*dz)
the number of cycles is nc = ceil((1/2)*kl*dz/Pi) with kl = 91.563
kmax=  kmin+nc*2*Pi/dz

Each evaluation of this integral using a call to the NAG routine (_d01akc ) eats up 

I'm using the printf command and i want my output a decimal number to get it published

s1=0.5555552487 and f(s1)=0.2542365647 what should i do?

hi;

when run file , maple give wrong result, plz help me

my file is attached

thx

 EKM3-cylindrical_pan.mw

I have solved  the following BVP system in Maple

restart; e := 12; Sc := .1; Pr := .7; Gr := 1; M := .1; sys := diff(f(t), `$`(t, 3))+2*f(t)*(diff(f(t), `$`(t, 2)))-(diff(f(t), t))^2-M^2*(diff(f(t), t))+Gr*theta(t) = 0, diff(theta(t), `$`(t, 2))+2*Pr*f(t)*(diff(theta(t), t)) = 0, diff(phi(t), `$`(t, 2))+Sc*(2*f(t)*(diff(phi(t), t))-phi(t)) = 0; bcs := f(0) = .1, (D(f))(0) = 0, theta(0) = .5, phi(0) = 1, (D(f))(e) = 0, theta(e) = 0, phi(e) = 0; sol := dsolve({bcs, sys}, numeric); plots[odeplot...

For educational purposes, I would like to hide some code (procedures) in a maple file. I put these procedures in a subsection which I can collapse. But everytime you execute the worksheet, the subsection is automatically expanded and I have to collapse the subsection again.
Is there a way to avoid the automatic expansion of subsection while executing the worksheet?