Hi, I am trying to solve a maximization problem using the following Lagrangean: L:=v*K*a^gamma*k^(1-gamma)-1/2*a^2-1/2*k^2 -1/2*r1*v^2*sigma^2-1/2*r2*(1-v)^2*sigma^2 -mu*(v*K*gamma*a^gamma*k^(1-gamma)/a-a) -lambda*((1-v)*K*(1-gamma)*a^gamma*k^(1-gamma)/k-k) I first computed the partial derivatives for a, k, v, mu and lambda (La:=diff(L,a)=0, Lk:=diff(L,k)=0, ...). (K, gamma, r1, r2 and sigma are constants) Then I applied solve ({La,Lk,...},{a,k,v,lambda,mu}) and did not get any result (or message) back. I already tried to solve this problem by step-by-step calculations, but was not successful in getting solutions for a, k and v (stated in terms of gamma, r and sigma). When I make the problem simpler (by setting r1:=0 and r2:=0) I get a solution.

Can Maple be used to create solids of revolution by rotating intersecting plane curves or a single curve bounded by certain intervals about lines or the cartesian axes ? If yes can someone help with the steps to accomplish this ?

hi how i get all the roots? > fsolve('P2(x)=0',x=0..1); 0.7653758275 > fsolve('P2(x)=0',x=0..0.6); 0.1740920943 > fsolve(P2); 4.274188187 i need the smallest root and P2(x) is numerical solved it's not continius

pls check the attached worksheet ... any reason for that behaviour?

PS: not attached due to some oddity of the forum software,

"Error attaching file 102_strangeRoot.mws: exceeds maximum file size"

so take it from http://www.axelvogt.de/maplekram/strangeRoot.mws

(or *.zip if perfered, it is below 4 kB anyway ...)

I'm trying to figure out how to fit the data below to a curve of the form:
a+b*sin[(pi*x)/6-c]
How do can I do this in Maple 10?
1,43.2
2,47.2
3,53.7
4,60.9
5,70.5
6,82.1
7,88.0
8,86.0
9,77.4
10,66.0
11,51.5
12,44.1

Is it possible with Maple 10 to change a document so that in the middle, there are two columns? One might imagine having some text describing how data was gathered where the text is spread over the usual six inches or so. Reading on down the document, the style might change to two columns: a graph of the data in the right column and comments on the data in the left. After this, the document could revert to the normal six inch, full line text.
Did I not see something such as this in a demo at the Maple Summer Conference?

I want to display a 3dplot of a current

iL_RMS_3dplot_Ue_const:=plot3d(iL_RMS_func_Ue_const, 0..2, 0..1, axes=normal, title="RMS current through the inductor with Ue constant and Ua variable", labels=["Ia_IaLGmax", "Ua_Ue", "iL_RMS_Ue_const"]);

together with an output curve. The output curve (which is two dimensional) should have as the z value in a 3d plot the value of iL_RMS_3dplot_Ue_const at that place. I have uploaded an example to http://www-users.rwth-aachen.de/Kai.Schueler/3dplot_example.jpg (the blue curve is the output curve). Is it possible, without recalculating the z value with a function to call, to extract it from the iL_RMS_3dplot_Ue_const plot data structure with an op command or something similar. The 3dplot command for the output curve should look somehow like this:

I don't understand, why maple has problems with the following command:

Output_Curve_test:=plot3d([0,0,iL_RMS_func_Ue_const(Iatest,U_test)],Iatest=0..2,U_test=0..1);

It seems that Maple doesn't call the funciton iL_RMS_func_Ue_const with the function variables Iatest and Utest. Instead it just uses the variable names Iatest and Utest. I get the error message:

Error, (in iL_RMS_AVG) cannot determine if this expression is true or false: 46.87500000*Iatest <= .1250000000*(1500-1500*U_test)*U_test

How do I get Maple to use the actual Values of Iatest und U_test? Below are the functions which are related to the problem.

Hi, All,
Please look at this example. I am using Newton's method to solve for the equation x^2-a=0.
> restart:
> Digits := 30:
> f :=x-> x^2-a:
> df :=x-> 2*x:
> for a from 1 to 6 do
> x0 := trunc(sqrt(a)):
> for j from 1 to 5 do
> x1 := x0 - f(x0)/df(x0);
> x0 := x1;
> rsd := f(x0);
> end:
> printf(`%3.10f\t %3.10e\n`, x0,rsd):
> end:
Now, I want to plot the solution x and the residue rsd against a in 2 separate figures, how to write the code?
David

Hello everyone, I was wondering what is the syntax to put a boundary condition like http://stw.ryerson.ca/~m4yip/bc.png in the dsolve command? I can do simple ones like v(0) = 0, but I don't know how to input boundary conditions involving a derivative ... Onto my second question: Is it possible to include a delta function within a Maple differential equation? I'm interested in the differential equation of EI v'''' + P v'' = Q

I've been exploring the link between Maple and Matlab, and there seems to be no documentation on how to use Maple to call commands (or .m files) in Matlab. On the other hand, it doesn't state that this is not possible. Has anyone else looked into this?

I am trying to create a procedure like this:
myproc := proc(x)
Digits := x;
sqrt(x)
end proc;
and then I want to plot the results like this:
plot(myproc(x), x=2..20);
However, I get an error: "Invalid assignment to Digits".
The following works fine, however:
seq(myproc(x), x=2..20);
How can I solve this?
Thanks.

How do I get maple to do polynomial long division?

It's obvious how to use Maple 10 for creating a document, but has anyone used Maple 10 for presentations? I'm going to present in a seminar in a couple of weeks and have some Maple demonstrations and graphics to show the group. I'd be interested in hearing tips, suggestions, experiences and ideas for presenting mathematical content to a group from within Maple. (I'm familair of course, with Powerpoint and LaTeX slides, and the LaTeX Prosper package. I'd use LaTeX/Prosper except that then I have to embed the Maple in the presentation. It would be more convenient if I can use Maple solely.)

Dear All,
First of all thank alec and Prof Doug. Please refer to post http://beta.mapleprimes.com/node/362#comment-729
My problem is to solve the following ODE with constrain given by the Eq:
> DEq := diff(f(x),x) = c1*sqrt( g(f(x)) - g(f(x=x1)) );
> Eq := ug - f(x=0) = c2*sqrt( g(f(x=0)) - g(f(x=x1)) );
> g(f(x)) := exp(f(x)-xn) +exp(-f(x)) +(1-exp(-xn))*f(x);
where ug is the independent variable and others are parameters
> c1 :=-257753830.552993398770016153893;
> c2 :=1.40376701562707174099362643812;
> xn := 35.4196679725860384893589519100;
> x0 := 0;
> ug := 40;
I have tried to followed alec's suggestion but the solution sometime still cannot be obtained. The independent variable ug has to be varied from -50 or less to + 100 or more. The intial guess in solving f0 (f0 := fsolve(g('fx0'(f0))-gfx0(f0),f0=30);) has to be input manually, this is a big drawback.