sigma_zz expression not right

[rpt26]

I think your expression for sigma_zz is off. If I've parsed the code right then it currently calculates:

but this increases with increasing y forever, which can't be right infinitely far from the dislocation. I think this spills over into the pressure then as well.
I think it should calculate something like:

which might be:
`def sigmaZZ(x_axis,y_axis,Ox,Oy,D,sign):

x,y=np.meshgrid(x_axis,y_axis)

sZZ=-2*nu*D*(((y-Oy)*((x-Ox)**2))/((x-Ox)**2+(y-Oy)**2)**2)

return sZZ`

[rpt26]

BTW, I should say I really like the visualisations, I'll be adding them into my lecture course.

[JQFonseca]

Thank you very much for this. Yes, these were wrong. I am not sure what I was doing there. Thank you so much for taking the time to leave this note, and sorry I only just noticed it!