Write the MATLAB code for the following equation (if R,Y axis B=0.1,r=1 M=5,15,40 Three plot in same code Ci= cosine integral function si= sine integral function R=ln(2)2Bi=0Mrr,ibr,i(ci(r,i)cos(r,i)si(r,i)sin(r,i)) If the br,i br,i=r(Mr)(1)i(Mri) and r,i=(i+r)