64 } printf("\n\n Nghiem phtrinh: %.3f",c); } else if (f(a)*f(b)>0) printf(" ( %f, %f) khong phai la khoang nghiem",a,b); else if (f(a)==0) printf(" \n Nghiem phtrinh: %.3f",a); else printf(" \n Nghiem phtrinh: %.3f",b); printf("\n\n Ban tiep tuc ko(c/k)?"); tt=getch(); if (tt!='c') break;} } void nhap(float *a, int n) { int i; printf("\n Nhap he so cua phuong trinh:\n"); for (i=0;i<=n; ++i) { printf(" a[%d]=",i); scanf("%f",a+i); } } /* ham tinh gia tri da thuc*/ float f(float x) { float p=d[0]; int i; for(i=1;i<=n;i++) p=p*x+d[i]; return p; } 3. PP tiếp tuyến # include "conio.h" # include "math.h" # define eps 1e-3 float f(float x); float fdh(float x); main() { float a,b; char tt; while (1) { printf("\nNhap xap xi ban dau: "); scanf("%f",&a); /*b=a-f(a)/fdh(a); 65 printf("\n%.3f %.3f %f",a,-f(a)/fdh(a),b);*/ do { b=a; a=b-f(b)/fdh(b); printf("\n%.3f %.3f %f",b,-f(b)/fdh(b),a); } while (fabs(a-b) >= 1e-3 ); printf("\nNghiem phtrinh: %.3f",a); printf("\nTiep tuc ko(c/k)?"); tt=getch(); if (tt=='k' || tt=='K') break;} } float f(float x) { return exp(x)-10*x+7; } float fdh(float x) { return exp(x)-10; } 4. Giải hệ phtrình đại số tuyến tính bằng PP Gauss # include <stdio.h> # include "conio.h" # include "math.h" void nhap(float *a, int n,int m); void xuatmt(float *a, int n,int m); main() { float a[10][10]; float x[10],m,s; char tt; int n,i,j,k; while (1) { printf("\n Nhap n= "); scanf("%d",&n); printf("\n Nhap he so cua he phuong trinh:\n"); for (i=1;i<=n; ++i) for (j=1;j<=n+1;++j) { printf(" pt[%d%d]=",i,j); scanf("%f",&m); a[i][j]=m; } 66 for (i=1;i<=n; i++) { printf("\n"); for (j=1;j<=n+1;j++) printf("%.3f ",a[i][j]); } /* bien doi A ve ma tran tam giac tren */ for(i=1;i<n;i++) for(j=i+1;j<=n;j++) { m=-a[j][i]/a[i][i]; for(k=i;k<=n+1;k++) a[j][k]+=a[i][k]*m; } printf("\n"); for (i=1;i<=n; i++) { printf("\n"); for (j=1;j<=n+1;j++) printf("%.3f ",a[i][j]); } /* tim nghiem theo qtrinh nguoc */ for(i=n;i>=1;i ) { s=a[i][n+1]; for(k=i+1;k<=n;k++) s-=a[i][k]*x[k]; x[i]=s/a[i][i]; } printf("\nNghiem he phtrinh:"); for(i=1;i<=n;i++) printf("%.3f ",x[i]); printf("\n\n Ban tiep tuc ko(c/k)?"); tt=getch(); if (tt!='c') break;} } /* Ham nhap mang a(m,n)*/ void nhap(float *a, int n,int m) { int i,j; printf("\n Nhap he so cua he phuong trinh:\n"); for (i=1;i<=m; i++) for (j=1;j<=n;j++) { printf(" pt[%d%d]=",i,j); scanf("%f",a+i*n+j); } } /* Ham xuat mang a(m,n)*/ void xuatmt(float *a, int n,int m) { int i,j; 67 for (i=1;i<=m; i++) { printf("\n"); for (j=1;j<=n;j++) printf("%.3f ",*(a+i*n+j)); } } 68 TÀI LI ỆU THAM KHẢO [1] Đặng Quốc Lương, Phương pháp tính trong kỹ thuật, Nhà xuất bản xây dựng Hà nội, 2001 [2] Phan Văn Hạp, Giáo trình Cơ sở phương pháp tính tập I,II. Trường ĐH Tổng hợp Hà nội, 1990 [3] Cao quyết Thắng, Phương pháp tính và Lập trình Turbo Pascal. Nhà XB giáo dục, 1998 [4] Tạ Văn Đĩnh, Phương pháp tính. Nhà XB giáo dục, 1994 [5] Dương Thủy Vỹ, Phương pháp tính. Nhà XB khoa học & kỹ thuật, 2001 [6] Phan Văn Hạp, Bài tập phương pháp tính và lập chương trình cho máy tính điện tử . Nhà XB đại học và trung học chuyên nghiệp, 1978 [7] Ralston A, A first course in numberical analysis. McGraw – Hill, NewYork, 1965 . Phương pháp tính. Nhà XB giáo dục, 1994 [5] Dương Thủy Vỹ, Phương pháp tính. Nhà XB khoa học & kỹ thuật, 2001 [6] Phan Văn Hạp, Bài tập phương pháp tính và lập chương trình cho máy tính điện. nhap(float *a, int n,int m); void xuatmt(float *a, int n,int m); main() { float a [10] [10] ; float x [10] ,m,s; char tt; int n,i,j,k; while (1) { printf(" Nhap n= "); scanf("%d",&n);. pháp tính trong kỹ thuật, Nhà xuất bản xây dựng Hà nội, 2001 [2] Phan Văn Hạp, Giáo trình Cơ sở phương pháp tính tập I,II. Trường ĐH Tổng hợp Hà nội, 1990 [3] Cao quyết Thắng, Phương pháp tính