Yes, you're right. Adding the quotes worked.

This seems like a big fat bug to me.

Yes, you're right. Adding the quotes worked.

This seems like a big fat bug to me.

I had previously tried applying evalf to rr and cb, and that doesn't work. I get the same error.

Using is(cb > rr)=true eliminates the error, but the plot still isn't working. Instead of the curve I'm expecting I get a flat line where y=0 for all b (0..Pi/2). That's not right.

If I call occultRatio(0, 2, 2) by itself I get the result 1/6.

However if I call plot(occultRatio(b, 2, 2), b=0..Pi/2) I get a flat line at zero. Not right.

If I call occultRatio(b, JupiterRadius, JupiterRadius) outside of plot where JupiterRadius has units [km] then I get this error:

Error, (in Units:-Standard:-+) the units `1` and `km` have incompatible dimensions

That makes no sense because the units are supposed to cancel out in the expression in rr (which is where I assume this is occurring).

I had previously tried applying evalf to rr and cb, and that doesn't work. I get the same error.

Using is(cb > rr)=true eliminates the error, but the plot still isn't working. Instead of the curve I'm expecting I get a flat line where y=0 for all b (0..Pi/2). That's not right.

If I call occultRatio(0, 2, 2) by itself I get the result 1/6.

However if I call plot(occultRatio(b, 2, 2), b=0..Pi/2) I get a flat line at zero. Not right.

If I call occultRatio(b, JupiterRadius, JupiterRadius) outside of plot where JupiterRadius has units [km] then I get this error:

Error, (in Units:-Standard:-+) the units `1` and `km` have incompatible dimensions

That makes no sense because the units are supposed to cancel out in the expression in rr (which is where I assume this is occurring).

Right. I guess the question should be why does using units cause this to fail?

I plan to use functions like this in a much larger system where units matter.

Can I still use units but isolate the function in a way that allows me to use it?

Right. I guess the question should be why does using units cause this to fail?

I plan to use functions like this in a much larger system where units matter.

Can I still use units but isolate the function in a way that allows me to use it?