进化版雅可比(有点小错误)

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了进化版雅可比(有点小错误)相关的知识,希望对你有一定的参考价值。

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include<math.h>
using namespace std;

int main() {
    freopen("in.txt","r",stdin);
    int a;
    int i,j,s,k=0;
    cin>>a;
    double x[a]= {0},y[a];
    double fc[a][a];
    for(i=0; i<a; i++)
        for(j=0; j<a; j++) {
            cin>>fc[i][j];
        }
    double g[a];
    for(i=0; i<a; i++)
        cin>>g[i];
    double eps=0.0000001;
    double temp;
    int flag[a];
    for(i=0; i<a; i++)
        flag[i]=1;
    do {
        for(i=0; i<a; i++) {
            s=0;
            flag[i]=0;
            for(j=0; j<a; j++) {
                s+=fc[i][j]*x[j]*flag[j];
            }
            flag[i]=1;
            y[i]=(g[i]-s)/fc[i][i];
        }
        for(i=0; i<a; i++) {
            if(temp<fabs(y[i]-x[i])) {
                temp=fabs(y[i]-x[i]);
            }
        }
        for(i=0; i<a; i++)
            x[i]=y[i];
        k++;
        if(k==100) {
            printf("ERROR!\n");
            break;
        }
    } while(temp>eps);
    if(k!=100) {
        for(i=0; i<a; i++) {
            cout<<x[i]<<endl;
        }
    }

    return 0;
}

 

以上是关于进化版雅可比(有点小错误)的主要内容,如果未能解决你的问题,请参考以下文章

CSS网页布局中易犯的10个小错误

一个小错误记录

php 编译小错误

小错误提示

HttpRunnerManager--部署小错误

Oracle中的字符串缓冲区太小错误