Convert this C program into matlab. Please help. Stuck very badly as con[i][j] matrix is showing decimal nos. like 23.33(which is obviously not accepted in indexing in matlab).That error is showing in MG_matrix[con[eltn][i] - 1][con[eltn][j] - 1]

for (ip=0;ip<9;ip++)
{
Ms3_matrix[i][j]=Ms3_matrix[i][j]+w3(epsilon3)*J*NTM_matrix[i][ip]*N_matrix[ip][j];   
Mte3_matrix[i][j]=Ms3_matrix[i][j];                                                 
}
 for (i = 0; i <27; i++)
 {
	    for (j = 0; j < 27; j++)
          {
                Mte6_matrix[i][j]=0;
}
}
 for (i = 0; i <27; i++)
         {
	    for (j = 0; j < 27; j++)
          {
                Mte_matrix[i][j]=0;
/*  *****************************Total Element Mass Matrix or Local Mass matrix(Kte_matrix)*************************** */
		Mte_matrix[i][j]= Mte_matrix[i][j]+Mte3_matrix[i][j]+Mte6_matrix[i][j];
	      }
}
for (i = 0; i <30; i++)
	{
	  for (j = 0; j <27; j++)
        con[i][j]=58*(j/3)+j+2*i+1
    }
}
for (i = 0; i <27; i++)

	{

	  for (j = 0; j <27; j++)

	    MG_matrix[con[eltn][i] - 1][con[eltn][j] - 1] = Mte_matrix[i][j];	(ERROR HERE)

	}

[KSSV]

Virajan Verma this is bad.....how do you expect people to do your work? If you have made any attempts..show us that..we are here to help/ guide you in MATLAB.

[Virajan Verma]

I have made several attempts and failed. I am on this very line for the past 10 days and not getting any solution. That is why i am asking.

[Virajan Verma]

Attempted to access MG_matrix(26.3333,26.3333); index must be a positive integer or logical.

Error in practisee2 (line 574)

MG_matrix(conn(eltn,i) -1,conn(eltn,j)-1)= Mte_matrix(i,j);

This is the error i am getting.

[KSSV]

That case post the MATLAB code.

[Virajan Verma]

for ip=1:9         %%DOUBT
Ms3_matrix(i,j)=Ms3_matrix(i,j)+w3(epsilon3)*J*NTM_matrix(i,ip)*N_matrix(ip,j);   
Mte3_matrix(i,j)=Ms3_matrix(i,j);    
end
for i=1:27
     for j=1:27
       Mte_matrix(i,j)=0;
     end
 end
 
  for i=1:27
     for j=1:27
       Mte6_matrix(i,j)=0;
     end
 end
for ip=1:27
                    Mte_matrix(i,j)= Mte_matrix(i,j)+Mte3_matrix(i,j)+Mte6_matrix(i,j);
        end
        for i=1:30
    for j=1:27
        conn(i,j)=58*(j/3)+j+2*i+1
    end
end
     eltn=30;
for i=1:27
     for j=1:27
  

	    MG_matrix(conn(eltn,i) -1,conn(eltn,j)-1)= Mte_matrix(i,j);

[Virajan Verma]

Please help me get over this. M stuck in this for the past 10 days. I really need this fixed ASAP

Answer 1

[Torsten]

for i=1:30
  for j=1:27
    conn(i,j) = 58*floor((j-1)/3)+(j-1)+2*(i-1)+1
  end
end

for i=1:27
  for j=1:27
     MG_matrix(conn(eltn,i) ,conn(eltn,j))= Mte_matrix(i,j);
  end
end  

Answer 2

[Virajan Verma]

@Torsten. How can i thank you. U have really solved this obstacle i was on since many days. Really appreciate ur skills and thanks for ur time.