if 'y' in vari:
                        Ss3 = simps(S * st ** 3, x=theta, axis=0)
                        Ss4 = simps(S * st ** 4, x=theta, axis=0)
                        vec.append(ky ** 4 * Ss4)
                        mtext.append('myyyy')