Logo 
Search:

C Programming Articles

Submit Article
Home » Articles » C Programming » Numerical MethodsRSS Feeds

PICARD'S METHOD

Posted By: Faith Hughes     Category: C Programming     Views: 4407

Write a program of PICARD'S METHOD.

Code for PICARD'S METHOD in C Programming

#include<stdio.h>
#include <math.h>
#include<conio.h>
// dy/dx = 1 + xy#define Y1(x)  (1 + (x) + pow(x,2)/2)
#define Y2(x)  (1 + (x) + pow(x,2)/2 + pow(x,3)/3 + pow(x,4)/8)
#define Y3(x)  (1 + (x) + pow(x,2)/2 + pow(x,3)/3 + pow(x,4)/8 + pow(x,5)/15 + pow(x,6)/48)
void main()
{
  double y1[20],y2[20],y3[20],a,n,h,i;
  int j;
  clrscr();
  printf("\nEnter the value of range: ");
  scanf("%lf %lf",&a,&n);
  printf("\n\nEnter the h: ");
  scanf("%lf",&h);
  for(i=a,j=0;i<=n;i=i+h,j++)
  {

   y1[j]=Y1(i);
   y2[j]=Y2(i);
   y3[j]=Y3(i);
  }
  printf("\nX |");
  for(i=a;i<=n;i=i+h)
  {
   printf("  %.3lf",i);
  }
  printf("\n--------------------------------------------------------------------------------");
  printf("\n\nY1|");
  for(i=a,j=0;i<=n;i=i+h,j++)
  {
   printf("  %.3lf",y1[j]);
  }
  printf("\n\nY2|");
  for(i=a,j=0;i<=n;i=i+h,j++)
  {
   printf("  %.3lf",y2[j]);
  }
  printf("\n\nY3|");
  for(i=a,j=0;i<=n;i=i+h,j++)
  {
   printf("  %.3lf",y3[j]);
  }
getch();
}

/*
OUT PUT
---------

Enter the value of range: 0 1


Enter the h: 0.1


X | 0.000 0.100 0.200 0.300 0.400 0.500 0.600 0.700 0.800 0.900 1.000

--------------------------------------------------------------------------------

Y1| 1.000 1.105 1.220 1.345 1.480 1.625 1.780 1.945 2.120 2.305 2.500


Y2| 1.000 1.105 1.223 1.355 1.505 1.674 1.868 2.089 2.342 2.630 2.958


Y3| 1.000 1.105 1.223 1.355 1.505 1.677 1.874 2.103 2.369 2.680 3.046


*/
  
Share: 


Didn't find what you were looking for? Find more on PICARD'S METHOD Or get search suggestion and latest updates.

Faith Hughes
Faith Hughes author of PICARD'S METHOD is from London, United Kingdom.
 
View All Articles

 
Please enter your Comment

  • Comment should be atleast 30 Characters.
  • Please put code inside [Code] your code [/Code].

 
No Comment Found, Be the First to post comment!