restart

u := proc (x, y) options operator, arrow; u[x](x)*(y/delta(x)-2*y^2/delta(x)^2+y^3/delta(x)^3+B*H*(1-3*y^2/delta(x)^2+2*y^3/delta(x)^3)/delta(x))/(3*B*H/delta(x)+2) end proc

int(u(x, y)^2, y = 0 .. delta(x))

I tried the following partial integration in MAPLE. But Maple says:

Error, (in u[x]) too many levels of recursion.

Analytically it is possible to evaluate the integral in the above lomits. Where lies the problem?

I wrote the following MAPLE lines to determine the three constants a, b and c from three boundary conditions. Whatever a, b and c I get as result are correct. Now I want to substitute back the values of obtained a, b and c in the expression for Te to obtain the resulting temperature distribution. I used both algsubs(a, b, c, Te) and subs(a, b, c, Te) one by one to substitute a, b and c in the expression for Te and get the result. But both way it is showing...

Hi!

I have seven curves (expressions) R:=f(x) for different value A.

1. How to Plot all curves on one 3D spece plot?
2. Hot to get (fit) surface equation R:=f(x,A) from all 2D curve?

surf_eq.mw 

wzel

Hello,

I have four old Maple 10 Worksheets sitting around my harddrive. I'd like to view them again, but my Maple 10 license won't be activated again, they tell me to buy a new 13 or 14 one.

Does anyone still have a Maple 10 (or compatible) version and would convert four small worksheets into PDF for me?

Thanks!

I wrote the following Maple lines:

I wrote the following lines in MAPLE:

Sorry

 endv=100 not 5 (as shown below)

I have used dsolve to obtain the solution of the inserted two systems below. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is in approximating the first system solution to be used in the second one. I would appreciate any advice on how to  obtain accurate results for the second system

I have used dsolve to obtain the solution of the inserted two systems below. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is in approximating the first system solution to be used in the second one. I would appreciate any advice on how to  obtain accurate results for the second system

 restart; endv := 5; sys0 := diff(f(x), `$`(x, 2...

Daer all! I want to solve the following equation

> H(z)*(1-m*H(0)^2*(1+z)^3/H(z)^2)-2*e*sqrt(1-c^2)*sqrt(r)*H(0) = 0;

where H(0) denotes the value of H(z) at z=0 and other constants are

> m := 0.211: c := 0.80: r := 0.338: e := 1:

If I want to use the fsolve command as

> eq := z->H*(1-m*H(0)^2*(1+z)^3/H^2)-2*e*sqrt(1-c^2)*sqrt(r)*H(0) = 0;

> Y := z-> fsolve(eq(z), H=1):  

what should I imply instead of H(0)? 

Suppose I have the following expression:

u=a+bx

Now I want to determine a and b using the conditions:

at x=h, u=4;

at x=6, du/dx=4.

Can I do the same using MAPLE? If possible how can I do it. Where to look for it?

Thanks in advance.

Hi I teach in a high school. Basically I have found that if I ask a question like: Factorize fully (x⁴-y⁴), and use a MathContainer for both the question and the answer (to be filled by students), I can't find a way to see whether the students have fully factorized the expression. Suppose (x⁴-y⁴) is in MathContainer1 and MathContainer2 is to be filled in by students. Using the code:

if GetProperty('MathContainer1','expression'...

I used dsolve to obtain the solution of the attached two systems of ODEs, where the second system depends on the solution of the first one. The obtained results for the first system is accurate, where the obtaind solution of the second one is not accurate. I think the problem is using cubic spline approximation. I would appreciate any advice on how to  obtain accurate results for the second system

Dear guys! I have an equation as

> eq := z-> H*(1-m*h^2*H^(-2)*(1+z)^(3))-2*e*sqrt(1-c^2)*sqrt(r)*h=0:

where

> m := 0.211: h := 0.741: c := 0.80: r := 0.338: e := 1:

At this equation H is an unknown function of z, h is the value of H at z=0 and I want to obtain H(z) actually(for example the graph of H(z) for z=-1..5). So I did as following

> Y := z-> fsolve(eq(z), H=h):    **

plot(eval(H, H='Y(z)'), z=-1 .. 5);

hi

I need to plot the phase portrait of a hamiltonian of 4-dimensional phase-space, and I have an adiabatic invariant, i.e. one of the momenta is conserverd. 

I don't know how to replace a constant for the conserved momentum and plot the 2-dimensional phase space.

I am writiing 

sys := {diff(Y1(t), t) = -(diff(K, y1)), diff(Y2(t), t) = -(diff(K, y2)), diff(y1(t), t) = diff(K, Y1), diff(y2(t), t) = diff(K, Y2)}

and using the DEplot tool as