tipo_de_macro=R2
id='x0' expresi&oacute;n='0' evaluar='una-sola-vez' constante='si'
id='y0' expresi&oacute;n='0' evaluar='una-sola-vez' constante='si'
id='t0' expresi&oacute;n='0' evaluar='una-sola-vez' constante='si'
id='dt' expresi&oacute;n='0.01' evaluar='una-sola-vez' constante='si'
id='T' expresi&oacute;n='30' evaluar='una-sola-vez' constante='si'
id='N' expresi&oacute;n='ent(abs(T)/abs(dt))'
id='X' tamao='N+1' evaluar='siempre' vector='si'
id='Y' tamao='N+1' evaluar='siempre' vector='si'
id='fx(x,y,t)' expresi&oacute;n='0'
id='fy(x,y,t)' expresi&oacute;n='0'
inicio='X[0]=x0;Y[0]=y0;i=0;t=t0;dt=(T>0)?abs(dt):-abs(dt)' hacer='p=X[i];q=Y[i];k1x=fx(p,q,t)*dt;k1y=fy(p,q,t)*dt;t=t+dt/2;k2x=fx(p+k1x/2,q+k1y/2,t)*dt;k2y=fy(p+k1x/2,q+k1y/2,t)*dt;k3x=fx(p+k2x/2,q+k2y/2,t)*dt;k3y=fy(p+k2x/2,q+k2y/2,t)*dt;t=t+dt/2;k4x=fx(p+k3x,q+k3y,t)*dt;k4y=fy(p+k3x,q+k3y,t)*dt;i=i+1;X[i]=p+(k1x+2*k2x+2*k3x+k4x)/6;Y[i]=q+(k1y+2*k2y+2*k3y+k4y)/6;' mientras='(i<=N)&&(t-t0<T)' id='CALCULAR' evaluar='siempre' algoritmo='si'