EULER MEJORADO

//Torres Martinez Gustavo Alan
//Grupo:4CM1
//Programa que realiza la solucion de una diferencial ordinaria por el metodo
//de EULER MEJORADO.

#include<stdio.h>
#include<conio.h>
#include<dos.h>
#include<math.h>
#include<stdlib.h>
#include<iostream.h>

void main(void)
{
int j=1;
double nh,h,nx,x=0,y0=2,yp,yc,i,era;
clrscr();
cout<<"La funci¢n es: f(x,y)=xy; sujeta a y0=2.";
cout<<"\n\nDame el valor del incremento h:";
cin>>nh;
cout<<"Dame el valor para x:";
cin>>nx;
for(i=0;i<nx;i=i+nh)
{
yp=y0+(x*y0*nh);
yc=y0+(nh/2)*((x*y0)+((x+nh)*yp));
era=((yc-y0)/yc)*100;if(era<=0.5){era=100;}
x=x+nh;
y0=yc;
printf("\nY%dp=%2.9f   Y%dc=%2.9f   Era(%d)=%2.2f%",j,yp,j,yc,j,fabs(era));
j++;
}
cout<<endl<<endl;
textcolor(10+128);
cprintf("             Presiona cualquier tecla....");
getch();
}