Giáo trình tính toán khoa học - Chương 4 pot

... -3 .41 42 -1 .0990 -0 .5858 stopcrit = 1.7764e-016 ========================== V = -0 . 546 8 -0 .6533 0 .44 83 -0 .2706 -0 .44 83 -0 .2706 -0 . 546 8 0.6533 -0 .44 83 0.2706 -0 . 546 8 -0 .6533 -0 . 546 8 ... V = 0.2706 0 .44 83 0.6533 0. 546 8 -0 .6533 -0 . 546 8 0.2706 0 .44 83 0.6533 -0 . 546 8 -0 .2706 0 .44 83...

Ngày tải lên: 11/07/2014, 09:20

... a s {0,1,2,3, ,9}. Chẳng hạn nếu: a = -3 , 141 6 = -( 3.10 0 + 1.10 -1 + 4. 10 -2 + 1.10 -3 + 6.10 -4 ) thì a 0 =3, a -1 = 1, a -2 = 4, a -3 = 1, a -4 = 6. Nếu  a  0,5.10 s thì a s ... 8 Chương 1 MỞ ĐẦU VỀ TÍNH TOÁN KHOA HỌC 1.1 MỞ ĐẦU 1.1.1 Khái niệm về Tính toán khoa học Tính toán khoa học (TTKH) đã trở thàn...

Ngày tải lên: 11/07/2014, 09:20

... -0 .000871189999 14 -0 .000871059200 54 0.0206709 746 2091 0.02066820221695 -0 .2066 042 6208 645 -0 .20657915670952 1.07206 045 679731 1.071 941 0966 044 6 -3 .178 345 4257 941 8 -3 .17801820352051 5.5831276512 646 3 ... inaccurate. RCOND = 2.055969e-018. ans = 1.0e+016 * -0 .45 04 0.9007 -0 .45 04 0.9007 -1 .80 14 0.9007 -0 .45 04 0.90...

Ngày tải lên: 11/07/2014, 09:20

... on; Thí dụ 10. >> peaks z = 3*(1-x).^2.*exp (-( x.^2) - (y+1).^2) - 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) - 1/3*exp (-( x+1).^2 - y.^2) 244 x2 = (R+r*cos(phi))*cos(2*j/n*pi); ... tạo thành dãy khuôn hình trình diễn trong MOVIE. Thí dụ 14. for j=1:n <plot_command> M(j) = getframe; end movie(M);  Các chương trình minh...

Ngày tải lên: 11/07/2014, 09:20

... array 47 (xem phần sai số qui tròn trong chương 1). % Chuong trinh kiem tra anh huong cua sai so qui tron can2 = [ 1 .4 1 .41 4 1 .41 421 1 .41 4213563 ]; format long; for t= can2 a = (t-1)^10; ... inline('x^P1+P3',3) F4 = Inline function: F4(x,P1,P2,P3) = x^P1+P3 2.5 MỘT SỐ LỆNH BỔ TRỢ 2.5.1 Các lệnh tính toán thời gian máy tính Bảng 2-1 0 Danh sá...

Ngày tải lên: 11/07/2014, 09:20

... 3,9 242 2 ,45 63 - 0,22 24 -1 ,3215 0,5506 Giải. Đầu tiên ta lập bảng tính toán: x i y i 2 i x 3 i x 4 i x i i x y 2 i i x y 1 1,5 2 2,5 3 3,5 1,2 341 3,9 242 2 ,45 63 - 0,22 24 - 1,3215 ... 1,2 341 8,82 94 9,8252 -1 ,3900 -1 1,8935 6, 744 8 13,5 6,6213 34, 75 97,875 292,1875 9 ,43 96 13,3501 - Giải hệ phương trình: 292,...

Ngày tải lên: 11/07/2014, 09:20

... 20.0000 -5 .0000 -4 5.0000 30.0000 120.0000 60.0000 -1 5.0000 -1 35.0000 90.0000 360.0000 -4 0.0000 10.0000 90.0000 -6 0.0000 -2 40 .0000 Fy = 8 16 4 -2 0 28 -6 -1 2 -3 15 -2 1 -2 0 -4 0 -1 0 50 -7 0 ... 2 4 1 -5 7 6 12 3 -1 5 21 -4 -8 -2 10 -1 4] >> [px py] =gradient(A) px = 2.0000 -0 .5000 -4...

Ngày tải lên: 11/07/2014, 09:20

... -2 M 2 3 5 -3 2 -5 -5 2 -9 1 0 2 0 1 0 -1 1 -1 2 5/2  k M -5 0 4 0 -5 -8 -1 4 3 0 -2 A 3 A 4 A 8 -4 -2 M 2 3 1 -3 2 1 -5 2 1 1 0 0 0 1 0 -1 ... 2 2 0 1 1 3 -5 -9 2 0 -1 5 -5 0 4 0 -5 -6 -8 -1 4 3 0 -2 -1 8 >> H=Xoay(H,1,3) H = -3 -5 1 0...

Ngày tải lên: 11/07/2014, 09:20

... sqrt((xx-E)^2 +yy^2); yp(1) = y(2); yp(2) = 2*y (4) +y(1 )- E*(y(1)+M)/r1^3-M*(y(1)-E)/r2^3; yp(3) = y (4) ; yp (4) = -2 *y(2) +y(3 )- E*y(3)/r1^3-M*y(3)/r2^3; yp=yp.'; - Sau đây là chương trình ... options=odeset('reltol',1e-5,'abstol',[1e-5 1e-5 1e-5]); 226 BÀI TẬP A. Cài đặt chương trình và lập hàm. 1. Cài đặt chương trình tìm...

Ngày tải lên: 11/07/2014, 09:20

... entries. Sau đây là chương trình tính toán: h m , g m mm gh * h(t), g(t) ( ) ( ) h t g t dt      246 Hãy chú ý quan hệ c k =c -k * . Dấu * dùng để chỉ phép tính lấy liên hợp của ... phân trên rất khó tính, nhưng biểu diễn phổ của tín hiệu thời gian của chúng rất có lợi cho việc tính toán. Ta có: ( ) ikt k k h t h e   , ( ) ilt l l g t g e   . 245...

Ngày tải lên: 11/07/2014, 09:20