function[] = prob9()

    f = @(t,x) 0.023*x-(1.5*10^(-7))*t^2;
    Euler(f, 50, 0, 5, 2120);