I am not getting the plot for thr following code.Can anyone please identify the mistake?

r=0.6;
P=82.7*10.^6;
a=150*10.^-6;
W=pi*a*2*P;
v=0.3;
for h = [1*10.^-6,4*10.^-6]; 
stress=(3*W*v)/(4*pi*h.^2);
end
plot(h,stress)

Answer 1

[Star Strider]

You need to index ‘h’ and ‘stress’, but it’s easiest to forget the loop and just vectorise it:

r=0.6;
P=82.7*1E^6;
a=150*1E-6;
W=pi*a*2*P;
v=0.3;
h = [1*1E-6,4*10E-6]; 
stress=(3*W.*v)./(4*pi*h.^2);
plot(h,stress)

You probably want to use the linspace function to define ‘h’:

r=0.6;
P=82.7*1E-6;
a=150*1E-6;
W=pi*a*2*P;
v=0.3;
% h = [1*1E-6,4*1E-6]; 
h = linspace(1*1E-6,4*1E-6); 
stress=(3*W.*v)./(4*pi*h.^2);
plot(h,stress)

[Daniel M]

It appears the formula for stress is changed as you have used the dot operator??

[Star Strider]

Yes. The formula is not really changed, simply made compatible with vector variables.