function f_np1 = f_1st(B,A,M)
f_np1 = conjgrad_trig(A,B);
f_np1 = f_np1 .* M;




