I need change this code in the " Matlab"

import numpy as np import matplotlib.pyplot as plt L 59 D 084 R 10 S 10000 r max10 N-100 plt.figure0 flux bS'np.exp(-L) 4 np.pirD) flux a S/4 np pi*D)*np.sinh((R-psnh(R/L) plt.plotir, flux_a, color 'red,labela)Sphere) pltplot(r,flu pltxlabelr) pltylabellux") pltlegendo plt.savefig(HW3_q2.eps) pltshow0