%% This function is used by pwmbp1.m 
%alpha + delta for pwm
function y = f(x)

global m n b w ma

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