Ok. You knew you had to use the .^ operator to raise the elements to a power, element-wise. There are also .* and ./ operators, for the same reason. Use them.
And then, we see you still did not even take your own advice, even though you said you knew how to raise sin(x) to a power. You still wrote the wrong expression.
y = sin(x).^3./((1-cos(x)).^5)
__ _
Note the two differences in what I wrote.
You should also learn to use linspace, NOT the colon operator for this. Why?
x = [0:.01:pi];
Look at the final element of x. Is it pi? No. In fact, it is wrong by a fair amount, that in some problems will be significant.
pi - x(end)
ans =
0.00159265358979299
However, if you do this:
x = linspace(0,pi,100);
pi - x(end)
ans =
0
Now it hits the endpoint exactly.
