thanx once again iyami,
actually the corrections which u told like
dv_y/dt = 1/ρ (〖dσ〗_xy/dx + 〖dσ〗_yy/dy+ 〖dσ〗_yz/dz) to be as
dv_y/dt=1/ρ (〖dσ〗_yx/dx +〖dσ〗_yy/dy+〖dσ〗_yz/dz)
does'nt affect coz both the stress matrix r same xy=yx
as i have defined initially in my code
tau_xx = zeros(yDim,xDim,zDim);
tau_yy = zeros(yDim,xDim,zDim);
tau_zz = zeros(yDim,xDim,zDim);
tau_xy = zeros(yDim-1,xDim-1,zDim);
tau_yz = zeros(yDim-1,xDim,zDim-1);
tau_xz = zeros(yDim,xDim-1,zDim-1);
ux = zeros(yDim,xDim-1,zDim);
uy = zeros(yDim-1,xDim,zDim);
uz = zeros(yDim,xDim,zDim-1);
but i will try once again using your advice.
i m not getting any clue why is the problem arising ....
---------- Post added at 09:32 ---------- Previous post was at 09:13 ----------
i got the bug... something like as u asked me to check ...
thanx man ... i really appreciate your help..
the code is converging now