I am solving for sigma_y and I used the solve function and symbolic math. In the end I need fprintf to display the formula in the form 'v*sigz', but it doesn't do this with the current fprintf. The pictures are there to better explain what I am doing.
syms epsy E sigy v sigx sigz
epsy = 0
sigx = 0
eqn = epsy == (1/E)*(sigy-(v*(sigx+sigz)))
sigy = solve(eqn,sigy)