maple矩阵求逆

maple矩阵求逆

with(LinearAlgebra):
a:=RandomMatrix(3,3);
E:=Matrix(3,3,shape=identity);
A:=;
B:=RowOperation(%,[2,1],(-A[2,1])/(A[1,1]));
C:=RowOperation(%,[3,1],(-A[3,1])/(A[1,1]));
d:=RowOperation(%,[1,2],(-C[1,2])/(C[2,2]));
F:=RowOperation(%,[3,2],(-d[3,2])/(d[2,2]));
g:=RowOperation(%,[1,3],(-F[1,3])/(F[3,3]));
G:=RowOperation(%,[2,3],(-F[2,3])/(F[3,3]));
c:=RowOperation(%,1,(1/G[1,1]));
f:=RowOperation(%,2,(1/G[2,2]));
H:=RowOperation(%,3,(1/G[3,3]));
a逆:=SubMatrix(%,1..3,4..6);
a对角:=SubMatrix(G,1..3,1..3);

a:=[-2,-1,1]; b:=[2,2.5,3.5];
f:=interp(a, b, x); with(plots):
p:=pointplot([[-1,2.5],[-2,2],[1,3.5]]):
g:=plot(f,x=-3..3,color=red):
display(p,g);

with(LinearAlgebra);
a:=RandomVector(3);
b:=RandomVector(3);
c:=RandomVector(3);
s:=GramSchmidt([a,b,c],normalized);
a1:=s[1];b1:=s[2];c1:=s[3];
evalb(DotProduct(a1,b1)=0);
evalb(DotProduct(a1,c1)=0);
evalb(DotProduct(b1,c1)=0)；

1
restart;
f:=x->sqrt(1+x);
g:=convert(series(f(x),x=0,10),polynom);
a:=evalf(sqrt(2),15);
y:=subs(x=1,g);
n:=10:
while abs(a)>=1e-7 do
a:=a-y;g:=convert(series(f(x),x=0,n+1),polynom);y:=subs(x=1,g)
end do;
y;
with(share):
readshare(FPS,analysis):
g:=FormalPowerSeries(sqrt(1+x),x=0);
sqrt(2):=evalf(subs(x=1,g),8);#?
sqrt(3):=evalf(subs(x=2,g),8);
2
p:=randpoly(x,degree=6);
q:=randpoly(x,degree=2);
s:=quo(p,q,x,'r');
r;
# check
evalb(p=simplity(s*q+r));#?
3
with(LinearAlgebra):
a:=RandomMatrix(3,3);
E:=Matrix(3,3,shape=identity);
A:=<a|E>;
B:=RowOperation(%,[2,1],(-A[2,1])/(A[1,1]));
C:=RowOperation(%,[3,1],(-A[3,1])/(A[1,1]));
d:=RowOperation(%,[1,2],(-C[1,2])/(C[2,2]));
F:=RowOperation(%,[3,2],(-d[3,2])/(d[2,2]));
g:=RowOperation(%,[1,3],(-F[1,3])/(F[3,3]));
G:=RowOperation(%,[2,3],(-F[2,3])/(F[3,3]));
c:=RowOperation(%,1,(1/G[1,1]));
f:=RowOperation(%,2,(1/G[2,2]));
H:=RowOperation(%,3,(1/G[3,3]));
a逆:=SubMatrix(%,1..3,4..6);
a对角:=SubMatrix(G,1..3,1..3);
4
restart;
with(LinearAlgebra);
a := RandomVector(3);
b := RandomVector(3);
c := RandomVector(3);
s := GramSchmidt([a, b, c], normalized);
a1 := s[1];
b1 := s[2];
c1 := s[3];
evalb(DotProduct(a1, b1) = 0);
evalb(DotProduct(a1, c1) = 0);
evalb(DotProduct(b1, c1) = 0);
5
a:=[-1,-2,1];
b:=[2.5,2,3.5];
f:=interp(a,b,x);#
a:=[-1,-2,1,2];
b:=[2.5,2,3.5,7];
f1:=interp(a,b,x);
with(plots):
g1:=pointplot([[-1,2.5],[-2,2],[1,3.5]],axes=boxed):
g2:=plot(f,x=-3..3):
g3:=plot(f1,x=-3..3):
display(g1,g2,g3);