In many engineering situations, we are interested in finding the correlation ρ between different quantities x and y based on the values xi and yi of these quantities measured in different situations i. Measurements are never absolutely accurate; it is therefore necessary to take this inaccuracy into account when estimating the correlation ρ. Sometimes, we know the probabilities of different values of measurement errors, but in many cases, we only know the upper bounds Δxi and Δyi on the corresponding measurement errors. In such situations, after we get the measurement results Xi and Yi, the only information that we have about the actual (unknown) values xi and yi is that they belong to the corresponding intervals [Xi - Δxi, Xi + Δxi] and [Yi - Δyi, Yi + Δyi]. Different values from these intervals lead, in general, to different values of the correlation ρ. It is therefore desirable to find the range [ρ-, ρ+] of possible values of the correlation when xi and yi take values from the corresponding intervals. In general, the problem of computing this range is NP-hard. In this paper, we provide a feasible (= polynomial-time) algorithm for computing at least one of the endpoints of this interval: for computing ρ+ when ρ+ > 0 and for computing ρ- when ρ- < 0.