please i want to write the maple programme  for solving y^/ =y, y(0)=1, 0≤x≤1 with h=0.1 using the following three schemes simultaneously:

1. (7/12)y_n+1 - y_n+1/2 +5/12y_n =(h/12)f_n+2

2. y_n+2 - (4/3)y_n+1 +1/3y_n =(2h/3)f_n+2

3. y_n+1 -y_n =hf_n+1/2

How do you enter h in a Euler's method problem? I need to change h=0.1 to h=0.05

I'm using the 2D input in Maple 14.

So far I have this:

> DE2 := diff(y(t), t) = 3+t-y(t);

> Euler(DE2, y(0) = 1, t = .4);

> Euler(DE2, y(0) = 1, t = .4, output = Error);

> Euler(DE2, y(0) = 1, t = .4, output = information);

When I learning "Maple introductory programming guide", there is a little procedure attracting me:

_________________________________

NewName := proc ()

    global N;

    N := N+1;

    while assigned(cat(C, N)) do

        N := N+1;

    end do;

    cat(C, N);

end proc

_________________________________

Dear primers,

I have a curve that shows 4 maxima (green curve in the image, red curve is the derivative). Now i simply want the coordinates of these maxima without defining a threshold. I tried the "ExtremePoints" function, but it returns an error ("...cannot determine discontinuities") and the "extrema" function which gives me only one max and min value.It would be great, if somebody could give me a hint...

Normally after Maple calculations i obtain results like this one:

(-1+du1(t))/(C*L)

when i would very much like to obtain it with a slight difference:

(du1(t)-1)/(C*L)

Is there a way i can achieve this?

Thanks!

Hi all,

When I solved:

restart: mu:=solve(h*sin(mu*b)+mu*cos(mu*b)=0,mu,AllSolutions=true);

where b and h are constant, the solution is:

                                                    ...

Let us consider the simple single-variable function,y=f(x)=x*x(i,e,x^2).Now suppose that we don't know our func is x*x,instead we are given:

                                x                        ...

Hello,

I'm currently working on a project where I need to create a 3D graph with close to 5500 data points. For smaller datasets, I have been able to create the graphs using the with (plots): pointplot3d command. An example of syntax used: 

with (plots): pointplot3d({[23.8, 2399, 24.2],[23.6, 3030, 24.0],[23.6, 3876, 24.2]},axes=normal,symbol=box);

However, when I try to use the same method to graph the larger data set, the program freezes ...

whats the difference between A.B and A&*B given that both A and B are matrices.

Hello,

I would like to create a model of a sinngle pendulum with Maple 13

I began but I cannot manage to animate it

Here is what i have already written:

pendulum:=proc(l,initialAngle,initialSpeed)

local eq,sol,angledata,x,y,plotseq

eq:=diff(θ(t),t,t) + 9.81/l sin(θ(t)) = 0;

sol:=dsolve({eq,θ(0)=initialAngle,D(θ)(0)=initialSpeed},θ(t),numeric);

I am trying to get the equation of a line. I have a buch of points and need the equation of a line going through those points. I am making computer program that can finde the y value of any point on the line with the equation I'm trying to get and an X value.

So, I have some points,..

I can get the graph,..but not the equation for that graph...

Is this posable,...?? I can't seem to find out how...

Thank you!

(I put in the points Im using just as are ference)

Hi all,

I don't find an easy way to do what follows.

I've a list of names, with assumptions on them, eg:

assume(a,real);
assume(b,real):

parameters:=[a,b]

then I try do to something like

i:=0;
for what in parameters do
temp:=i;

assign(cat('crlb_', what), temp);
i:=i+1;
end do

because I want to create variables which keep a trace of the parameters they describe in the last part of their name.

ifplease how can i factorise in powers of (x-x_n), the expression  y(t):=(1/3h²)[(-x+x_n+h)(-x+x_n+3h)]y_n+(1/3h²)[(-x_n+x)(-x+x_n+4h)]y_n+1 -(1/3h)[(-x_n+x)(-x+x_n+3h)]_fn+2?

i was trying to familiarize with maple 13 but encountered this. pls how may i go out from here?

> with(Student[NumericalAnalysis]);
> DE1 := d*y(t)/dt = y(t)-t^2+1;
                           d y(t)           2   
        ...

I need to differentiate a postion vector with respect to time and extract the theta dots in the following http://img543.imageshack.us/img543/3079/positionvector.jpg

Many Thanks