The error is on f1(x)
definition. It has boolean result.
The function
f1(x)= (x<=R0 ? rho1 : rho2) && (x>=R ? rhom : rho2)
for x=40 returns
f1(x)= (true) && (false) -> false
which is boolean.
If you would like to plot rho1
for x<=R0
, rho2
for x in (R0,R)
and rhom
for x=>R
then you have to use:
f1(x)= (x<=R0 ? rho1 : (x>=R ? rhom : rho2) )
p [20:80] f1(x)