Taxicab.f08

program Taxicab
implicit none

    integer (kind = 16) :: taxicabnumber = 48988659276962496
    integer (kind = 16) :: i,j,i3,temp,erotus,raja
    print *,"Ta(3) = 87539319 löytyi v.1957 tietokoneella"
    print *, "Ta(4) = 6963472309248 löytyi v.1991 tietokoneella"
    print *, "Ta(5) = 48988659276962496 löytyi v.1994 tietokoneella"
    print *, "Ta(6) = 24153319581254312065344 löytyi v.2008 tietokoneella"
    read *, taxicabnumber
    print *,"taxicabnumber = ",taxicabnumber
    raja = ceiling(taxicabnumber**(1.0/3.0),8)
    print *,"raja = ",raja
    temp = 0
    do i = 1,raja
        i3 = i*i*i
        erotus = taxicabnumber - i3
        j = ceiling(erotus**(1.0/3.0),8)
        if (i3 + j*j*j == taxicabnumber) then
            if (i == temp) then
                exit
            end if
            print *,"taxicabnumber = ",i,"^3 + ",j,"^3"
            temp = j
        end if
    end do



end program Taxicab