Hosted by uCoz
 

Robert M. Corless
Department of Applied Mathematics
University of Western Ontario
London, Canada

Copyright 2001 by Robert M. Corless
All rights reserved

Useful one-word commands 

Использование монокоманд

Solving equations 

Решение уравнений

Linear equations 

Линейные уравнения

> restart;
> with(LinearAlgebra):
> Yee := Matrix( 8,8,
[[e^(-1), e^(-1), e^(-1), e^(-1), e^(-1), e^(-1), e^(-1), 0],
[1, 1, 1, 1, 1, 1, 0, 1],
[1, 1, 1, 1, 1, 0, 1, 1],
[1, 1, 1, 1, 0, 1, 1, 1],
[1, 1, 1, 0, 1, 1, 1, 1],
[1, 1, 0, 1, 1, 1, 1, 1],
[1, 0, 1, 1, 1, 1, 1, 1],
[0, e, e, e, e, e, e, e]]);

Yee := _rtable[14287132]
> P,L,U,R,det := LUDecomposition(Yee,output=['P','L','U1','R','determinant']):
> Norm(Yee-P.L.U.R ,infinity);

0
> det;

-7
> P;

_rtable[14291808]
> L;

_rtable[12231876]
> U;

_rtable[14296580]
> L^(-1);

_rtable[12376464]
> U^(-1);

_rtable[14303684]
> P^(-1);

_rtable[14308928]
> Norm(U^(-1).L^(-1).P^(-1).Yee-IdentityMatrix(8,8),infinity);

0

С официального разрешения                    © 2002 Waterloo Maple, Inc