bài tập môn phương pháp tính cho sinh viên khoa công nghệ thông tin

Bài tập môn phương pháp tính cho sinh viên khoa công nghệ thông tin 1.bảng nội suy ayken và giá trị #include<stdio.h> #include<conio.h> #define max 100 void nhap (float a[max],int n) { int i; for (i=1;i<=n;i++) { printf("a[%d]=",i); scanf("%f",&a[i]); } } void xuat (float a[max],int n) { int i; for (i=1;i<=n;i++) printf("%f",a[i]); } main() { int i,j,n; float x[max],y[max],a[max][max],c; char tl; while (1) { printf("nhap so nguyen duong n="); scanf("%d",&n); printf("\n nhap he so x[i] \n"); nhap (x,n); xuat (x,n); printf("\n nhap he so y[i] \n"); nhap (y,n); xuat (y,n); printf("\n nhap he so c can tinh c="); scanf("%f",&c); for (i=1;i<=n;i++) { a[i][0]=x[i]; a[i][1]=y[i]; printf("%f",a[i][0]); printf("%f",a[i][1]); for (j=2;j<=i;j++) { a[i][j]=(a[j-1][j-1]*(x[i]-c)-a[i][j-1]*(x[j-1]-c))/(x[i]-x[j-1]); printf("%f",a[i][j]); } printf("\n"); } printf("\n vay gia tri can tinh la %.3f",a[n][n]); printf("\n ban co muon tiep tuc khong (c/k)?\n"); tl=getch(); if (tl=='k'||tl=='K') break; } } 2.phương pháp bình phương bé nhất #include<stdio.h> #include<conio.h> #include<math.h> void nhap (float a[10],int n) { int i; for (i=1;i<=n;i++) { printf("a[%d]=",i); scanf("%f",&a[i]); } } void xuat (float a[10],int n) { int i; for (i=1;i<=n;i++) printf("%.3f",a[i]); } void hoanvi (float *a,float *b) { float tmp; tmp=*a; *a=*b; *b=tmp; } void doidong (float a[10][10],int n,int p,int q) { int k; if (p<=n&&q<=n&&p!=q) for (k=1;k<=n+1;k++) hoanvi (&a[p][k],&a[q][k]); } float det (float a[10][10]) { int i,j,k; float m,s=1; for (i=1;i<=3;i++) { if (a[i][i]==0) { for (k=i+1;k<=3;k++) if (a[k][i]!=0) doidong (a,3,i,k); } for (j=i+1;j<=3;j++) { m=-a[j][i]/a[i][i]; for (k=i;k<=3;k++) a[j][k]+=a[i][k]*m; } s*=a[i][i]; } return s; } void giaihephuongtrinh (float b[10][10],float c[10][10],float d[10][10],float e[10][10]) { float p,q,r,l; p=det (b); q=det (c); r=det (d); l=det (e); if (p==0) printf("\n he phuong trinh vo nghiem \n"); else { printf("\n nghiem thu nhat cua he la %.3f",q/p); printf("\n nghiem thu hai cua he la %.3f",r/p); printf("\n nghiem thu ba cua he la %.3f",l/p); printf("\nvay ket qua can tim la:y=%.3f+%.3f*x+%.3f*x*x",q/p,r/p,l/p); } } main() { int n,i; printf("nhap n="); scanf("%d",&n); float x[10],y[10],b[10][10],c[10][10],d[10][10],e[10][10]; float x1=0,xx=0,y1=0,yy=0,xxx=0,xxxx=0,xxy=0,xy=0; nhap (x,n); xuat (x,n); printf("\n"); nhap (y,n); xuat (y,n); for (i=1;i<=n;i++) { x1+=x[i]; y1+=y[i]; xx+=x[i]*x[i]; yy+=y[i]*y[i]; xy+=x[i]*y[i]; xxx+=x[i]*x[i]*x[i]; xxxx+=x[i]*x[i]*x[i]*x[i]; xxy+=x[i]*x[i]*y[i]; } b[1][1]=n; b[1][2]=x1; b[1][3]=xx; b[2][1]=x1; b[2][2]=xx; b[2][3]=xxx; b[3][1]=xx; b[3][2]=xxx; b[3][3]=xxxx; c[1][1]=y1; c[1][2]=x1; c[1][3]=xx; c[2][1]=xy; c[2][2]=xx; c[2][3]=xxx; c[3][1]=xxy; c[3][2]=xxx; c[3][3]=xxxx; d[1][1]=n; d[1][2]=y1; d[1][3]=xx; d[2][1]=x1; d[2][2]=xy; d[2][3]=xxx; d[3][1]=xx; d[3][2]=xxy; d[3][3]=xxxx; e[1][1]=n; e[1][2]=x1; e[1][3]=y1; e[2][1]=x1; e[2][2]=xx; e[2][3]=xy; e[3][1]=xx; e[3][2]=xxx; e[3][3]=xxy; giaihephuongtrinh (b,c,d,e); getch(); } 3.đệ qui quay lui #include<stdio.h> #include<conio.h> #include<math.h> int count; int x[3]; void out (int x[]) { int i; printf("\n%d\t",++count); for ( i = 1; i <= 3; i++ ) printf("%d", x[i] ); } void thu(int i, int x[]) { int j; if (i==1) { for (j=1; j<=5; j++) { x[i]=j; thu(i+1,x); } } if (i==2) { for (j=1;j<=5;j++) { x[i]=j; thu(i+1,x); } } if (i==3) { for (j=1;j<=5;j++) { x[i]=j; out (x); } } } main() { thu (1,x); getch(); } 4.tính định thức ma trận #include<stdio.h> #include<conio.h> #include<math.h> [...]... can tren b="); scanf("%f",&b); printf("\n nhap can tren a="); scanf("%f",&a); printf("\n cac he so yi la:\n"); nhap (y,n); xuat (y,n); printf("\ngia tri cua tich phan la %.3f",tinh (y,n,a,b)); getch (); } 13 .phương pháp bình phương bé nhất #include #include void nhap (float a[10],int n) { int i; for (i=1;i . Bài tập môn phương pháp tính cho sinh viên khoa công nghệ thông tin 1.bảng nội suy ayken và giá trị #include<stdio.h>. can tinh la %.3f",a[n][n]); printf(" ban co muon tiep tuc khong (c/k)? "); tl=getch(); if (tl=='k'||tl=='K') break; } } 2 .phương pháp bình phương. for(j=1;j<=n;j++) printf(" %d",p[i][j]); } getch(); } 6.giải phương trình bằng phương pháp gauss #include<stdio.h> #include<conio.h> #include<math.h>