I am solving a set of differential equations using the following command, numericalsol:=dsolve({ODEsys, ics}, numeric, method=rkf45) and I'm plotting my numerical solution using the following command, plots[odeplot](numericalsol,[t,w[1](t)],0..20); But when the solution is plotted the graph looks jagged and not smooth at all. I'm assuming that I am not getting enough points in my numerical solution and that's whats causing the problem. Is there any way to get more points so that my plot will look smooth?

Hi everyone.

In 1-D input and in the Maple calculator x^y^z is treated as ambiguous requiring brackets while in 2-D input it is mysteriously (to me anyway) evaluated as x^(y^z).

What part of my elementary education am I missing here which instructs me to treat  x^y^z as (x^y)^z or so I had always thought?

Seriously I'd appreciate an expert explanation.

I have a set of datapoints, how do I match them to a curve and obtain intercepts and extrema? Ryan

How would I plot the direction field of y'=sin(y),   (pi,2) 

Im really new to maplesoft.

I did the following:

 

restart;
with(Statistics):
CRRA := t->((t)^(1-rho))/(1-rho);
X1 := RandomVariable(Geometric(q));
assume(t::integer, 0 < rho, rho <> 1, 0 < q, q < 1);
ExpectedValue(CRRA(X1));

As result I get a formula, where Maple defines _t0 (as far as I can see in italic font).

But why does Maple do so, and what does this mean?

 

Thx for answering.
 

how do i find the derivative of this? I dont know what to type into maple

C:=x->a(2*x^2+(V/x))+b(8*x+(V/x^2))

 

you have been hired to build a rectangular box with a square base. the volume of this box is to be 1000 cm3. you want to make this box at minimal cost. each of six faces costs [a] cents per square inch and gluing each of the 12 edges costs [b] cents per (linear) inch of length. find an equation that represents the total cost [c] as a function of x, the length of one side of the base. Suppose a=b=1, find the dimensions of the box of minimal cost. How much would this box cost to construct? Use the graph of c(x) to show that you have a minimum.

Problem Statement: In a catalytic reaction that follows a mechanism known as the Michaelis-Menten kinetics, the reactant (S) combines reversibly with a catalyst (E) to form a complex (ES) with a forward and reverse rate constant of k₁ and k_i, respectively. The complex then dissociates into product (P) with a reaction rate constant of k₂ and the catalyst is regenerated.

I would like to accesse the value of a slider in a Maplet, and use it in a computation in a procedure.

My problem is that using Get(Slider1) does not give me a real number to work with. What kind of object is the value of Slider1, and how do I convert it into a real number?

i am trying to solve a set of 12 nonlinear equations using maple but the fsolve command if returning the command text as output.

can any one please help?

 

Hi,

I have a question about integrating an analytic solution to a differential equation in Maple. I have solved a system of 3 odes:
sol := dsolve({initial, syst}, func, numeric)
where
syst := diff(OmegaLambda(x), x) = ode1, diff(OmegaK(x), x) = ode2, diff(lnH(x), x) = ode3.

This then gives me analytic solutions as a function of x to the 3 odes and I can plot, say the third solution, by defining
lnH2 := x-> rhs(op(4, sol(x)))
and plotting lnH2.

Hi all,

This is probably a rookie mistake, but my Maple code is not outputting as I'd expect. E is a list of eigenvectors, s is a list of sets containing coordinates and n is the number of vectors in the list - my code should loop through the elements, and determine the type of stationary point. However, the print function seems to be causing unneccessary trouble.

I'm trying to write a code to find out how many strings of 3 or more same numbers appear in the first 10000 digits of pi. This is what I wrote but I know it's very very inneficient, although it does seem to work I think...

A := NULL:
a := NULL:
for k from 0 to 9997 do
a := [floor(evalf[10000]((10^k)*Pi)) mod 10,
floor(evalf[10000]((10^(k+1))*Pi)) mod 10,
floor(evalf[10000]((10^(k+2))*Pi)) mod 10]:
if a[1]=a[2] and a[1]=a[3] then
A := A+1:
k := k+3
end if:
end do:
A;

 

I am looking for help with basic algorithms to approximate numerical solution to the N body problem. I have been trying to develop basic algorithms using Runge Kutta, Euler ect. Initially I specified to the 3 body case and the equations of motion. I am now trying to develop a more general algorithm. If anyone has knowledge in this are or can look through my annotated procedure to see where I have gone wrong I would appreciate it. I am very new to maple and would appreciate detailed exposition, "for dummy's", if you will. Thanks in advance

Maple should do this

1e123456789876 + 1e123456789876

                           float(infinity)

The answer is 2e123456789876.

Is there a way to get maple to actually do the calculation without doing any manipulation of the numbers before hand?  Maybe there is a package out there for working with extremely large numbers?