Unit 5.5 Picard's Method

#include<stdio.h>

#include<conio.h>

#include<math.h>

#define y1(x) (1+pow(x,2)/2-x)

#define y2(x) (1-pow(x,3)/6+pow(x,2)-x

#define y3(x) (1+pow(x,4)/24+pow(x,3)/3-pow(x/2)/2-x)

#define y4(x) (1+pow(x,5)/120-pow(x,5)/120-pow(x,4)/12+pow(x,3)/3-pow(x/2)/2-x)

void main()

{

double y1[20],y2[20],y3[20],y4[20],a,n,h,i;

int j;

clrscr();

printf("enter value of range");

scanf("%lf%lf",&a,&h);

printf("enter the h");

scanf("%lf",&h);

for(i=a,j=0;i<=h;i+i+h,j++)

{

y1[j]=y1[i];

y2[j]=y2[i];

y3[j]=y3[i];

y4[j]=y4[i];

}

printf("|nx|");

for(i=a;i<=n;i=i+h)

{

printf("%3f",i);

printf("y1|");

for(i=a,j=0;i<=n;i=i+n,j++)

{

printf("y=|");

for(i=a,j=0;j<=n;i+i+h,j++)

{

printf("%3f",y2[j]);

}

printf("y3|");

for(i=a,j=0;i<=n;i=i+h,j++)

{

printf("%3f",y3[j]);

}

getch();

}

}

}