How to plot a against varying phi value

I am sorry if this is an easy answer, I am a noob on this.

del=sin(beta)^2(1+s^2+I)-sin(phi)(sin(phi)+u*cos(phi))

tan(beta)=1/s

s=0.95,1.0,1.1,1.2

u=0.75

I=1

I have to plot del against phi

Phi should go from 0 degrees to 90 degrees

How do I go about doing this ?

Thanks in advance.

Answer 1

[Andrei Bobrov]

phi = (0:90)';
s = [0.95,1.0,1.1,1.2];
sd = sind(phi);
plot(phi,bsxfun(@minus,sin(1./s).^2.*(1+s.^2+1),sd.*(sd+.75*cosd(phi)));