All constants need to be defined. Then use:
t = pi/2:0.01:2*pi+(pi/2);
phi1 = (a^2+a^2+d^2+2*a*d*cos(pi-t)-d^2)./(2*b*sqrt(a^2+d^2+2*a*d*cos(pi-t)));
plot(t, phi1)
t is a vector, therefore, you need to use the ./ operator.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!