#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) {
   depth++;
   init_s((S->SI));
 }
}

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;
   }
   printf("All Diagonals =29 except\n");
   for(i=0;i<n;i++) {
    if(A[i*nmax+i]!=29.e0) {
       printf("(%d,%d),%f\n",i,i,A[i*nmax+i]);
    }
   }
   //free(A);
}