% This function is used by pwmbp1.m 

function y = f(x)

global m n b w ma

y=m(2*n-1)*x + b(2*n-1) - ma*sin(w*x);