Hi Simon,
The error you're encountering indicates that the function you're trying to find the root of does not cross the x-axis within the interval you've specified. In other words, the function Smax*F(b)*sqrt(b/0.4)-Kic must have different signs (one positive and one negative) at the two endpoints of your interval est_bc for fzero to work. This is a requirement for the bisection method, which is part of the algorithm that fzero uses to find zeros.
To resolve the error, you need to ensure that the interval you provide to fzero contains a sign change. Here are some steps you can take:
- Evaluate your function at the endpoints of your interval to make sure one is positive and the other is negative.
- If both function values have the same sign, you need to adjust your interval. You can either widen the interval or choose different starting points based on your knowledge of the problem.
- If possible, plot the function over a range of b values to visually inspect where the sign change might occur. This can help you choose a better interval for fzero.