af=linspace(0,pi/3,100);
bt=linspace(0,pi/6,100);
[AF,BT]=meshgrid(af,bt);
y=(tan(AF).*(0.5*(sin(AF)+(BT))-cos(AF+BT)))./((sin(BT)+0.5*cos(BT)).*(1+tan(AF).*cot(BT)).*cos(AF));
mesh(AF,BT,y)
clear;
af=pi/3;bt=pi/6;
y=(tan(af).*(0.5*(sin(af)+(bt))-cos(af+bt)))./((sin(bt)+0.5*cos(bt)).*(1+tan(af).*cot(bt)).*cos(af));
subs(y)
ans =