function yp=picard(t,y); %Computes dy/dt=picard(t,y) %Here, yp, t and y are scalars. %But Matlab can handle yp and y being arrays. %The right hand side should satisfy the Lipschitz Condition. %foft(t) ...