clear all; close all; a = .023; b = 1.5*10^(-7) c = 2120; f = @(t,x) 0.023*x-b*x.^2; [answer1, answer2] = Euler(f, 50, 0, 5, c); plot(answer1,answer2,'r');