I have a surfdata-plot which - by interpolation - goes through different circles with different inclinations with respect to each other. The surface is color-coded with zhue from violet to red. In order to see the valleys behind peaks I would like to have a filled 2D-plot of my zhue-surfdata colors with contours wherever the surface passes z-coordinate integers located somewhere below the surface (e.g. at z=-7). Unfortunately, I haven't found examples of contourplots for surfdata plots so far.

Download Worksheet_Surfdata-contourplot.mwWorksheet_Surfdata-contourplot.mw

Hi all
how can I solve an equation? Pleas, help!

I have a differential equation with one unknown.
I first solve the differential equation, and then use the boundary condition to find his unknown, however, get the error:
Error, (in sol_2) parameter 'X' must be assigned a numeric value before obtaining a solution.

How can I find X?

 F:=1:
 l:=1:
 a:=0.1:
 n:=4.5:

 B:=1.47*10^(-11):
 J_nc:=evalf(2*a*int(y^(1/n+1),y=0..a/2)):
 M_f:=z->piecewise(`and`(z>=0,z<=l),F*l,`and`(z>l,z<=2*l),F*l-F*(z-l),`and`(z>2*l,z<=3*l),F*l-F*(z-l)+F*(z-2*l));
 M_ed:=z->z-3*l:
 M_x:=z->M_f(z)+X*M_ed(z);
 eq := diff(V(z), `$`(z, 2)) = B*(abs(M_x(z))/J_nc)^n*signum:cond := V(0) = 0, (D(V))(0) = 0;
 sol_2 := dsolve({cond, eq}, numeric);
 sol_3:=solve({sol_2(3*l)=0},{X});

Please see attached document.

The 2 bits of data I am working with are a Constraint Function and Objective Function:

I need some pointers on how to do this. I keep getting solutions of

and I am ont receving an x value for some reason and leaving me unable to continue with the problem.

If anyone can send me some code I would be grateful. If anyone would like to send the correct full working code so I can see how to do this and review it I would be grateful.

 Code I had so far was:

When trying to plot "Example1" from this page www.rhitt.com/courses/227/su99/Maple/html/surfaces2.html I get an error message that Maple cannot make a plot structure from an object with the transformation proc( (x,y) -> [x,y,-.3] end proc) . How can I plot this example?

Helo Dears, I need help to generate a vector whose detail given below

`for `K:=1 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]]]]

`similarly for `K:=2 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[2,0]],[psi[2,1]],[psi[2,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]],[psi[2,0]],[psi[2,1]],[psi[2,2]],[psi[2,3]]]]

Simply I need to generate a general vector for arbitrary values of k and M as discus in above.

New_Wavelet_with_Operational_Matrix.mw

We're starting on indefinite integrals in my 1st year calculus class.

A quick example would be int(sin(x), x);=-cos(x)+C

Maple doesn't add the +C on the end of it's solution. Can someone explain or point me to a resource? I've tried searching but I can't find an answer. 

I have the following functions: and 

I want to differentiate e(a,A) wrt A and I keep getting the following error

Error, invalid input: ln expects its 1st argument, x, to be of type algebraic, but received [gamma/(c(A)^(sigma*phi)-1+gamma)].

Not sure how to proceed.

Thanks in advance for your help.

Hello there,

I have, of late, experienced some funny output when using the DifferentialGeometry package.

The command DGsetup retuns vector fields and 1-forms with stars enclosed in back primes attached to them. The same predicament occurs as a result of an operation (evalDG) but only when the vector or the form is not multiplied by a scalar.

This phenomenon does not occur with Maple 17!

Has anyone experienced such a phenomenon and how has they fixed the issue (providing it is an issue).

Many thanks in advance.

Kind regards,

Freddy Baudine

Here's the file:

DiffGeo_potential_issue.mw

hi friends

i have a problem in maple with an error

dx := diff(x(t), t, t) = -G*Mz*x(t)/(x(t)^2+y(t)^2)^(3/2):

dy := diff(y(t), t, t) = -G*Mz*y(t)/(x(t)^2+y(t)^2)^(3/2):

 G := 6.67*10^(-11); Mz := 6*10^24:
 IniC := x(0) = 7*10^6, (D(x))(0) = 0, y(0) = 0, (D(y))(0) = 9*10^3:
 Digits := 15:
 Ns := dsolve({dx, dy, IniC}, {x(t), y(t)}, numeric):
 dsnumsort(Ns(0), [x, y]);

>for i from 0 to 400 do T := 40*i; NsT := Ns(T); X[i] := rhs(NsT[C1]); Vx[i] := rhs(NsT[V1]); Y[i] := rhs(NsT[C2]); Vy[i] := rhs(NsT[V2]); MofI[i] := X[i]*Vy[i]-Y[i]*Vx[i] end do:

> with(plots);
> p1 := polarplot(6378*10^3, phi = 0 .. 2*Pi);
> p2 := plot([seq([X[i], Y[i]], i = 0 .. 327)], thickness = 2);

display({p1, p2}, labels = ['x', 'y'], scaling = constrained);

but i see and I I can't draw PLOT:

Error, invalid input: rhs received [t = HFloat(0.0), x(t) = HFloat(1.0), diff(x(t), t) = HFloat(0.0), y(t) = HFloat(0.0), diff(y(t), t) = HFloat(1.0), z(t) = HFloat(0.75), diff(z(t), t) = HFloat(0.0)][C1], which is not valid for its 1st argument, expr

can you helpe me?

I'm a student using Maple 18.

I'm still having some navigation issues with the Maple text editor. I realize that being in "math" mode or "text" mode causes Maple to interpret an <enter> key press differently.

However when I paste in some text, I get unexpected behavior. 

My professor will assign a list of problems, and I will paste them into Maple so it looks similar to this

1.) solve math problem 1

2.) solve math problem 2

....

Then I go to problem one, put my cursor at the end of the line and push <enter> hoping to make some room between problem 1 and 2 for my work. When I do this, the cursor jumps to the end of the list I had just pasted into the document.

How can I insert new blank space between the listed problems? I've had the same experience pasting in from MS word or just plain notepad.

Right now the only solution I've found is to copy the listed problems into maple 1 at a time.

Thanks in advance for any suggestions.

Brent

Hello
There is some kind of strange behaviour with copying/pasting maple code from IDE to some text editor. Namely, the code collapses into one long line in text editor. Is it OK with Maple IDE?

I found this problem in a book about mathematica.
I want repeat it with maple.
But even Digits=1400, the exact value of s@x=70 is still different with the result(0.633319) on that book.
How to resolve this problem?
I want some help.

code:
restart;
Digits:=1400:
p:=convert(series(cos(x),x,201),polynom):
s:=add(op(i,p),i=[1,2,3,-1,-2,-3]);
add(evalf(coeff(p,x,i))*70^i,i=[200,198,196])+add(evalf(coeff(p,x,i))*70^i,i=[2])+add(evalf(coeff(p,x,i))*70^i,i=[4]);
#plot(s,x=0..40,view=[0..40,-2..2],numpoints=6000,gridlines,style=point);
simplify(eval(s,x=70));
evalf(%);

Dear all
I am using Physics[Vectors] package of Maple 17. I want to define an orthogonal curvilinear coordinates through alpha and beta independent variables. To define unit vectors of alpha and beta, I have to apply the derivatives of position vector r_ with respect to alpha and beta, respectively. Please help me to define the unit vectors in directions alpha and beta as derivatives of position vector r_ with respect to alpha and beta, respectively. Please see the below Code:

> restart;

> with(Physics[Vectors]);

> Setup(mathematicalnotation=true);

> r_:=X(alpha,beta)*_i+Y(alpha,beta)*_j+Z(alpha,beta)*_k;

The unit vectors in directions alpha and beta should be defined as:

_alpha=diff(r_,alpha)
_beta=diff(r_,beta)

It is worthwhile to mention that the following expression governs:

> diff(r_,alpha).diff(r_,beta)=0;

Please guide me.

Best wishes

Ali