Write a C program on Numerical solution of ODE using Euler method by MistarAV

 #include<stdio.h>

#include<conio.h>

#define f(x,y) x+y

int main()

{

clrscr();

float x0, y0, xn, h, yn, slope;

int i, n;

printf("Enter Initial Condition\n");

printf("x0 = ");

scanf("%f", &x0);

printf("y0 = ");

scanf("%f", &y0);

printf("Enter calculation point xn = ");

scanf("%f", &xn);

printf("Enter number of steps: ");

scanf("%d", &n);

/* Calculating step size (h) */

h = (xn-x0)/n;

/* Euler's Method */

printf("\nx0\ty0\tslope\tyn\n");

printf("------------------------------\n");

for(i=0; i < n; i++)

{

slope = f(x0, y0);

yn = y0 + h * slope;

printf("%.4f\t%.4f\t%0.4f\t%.4f\n",x0,y0,slope,yn);

y0 = yn;

x0 = x0+h;

}

/* Displaying result */

printf("\nValue of y at x = %0.2f is %0.3f",xn, yn);

getch();

}





Output

Enter Initial Condition

x0 = 0

y0 = 1

Enter calculation point xn = 1

Enter number of steps: 10

x0 y0 slope yn

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

0.0000 1.0000 1.0000 1.1000

0.1000 1.1000 1.2000 1.2200

0.2000 1.2200 1.4200 1.3620

0.3000 1.3620 1.6620 1.5282

0.4000 1.5282 1.9282 1.7210

0.5000 1.7210 2.2210 1.9431

0.6000 1.9431 2.5431 2.1974

0.7000 2.1974 2.8974 2.4872

0.8000 2.4872 3.2872 2.8159

0.9000 2.8159 3.7159 3.1875

Value of y at x = 1.00 is 3.187

Comments

Popular posts from this blog

Google Apprenticeship 2023 scheme

SILENT CRYPTOMINER v2.0.1 - MINER FOR ETH ― ETC ― XMR & MORE BY MISTARAV

How to install VS Code in an Android Phone by Termux ?