#include <stdlib.h>
#include <stdio.h>
   struct S0 {
    char name;
    double x;
    float b[3];
    double y;
    struct S0 *SI;
   };
#define DEPTH 4
int depth=0;
void i_matrix(int n,int lda,float *A) {
 int i,j;
 for(i=0;i<n;i++) {
  for(j=0;j<n;j++) {
   if(i==j) {A[i*lda+j]=4.0f;}
   else     {A[i*lda+j]=(i+j)/(float)(n);}
 }
 }
}
void init_s(struct S0 S) {
 int i;
 S.name='x';
 S.x=2.5e1;
 for(i=0;i<3;i++) S.b[i]=2.f+(float)i;
 S.y=-7.e0;
 S.SI=(struct S0*)malloc(sizeof(struct S0));
 if(depth<DEPTH) {
   init_s(*(S.SI));
   depth++;
 }
}

int main(int argc, char **argv) {
   int i,n,nmax=20;
   float *A=(float*)malloc(sizeof(float)*nmax*nmax);
   struct S0 Structure;

   init_s(Structure);
   n=10;
   if(n<nmax) 
      i_matrix(n,n,A);
   for(i=0;i<n;i++) { 
    A[i*nmax+i]+=Structure.x;
   }
   //free(A);
}