Please solve on MATLAB ( computer software)

A steel shaft is to be force fitted into a fixed-ended cast-iron hub The shaft is subjected to a bending moment M, a torque T, and a vertical force P as shown. Suppose that at a point Q in the hubthe stress field is represented by T,--4.7, Tz,-6.45, = 11.8 (MPa) Try If we decide to limit the maxcimum shear stress in the element to 25 MPa, determine the allowed range of o, Note: The maximum shear stress or maxomum principal shear stress is equal to one-half the difference between the largest and smallest principal stresses Hint: White a user defined unchon, whch takes , asut sets up stress matrix calculates ma num shear stress r a, and retums y a. 25 as output then v,such that y = 0 Answer:-35.645 SoS 32.083 MPa o the value of 4.6 76.45 11.8