double *****init(Mesh *cm, slv *sv)
{
    int i, j, k, s;
    int sb;
    double *****K;
    memget(K, double ****, cm->nx);
    for(s=0; s<cm->nx; s++) {
        memget(K[s], double ***, sv->sect[s].nm);
        for (i=0; i<sv->sect[s].nm; i++) {
            sb = sv->sect[s].m[i].ne;
            memget(K[s][i], double **, sb);
            for (j=0; j<sb; j++) {
                memget(K[s][i][j], double *, tkb);
                for (k=0; k<tkb; k++) {
                    memget(K[s][i][j][k], double, tkb);
                }
            }
        }
    }
    return K;
}

void zero(slv *sv, Mesh *cm, double *****K1, double *****K2)
{
    int s, i, j, k;
    int val;
    int sb;
    sect *ss;
    for(s=0; s<cm->nx; s++) {
        ss=&(sv->sect[s]);
        val=ss->nm;
        for (i=0; i<val; i++) {
            sb=ss->m[i].ne;
            for (j=0; j<sb; j++)
                for (k=0; k<tkb; k++)
                {
                    memset(K1[s][i][j][k], 0, tkb*sizeof(double));
                    memset(K2[s][i][j][k], 0, tkb*sizeof(double));
                }
        }
    }
}