CODE BÀI TẬP LỚN MATLAB PHƯƠNG PHÁP TÍNH THẦY LỘC

THẦY LỘC CODE MATLAB PHƯƠNG PHÁP TÍNH, ĐÂY LÀ MỘT TRONG NHỮNG TÀI LIỆU RẤT QUAN TRỌNG VÀ HỮU ÍCH CHO CÁC BẠN SINH VIÊN ĐANG HỌC TẠI CÁC TRƯỜNG KỸ THUẬT BÁCH KHOA, KHOA HỌC TỰ NHIÊN. CÁC BẠN CÓ THỂ TỰ TÌM HIỂU SÁNG TẠO VÀ GIẢI CÁC BÀI TẬP NHƯNG ĐÂY CHẮC CHẮN SẼ LÀ MỘT TÀI LIỆU THAM KHẢO QUÝ GIÁ ĐỂ TIẾT KIỆM THỜI GIAN CHO TẤT CẢ CÁC BẠN SINH VIÊN ĐANG CẦN.

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);

A(i,i-1) = subs(px,x(i))/(h^2) - subs(qx,x(i))/(2*h);

A(i,i+1) = subs(px,x(i))/(h^2) + subs(qx,x(i))/(2*h);

B(i,1) = subs(fx,x(i));

end

C = inv(A)*B;

Xap = 'Xap xi gia tri cua ham y(0.1) = ';disp(Xap);disp(C(1)); Xap = 'Xap xi gia tri cua ham y(0.5) = ';disp(Xap);disp(C(5)); Xap = 'Xap xi gia tri cua ham y(0.9) = ';disp(Xap);disp(C(9)); end