Отчет к упражнению 1

x=1.6:0.01:2;

plot(x,-A.^(1./x).*log(A)./x.^2);

x=1.6;e=0.00001;

a=0.3544;r=0.1579;

n=(log(e*(1-a))-log(r))/log(a)

n =

9.7411

A=1+14*pi/30;x=1.6;

for i=2:1:10;

x=A^(1/x);

end

x

x =

1.7007

Отчет к упражнению 2

Записать систему уравнений в виде . Проверить выполнение условия сжатости матрицы . Решить систему уравнений с точностью 0,001.

Проверить решение подстановкой.

syms x1 x2 x3;

X=[x1;x2;x3];

A=[1 0.06 -0.02;0.03 1 -0.05;0.01 -0.02 1];

B=[2;3;5];

C=-A+eye(3)

C =

0 -0.0600 0.0200

-0.0300 0 0.0500

-0.0100 0.0200 0

a=[];

for k=1:3

a=[a;abs(sum(C(k,1:3)))];

end

alpha=max(a)

alpha =

0.0400

X0=B;

X1=B+C*X0

X1 =

1.9200

3.1900

5.0400

r=max(X1-X0)

r =

0.1900

e=0.001;

n=(log(e*(1-alpha))-log(r))/log(alpha)

n =

1.6428

X2=B+C*X1

X2 =

1.9094

3.1944

5.0446

A*X2

ans =

2.0002

2.9995

4.9998