Определить по заданным числам график функции и его формулу

program planer
implicit none
real :: f, l, s1, s2, x, x1, x2, y, y1, y2
 
s1 = 5      !Сторона с одного края
s2 = 3      !Сторона с другого края      
l  = 300    !Полная длина рейки
 
x1 = 1; x2 = l-1; y = 1
 
    do while (abs (y) > 0.1)
        y1 = f(s1, s2, l, x1)
        y2 = f(s1, s2, l, x2)
        x  = (x1 + x2) / 2.0
        y  = f(s1, s2, l, x)           
          if (y1 * y <= 0) then
            x2 = x
          else
            x1 = x
          end if
    end do
 
print '(a, f9.3)', "x = ", x
end program planer
 
!Функция
real function f (s1, s2, l, x)
implicit none
real :: alfa, s1, s2, s3, x, l, v1, v2, h1, h2
   h1 = l - x; h2 = x; alfa = l/(0.5*(s1-s2)); s3 = s1-2*x/alfa
   v1 = h1*(s2**2 + sqrt (s2**2*s3**2) + s3**2)/3.0
   v2 = h2*(s1**2 + sqrt (s1**2*s3**2) + s3**2)/3.0
   f = v1-v2
end function f

program planer
implicit none
real :: x, a1, a2, l, summa1, summa2, s, alfa
 
a1 = 5      !Сторона с одного края
a2 = 3      !Сторона с другого края      
l = 300     !Длина рейки
 
 alfa = l / (0.5*(a1-a2))
 x = 0; summa1 = 0; summa2 = 0
    
  do while (x < l)
    s = (a1 - 2*x/alfa)**2
    summa1 = summa1 + s * x
    summa2 = summa2 + s
    x = x + 0.1
  end do
 
print '(a, f5.1)', "x = ", summa1 / summa2  
end program planer