Đang tải... (xem toàn văn)
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 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