clangエラー c++ error

Question

clang:error:linker command failed with exit code 1 (use -v to see invocation)
という風に表示されるのですが，調べたところライブラリが足りていない時のエラーだと書いてありました．これは容量不足ということでしょうか？

どのように対処すれば良いですか？
いらないファイルを削除したりゴミ箱を空にしたりしましたがダメでした．

OSはmacで，Visual Studio Codeを使って書いていて，ターミナルでclang++ -std=c++11 newton.cppで実行しています．

一応コードも貼っておきます．

#include<stdio.h>
#include<math.h>
#include<stdlib.h>

#define EPS pow(10.0,-8.0) //epsilonの設定
#define KMAX 100 //最大反復回数
#define N 4 //４元方程式

double **dmatrix(int nr1,int nr2,int nl1,int nl2);
void free_dmatrix(double **a,int nr1,int nr2,int nl1,int nl2);
double *dvector(int i,int j);
void free_dvector(double *a,int i);
double *gauss(double **a,double *b);
double vector_norm1(double *a,int m,int n);

double f(double x,double y,double z,double u); //1行目の式をfとおいた．2行目以降も同様．
double g(double x,double y,double z,double u);
double h(double x,double y,double z,double u);
double i(double x,double y,double z,double u);

double f_x(double x,double y,double z,double u); //fをxで偏微分した式をf_xとおいた．以下同様．
double f_y(double x,double y,double z,double u);
double f_z(double x,double y,double z,double u);
double f_u(double x,double y,double z,double u);

double g_x(double x,double y,double z,double u);
double g_y(double x,double y,double z,double u);
double g_z(double x,double y,double z,double u);
double g_u(double x,double y,double z,double u);

double h_x(double x,double y,double z,double u);
double h_y(double x,double y,double z,double u);
double h_z(double x,double y,double z,double u);
double h_u(double x,double y,double z,double u);

double i_x(double x,double y,double z,double u);
double i_y(double x,double y,double z,double u);
double i_z(double x,double y,double z,double u);
double i_u(double x,double y,double z,double u);

void newton2(double x,double y,double z,double u);

int main(void){
    double x,y,z,u;
    printf("初期値を入力してください-->x0,y0,z0,u0\n");
    scanf("%lf %lf %lf %lf",&x,&y,&z,&u);
    newton2(x,y,z,u);
    return 0;
}

void newton2(double x,double y,double z,double u){
    int j,k=0;
    double *xk,*d,**J;

    d = dvector(1,N);
    xk = dvector(1,N);
    J = dmatrix(1,N,1,N); //N*Nのヤコビ行列

    xk[1]=x; xk[2]=y; xk[3]=z; xk[4]=u;

    do{
        d[1] = -f(xk[1],xk[2],xk[3],xk[4]);
        d[2] = -g(xk[1],xk[2],xk[3],xk[4]);
        d[3] = -h(xk[1],xk[2],xk[3],xk[4]);
        d[4] = -i(xk[1],xk[2],xk[3],xk[4]);

        J[1][1] = f_x(xk[1],xk[2],xk[3],xk[4]);
        J[1][2] = f_y(xk[1],xk[2],xk[3],xk[4]);
        J[1][3] = f_z(xk[1],xk[2],xk[3],xk[4]);
        J[1][4] = f_u(xk[1],xk[2],xk[3],xk[4]);

        J[2][1] = g_x(xk[1],xk[2],xk[3],xk[4]);
        J[2][2] = g_y(xk[1],xk[2],xk[3],xk[4]);
        J[2][3] = g_z(xk[1],xk[2],xk[3],xk[4]);
        J[2][4] = g_u(xk[1],xk[2],xk[3],xk[4]);

        J[3][1] = h_x(xk[1],xk[2],xk[3],xk[4]);
        J[3][2] = h_y(xk[1],xk[2],xk[3],xk[4]);
        J[3][3] = h_z(xk[1],xk[2],xk[3],xk[4]);
        J[3][4] = h_u(xk[1],xk[2],xk[3],xk[4]);

        J[4][1] = i_x(xk[1],xk[2],xk[3],xk[4]);
        J[4][2] = i_y(xk[1],xk[2],xk[3],xk[4]);
        J[4][3] = i_z(xk[1],xk[2],xk[3],xk[4]);
        J[4][4] = i_u(xk[1],xk[2],xk[3],xk[4]);

        d = gauss(J,d);
        for(j=1;j<=N;j++) xk[j] += d[j];
        k++;
    }while(vector_norm1(d,1,N)>EPS && k<KMAX);

    if(k==KMAX){
        printf("答えが見つかりませんでした\n");
    }else{
        printf("答えは x=%f, y=%f, z=%f, u=%f　です\n",xk[1],xk[2],xk[3],xk[4]);
    }

    free_dmatrix(J,1,N,1,N);
    free_dvector(d,1);
    free_dvector(xk,1);
}

double f(double x,double y,double z,double u)
{
    return(x*x*x + 2*y*y + z*z - 2*u*u - z + 2);
}
double g(double x,double y,double z,double u)
{
    return(y*y - x*y - 2*x*z - 5);
}
double h(double x,double y,double z,double u)
{
    return(3*x*y*z + 2*y*u - u + 3);
}
double i(double x,double y,double z,double u)
{
    return(y*y*z - x*x + 2*y*z - x - 2*z - u);
}

double f_x(double x,double y,double z,double u)
{
    return(3*x*x);
}
double f_y(double x,double y,double z,double u)
{
    return(4*y);
}
double f_z(double x,double y,double z,double u)
{
    return(2*z - 1);
}
double f_u(double x,double y,double z,double u)
{
    return(-4*u);
}

double g_x(double x,double y,double z,double u)
{
    return(-y-2*z);
}
double g_y(double x,double y,double z,double u)
{
    return(2*y-x);
}
double g_z(double x,double y,double z,double u)
{
    return(-2*x);
}
double g_u(double x,double y,double z,double u)
{
    return 0;
}

double h_x(double x,double y,double z,double u)
{
    return(3*y*z);
}
double h_y(double x,double y,double z,double u)
{
    return(3*x*z + 2*u);
}
double h_z(double x,double y,double z,double u)
{
    return(3*x*y);
}
double h_u(double x,double y,double z,double u)
{
    return(2*y - 1);
}

double i_x(double x,double y,double z,double u)
{
    return(-2*x-x);
}
double i_y(double x,double y,double z,double u)
{
    return(2*z*y + 2*z);
}
double i_z(double x,double y,double z,double u)
{
    return(y*y + 2*y -2);
}
double i_u(double x,double y,double z,double u)
{
    return -1;
}

Comments

[y_waiwai]

エラーメッセージはいらぬ省略翻訳しないでそのまま提示しましょう

必要な情報隠されては答えようがないです

[hoshi-takanori]

↓ これらの関数が宣言だけされてて定義がないからでは。

double **dmatrix(int nr1,int nr2,int nl1,int nl2);
void free_dmatrix(double **a,int nr1,int nr2,int nl1,int nl2);
double *dvector(int i,int j);
void free_dvector(double *a,int i);
double *gauss(double **a,double *b);
double vector_norm1(double *a,int m,int n);

[episteme]

clang++ -std=c++11 newton.cpp したってんならその↑とおりでしょね。