age = unifReal(45,90);
bpSystolic = unifReal(80,200);
bpSystolicError = bpSystolic + unifReal(-20,20);
hypTension= unifReal(0,1);
hypTensionErr = hypTension + unifReal(-0.05,0.05);
bmi=unifReal(20,40);
bmiErr = bmi+unifReal(-4,4);
prInterval =unifReal(100,250);
prIntervalErr = prInterval + unifReal(-10,10);
sigMurmur = unifReal(0.0,0.55);
if (sigMurmur <= 0.5)
then
   hFail = unifReal(0.0,0.52)
else
   hFail = unifReal (0.0,0.7)
end;
sigMurmurErr = sigMurmur + unifReal(-0.01, 0.01);
hFailErr = hFail + unifReal(-0.01, 0.01);

points = 0;
pointsErr= 0;
if (age <= 49.0)
then 
   points = points - 3
else
   if (age <= 54.0)
   then 
      points = points - 2
   else 
       if (age <= 59.0)
       then 
          points = points + 0
       else 
           if (age <= 64.0)
           then 
             points = points + 1
           else
              if (age <= 69.0) 
              then 
                points = points + 3
              else
                if (age <= 74.0)
                then
                  points = points + 4 
                else
                  if (age <= 79.0)
                  then 
                     points = points + 6
                  else
                     points = points + 7 
                  end 
                end
              end 
           end 
       end
   end
end;

pointsErr = points;

if (bpSystolic >= 160.0)
then 
     points = points + 1
end;
if (bpSystolicError >= 160.0)
then 
     pointsErr = pointsErr + 1
end;
if (hypTension >= 0.5)
then 
     points = points + 1
end;

if (hypTensionErr >= 0.5)
then 
     pointsErr = pointsErr + 1
end;

if (bmi >= 30.0)
then 
   points = points + 1
end;

if (bmiErr >= 30.0)
then 
     pointsErr = pointsErr + 1
end;
  
if (prInterval >= 160.0)
then 
     if (prInterval <= 199.0)
     then 
     	  points = points + 1
     else
	  points = points + 2
     end
end;


if (prIntervalErr >= 160.0)
then 
     if (prIntervalErr <= 199.0)
     then 
     	  pointsErr = pointsErr + 1
     else
	  pointsErr = pointsErr + 2
     end
end;

if (age <= 54.0)
then 
   if (sigMurmur >= 0.5)
   then	
   	points = points + 5
   end;
   if (sigMurmurErr >= 0.5)
   then 
   	pointsErr = pointsErr + 5
   end;
   if (hFail >= 0.5)
   then 
   	points = points + 10
   end;
   if (hFailErr >= 0.5)
   then 
   	pointsErr = pointsErr + 10
   end
   
else
   if (age <= 64.0)
   then
	if (sigMurmur >= 0.5)
   	then	
   	   points = points + 4
   	end;
   	if (sigMurmurErr >= 0.5)
   	then 
   	   pointsErr = pointsErr + 4
        end; 
  	if (hFail >= 0.5)
        then 
   	  points = points + 6
        end;
        if (hFailErr >= 0.5)
        then 
   	  pointsErr = pointsErr + 6
        end

    else 
       if (age <= 74.0) 
       then 
       	if (sigMurmur >= 0.5)
   	then	
   	   points = points + 2
   	end;
   	if (sigMurmurErr >= 0.5)
   	then 
   	   pointsErr = pointsErr + 2
        end; 
  	if (hFail >= 0.5)
        then 
   	  points = points + 2
        end;
        if (hFailErr >= 0.5) 
        then 
   	  pointsErr = pointsErr + 2
        end
       else
         if (age <= 84.0)
	 then 
	   if (sigMurmur >= 0.5)
   	   then	
   	     points = points + 1
   	   end;
   	   if (sigMurmurErr >= 0.5)
   	   then 
   	     pointsErr = pointsErr + 1
           end   
	 end
       end
   end

end;

estimateProb(points >= 10);
estimateProb(points - pointsErr >= 5);
estimateProb(pointsErr - points <= -5 )