SOMA[M_,ts_,tz_,n_]:=(Clear[ij,i,di];di=(tz-ts)/n;For[ij=0,ij<n+1,ij++,n1013[ts+di*ij]=(ts+di*ij)^(2/3)(1+(2Sin[M(ts+di*ij)])/(3M(ts+di*ij)));];Table[{ts+di*i,n1013[ts+di*i]},{i,0,n}])