MATLAB: Error :Nonscalar arrays of function handles are not allowed; use cell arrays instead.

errorgraphhelp

Hi everyone,
trying to plot the graph like in the picture but got Error using this code:
Es=200000;  %Mpa



Esh=8500; %Mpa
fy=500;  %Mpa
fsu=750; %Mpa
epssh=0.009;
epssu=0.075;
P=Esh*((epssu-epssh)/(fsu-fy));
epsy=fy/Es;
epscmv = linspace(0.1, 10, 500)*1E-3;
for i=1:numel(epscmv);
    epscm = epscmv(i);
    sigmaSteel(i)=@(epscm) Es*epscm .* (epscm<=epsy) + fy .* (epscm>epsy & epscm<=epssh) + fsu+(fy-fsu)*abs((epssu-epscm)./(epssu-epssh))^(1/P) .* (epscm>epssh & epscm<=epssu) + 0 .* (epscm>epsu);
end
plot(epscmv, sigmaSteel)
grid on
Thank you very much

Best Answer

The loop is not necessary.
This required a few corrections (a typographical error and to vectorise the exponentiation), and now appears to work: 
Es=200000;  %Mpa



Esh=8500; %Mpa
fy=500;  %Mpa
fsu=750; %Mpa
epssh=0.009;
epssu=0.075;
P=Esh*((epssu-epssh)/(fsu-fy));
epsy=fy/Es;
epscmv = linspace(0.1, 10, 500)*1E-3;
sigmaSteel=@(epscm) Es*epscm .* (epscm<=epsy) + fy .* (epscm>epsy & epscm<=epssh) + fsu+(fy-fsu)*abs((epssu-epscm)./(epssu-epssh)).^(1/P) .* (epscm>epssh & epscm<=epssu) + 0 .* (epscm>epssu);
figure
plot(epscmv, sigmaSteel(epscmv))
grid on
It does not exactly reproduce the plot in the 1.jpg attachment, although it appears to be reasonably close.
Related Question