开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > matlab例程 > 查看源码
G-S Methom.cpp
资源名称：数值分析课程设计.zip　[点击查看]
上传用户：zhdd911129
上传日期：2007-05-11
资源大小：722k
文件大小：1k
源码类别：

matlab例程

开发平台：

Matlab

G-S Methom.cpp：源码内容
  1. int ReachTOL(double x1[],double x2[],double tol,int n)
  2. {
  3. int i;
  4. double error=0;
  5. for (i=0;i<n;i++)
  6. {
  7. error+=(x1[i]-x2[i])*(x1[i]-x2[i]);
  8. }
  9. if(error<tol) 
  10. return 1;
  11. else
  12. return 0;
  13. }
  15. void GSI(double a[][N],double b[],double x0[],double tol,int max,int n)
  16. {
  17. double x1[N];
  18. int k,i,j;
  19. for (i=0;i<n;i++)                     //Initialazation
  20. {
  21. x1[i]=0;
  22. }
  23. for (k=0;k<max;k++)
  24. {
  25. for (i=0;i<n;i++)
  26. {
  27. for (j=0;j<n;j++)
  28. {
  29.     if(j!=i)
  30. {
  31.     x1[i]+=a[i][j]*x0[j];
  32. } 
  33. }
  34.     x1[i]=((-1)*x1[i]+b[i])/a[i][i];
  35. if(i<n-1)
  36. x0[i]=x1[i];
  37. }
  38. if(ReachTOL(x0,x1,tol,n))
  39. {
  40. for (i=0;i<n;i++)
  41. { 
  42.         x0[i]=x1[i];
  43. } 
  44. return;
  45. }
  46. x00[n-1]=x1[n-1];
  47. for (i=0;i<n;i++)
  48. {
  49.     x1[i]=0;
  50. }
  51. }
  52. for (i=0;i<n;i++)
  53. { 
  54.          x0[i]=x1[i];
  55. }
  56. return;
  57. }