# RK4. Group 1, initialize. f:=(x,y)->1-x-y: x0:=0:y0:=3:h:=0.2:Dots:=[x0,y0]:n:=2:# Group 2, repeat n times. RK4 method. for i from 1 to n do k1:=h*f(x0,y0): k2:=h*f(x0+h/2,y0+k1/2): k3:=h*f(x0+h/2,y0+k2/2): k4:=h*f(x0+h,y0+k3): Y:=y0+(k1+2*k2+2*k3+k4)/6: x0:=x0+h:y0:=Y:Dots:=Dots,[x0,y0]; od: # Group 3, display some dots and plot. P3:=Dots; plot([Dots]);