Materi Metode Numerik

Bahan di download disini

#include
#include

int main()
{
    int d,b,k,n;

    float m[4][5],ms,x[4];
    m[1][1]=3;m[1][2]=2;m[1][3]=1;m[1][4]=6;
    m[2][1]=4;m[2][2]=1;m[2][3]=2;m[2][4]=4;
    m[3][1]=1;m[3][2]=1;m[3][3]=1;m[3][4]=2;
    for (d=1;d<=3;d++)
    {
        for (b=d+1;b<=3;b++)
        {
            ms=m[d][d]/m[b][d];
            for (k=1;k<=4;k++)
            {
                    m[b][k]=m[b][k]*ms-m[d][k];
            }
        }
    }
    for (b=1;b<=3;b++)
    {
        for (k=1;k<=4;k++)
        {

            printf("%f ",m[b][k]);
        }
        printf("\n");
    }

for (b=3;b>0;b--)
{
    x[b]=m[b][4];
    for(k=b+1;k<=3;k++)
    {
        x[b]=x[b]-m[b][k]*x[k];
    }
    x[b]=x[b]/m[b][b];

}
for (b=1;b<=3;b++)
{
    printf("x(%d)=%f",b,x[b]);
}
    return 0;

}