ĐỀ PHƯƠNG PHÁP TÍNH THẦY LỘC

CODE BÀI TẬP LỚN PHƯƠNG PHÁP TÍNH THẦY LỘC
CÂU 1:
function Cau1
clc;
format short
syms x
mn = input('Nhap 2 so cuoi MSSV: ');
M = (mn+12)/10;
%%%%%%%%%%%%%%
f = exp(x)+2*x^2+sin(x)/M-10;
a = 1;
b = 2;
if double(subs(f,x,a)*subs(diff(diff(f,x),x),x,a))>0
X = a;
else
X = b;
end
if double(subs(diff(f,x),x,a))<double(subs(diff(f,x),x,b));
min = subs(diff(f,x),x,a);
else
min = subs(diff(f,x),x,b);
end
for i = 1:2 %chay nghiem tu x1 -> x2

clc;
mn = input('Nhap 2 so cuoi MSSV: ');
M = (mn + 12)/10;
%%%%%%%%%%%%%%%%%%%%%%%%
a = [12+M 2 -3 4 5; 4 13+M 4 -2 -6; 3 -3 14+M 2 -5; 2 -2 4 15+M -3; 5 -4 5 -3 16+M];
b = [9; 8; 7; 6; 5];
x = [1.5; 0.3; 3.4; 1.4; 5.6];
i=0;A=x(1);B=x(2);C=x(3);D=x(4);E=x(5);
while i<3 %tim nghiem x3
X=(b(1)-a(1,2)*B-a(1,3)*C-a(1,4)*D-a(1,5)*E)/a(1,1);
Y=(b(2)-a(2,1)*A-a(2,3)*C-a(2,4)*D-a(2,5)*E)/a(2,2);
Z=(b(3)-a(3,1)*A-a(3,2)*B-a(3,4)*D-a(3,5)*E)/a(3,3);
T=(b(4)-a(4,1)*A-a(4,2)*B-a(4,3)*C-a(4,5)*E)/a(4,4);
K=(b(5)-a(5,1)*A-a(5,2)*B-a(5,3)*C-a(5,4)*D)/a(5,5);
A = X; B = Y; C = Z; D = T; E = K;
i=i+1;
end
x1 = A, x2 = B, x3 = C, x4 = D, x5 = E,
end
CÂU 4:
function GaussSeidel
clc;
mn = input('Nhap 2 so cuoi MSSV: ');
M = (mn + 12)/10;
a = [12+M 2 -3 4 5; 4 13+M 4 -2 -6; 3 -3 14+M 2 -5; 2 -2 4 15+M -3; 5 -4 5 -3 17+M];
b = [9; 8; 7; 6; 4];
x = [0.1; 0.3; 0.4; 0.5; 0.9];
i=0;A=x(1);B=x(2);C=x(3);D=x(4);E=x(5);
while i<3 %tim nghiem x3
A=(b(1)-a(1,2)*B-a(1,3)*C-a(1,4)*D-a(1,5)*E)/a(1,1);

end
c = inv(A)*B;
for i=1:n-1
b(i)=(Y(i+1)-Y(i))/h(i)-h(i)*(c(i+1)+2*c(i))/3;
d(i)=(c(i+1)-c(i))/(3*h(i));
end
t = 1.4;
I = 0;
for i =1:n-1
if t>= X(i) && t<X(i+1)
I = Y(i) + b(i)*(t-X(i)) + c(i)*(t-X(i))^2 + d(i)*(t-X(i))^3;
end
end
fprintf('Xap xi gia tri cua ham tai x = %.1f',t);
disp(I);
t = 2.5;
I = 0;
for i =1:n-1
if t>= X(i) && t<X(i+1)
I = Y(i) + b(i)*(t-X(i)) + c(i)*(t-X(i))^2 + d(i)*(t-X(i))^3;
end
end
fprintf('Xap xi gia tri cua ham tai x = %.1f',t);
disp(I);
end
CÂU 6:
function Cau6
clc;
format short ;
syms x m real

I = Y(i) + b(i)*(t-X(i)) + C(i)*(t-X(i))^2 + d(i)*(t-X(i))^3;
end
end
fprintf('Xap xi gia tri cua ham tai x = %.1f',t);
disp(I);
t = 3.0;
I = 0;
for i =1:n-1
if t>= X(i) && t<X(i+1)
I = Y(i) + b(i)*(t-X(i)) + C(i)*(t-X(i))^2 + d(i)*(t-X(i))^3;
end
end
fprintf('Xap xi gia tri cua ham tai x = %.1f',t);
disp(I);
end
CÂU 7:
function Cau7
clc;
format short
syms x m
M = input('Nhap M: ');

%%%%%%%%%%%%%%%%%%%%
u = [1.2 1.3 1.4 1.5 1.7];
y = [2*M 2.5 5 4.5 5.5];
gx = sqrt(x.^2+1);
hx = cos(x);
A = zeros(2); B = zeros(2,1);
for i = 1:5
A(1,1)=A(1,1)+subs(gx.^2,u(i)) ;

end
for i = 1:(size(c,2)-1)
d(i) = (c(i+1)-c(i))/(p(i+4)-p(i));
x4 = x4*(x-p(i));
end
for i = 1:(size(d,2)-1)
e(i) = (d(i+1)-d(i))/(p(i+5)-p(i));
x5 = x5*(x-p(i));
end
A = subs(diff((q(1)+a(1)*x5+b(1)*x4+c(1)*x3+d(1)*x2+e(1)*x1),x),x,0.5);
disp('gia tri xap xi dao ham cap 1 cua ham tai x=0.5 la:');disp(double(A));
end
CÂU 9:
function Cau9
clc;
syms x y;
mn = input('Nhap 2 so cuoi MSSV: ');
M = (mn+12)/10;
%%%%%%%%%%%%%%%%
a= 2;
b = 62;
n =120;
h= (b-a)/n;
m=n/2;
f= (2*M*x^2+x+1)/(7*x^4+x+6);
y= subs(f,x,a)+subs(f,x,b)+4*subs(f,x,a+h);
for i = 1 :m-1
y= y+ 2*subs(f,x,a+2*i*h)+4*subs(f,x,a+(2*i+1)*h);
end
y=double(y*h/3);

disp('Xap xi: '), y(n),
end
CÂU 11:
function Cau11
clc;
format short
syms x real
mn = input('Nhap 2 so cuoi MSSV: ');
M = (mn+12)/10;
a = 0; b = 1;
h = 0.1; ya = 1; yb = 1.2;
n = (b-a)/h;
px = x + 2*M;
qx = x.^3;
rx = -30;
fx = -x.*(x+1);
for i = 1:n
x(i) = a + i*h;
end
n = round(n);
A = zeros(n-1);
B = zeros(n-1,1);
A(1,1) = subs(rx,x(1)) - 2*subs(px,(x(1)))/(h^2);
A(1,2) = subs(px,x(1))/(h^2) + subs(qx,x(1))/(2*h);
A(n-1,n-1) = subs(rx,x(n-1)) - 2*subs(px,x(n-1))/(h^2);
A(n-1,n-2) = subs(px,x(n-1))/(h^2) - subs(qx,x(n-1))/(2*h);
B(1,1) = subs(fx,x(1)) - subs((px/(h^2)-qx/(2*h)),x(1))*ya;
B(n-1,1) = subs(fx,x(n-1))-subs((px/(h^2)+qx/(2*h)),x(n-1))*yb;
for i=2:n-2
A(i,i) = subs(rx,x(i)) - 2*subs(px,(x(i)))/(h^2);