Adri van der Meer

Adri vanderMeer

1355 Reputation

18 Badges

15 years, 144 days
University of Twente (retired)
Enschede, Netherlands

My "website" consists of a Maple Manual in Dutch

MaplePrimes Activity


These are answers submitted by Adri van der Meer

Assuming that u is known, and defined as

u :=  x ->...

you can do:

s := dsolve( {DE,BC }, x(t), numeric, output = listprocedure );
U := u @ subs(s,x(t));
plot( U, 0..1 );

There is an unassigned parameter l in the ODE.

If you change the command in:

sol := dsolve({bc, eqn}, numeric, parameters=[l]);

you get the errormessage:

Error, (in dsolve/numeric) cannot numerically solve a parametric boundary value problem

Additionally: there are 5 ic/bc's and that results in an inconsistent system of equations if you try to solve for the constants in the general solution.
Numerical solutions can be found for 4 initial conditions.

Of course this is possible, but no very easy.

Use the option output=list:

with(LinearAlgebra):
A := <a11,-a11,a13; 0,0,0; 0,0,0>;
eig := Eigenvectors(A, output=list);

s := sort( eig, (x,y)->type(x[1],symbol) ); # if you want the nonzero element(s) first
s := sort( eig, (x,y)->x[1]=0 ); # if yoy want the zero(s)  first
V := Matrix( map( op, s[1..2,3] ) ); # the matrix of eigenvectors

First 25 26 27 Page 27 of 27