Gett ing started
Every Maple worksheet should begin by re-initializing the Maple "kernel" and loading the additional packages that we are most likely to use.
> restart;
> with( plots ):
> with( DEtools ):
>
DO NOT CHANGE ANYTHING IN THE FOLLOWING INPUT REGION - YOU HAVE BEEN WARNED !!!
> Euler := proc( MODEL, IC, DOMAIN, N )
> local a, b, dt, dy, f, i, tt, yy, LISTpts, vars;
> vars := (op(2,lhs(MODEL) ), op([1,0],lhs(MODEL) ) );
> f := unapply( rhs(MODEL), vars );
> a := op(1,rhs(DOMAIN)); b := op(2,rhs(DOMAIN));
> dt := (b-a)/N;
> tt := a;
> yy := rhs(IC);
> LISTpts := [ tt, yy ];
> for i from 1 to N do
> dy := evalf( f(tt,yy) ) * dt;
> yy := yy + dy;
> tt := tt + dt;
> LISTpts := LISTpts, [ tt, yy ];
> od; RETURN( [ LISTpts ] );
> end: