please make a script to solve the following equation 2cos(2π200t+π/3)+bcos(2π200t+ψ)=cos(2π200t+π2/3) with b>0 and -π<ψ<π find b and ψ in matlab

main.m
syms b psi t
eqn = 2*cos(2*pi*200*t+pi/3) + b*cos(2*pi*200*t+psi) == cos(2*pi*200*t+pi/3);
[b, psi] = solve(eqn, 'b', 'psi', 'ReturnConditions', true);
152 chars
4 lines