28 Purva Arun Kumbhalkar

1 #include <stdio.h> #include <stdlib.h> #include <math.h> #define f(x) ((x*x*x)-18) int main(){ float a=0,b=0,error=0,m,mold; int i=0; printf("Input Interval: "); scanf("%f %f",&a,&b); if((f(a)*f(b))>0){ printf("Invalid Interval Exit!"); //to test whether search interval exit(1); //is okay or not } else if(f(a)==0 || f(b)==0){ printf("Root is one of interval bounds. Root is %f\n",f(a)==0?a:b); exit(0); } printf("Ite\ta\t\tb\t\tm\t\tf(m)\t\terror\n"); do{ mold=m; m=(a+b)/2; printf("%2d\t%4.6f\t%4.6f\t%4.6f\t%4.6f\t",i++,a,b,m,f(m)); if(f(m)==0){ printf("Root is %4.6f\n",m); }else if ((f(a)*f(m))<0){ b=m; }else a=m; error=fabs(m-mold); if(i==1){ printf("----\n"); }else printf("%4.6f\n",error); }while(error>0.00005); printf("Approximate Root is %4.6f",m); return 0; }

#include<stdio.h> #include<conio.h> #include<math.h> #define SIZE 10 int main() { float a[SIZE][SIZE], x[SIZE],x_new[SIZE]; float temp, lambda_new, lambda_old, error; int i, j,n, step=1; /* Inputs */ printf("Enter Order of Matrix: "); scanf("%d", &n); printf("Enter Tolerable Error: "); scanf("%f", &error); printf("Enter Coefficient of Matrix:\n"); for(i=1;i<n;i++) /* Reading Matrix */ for(j=1;j<=n;j++) { printf("a[%d] [%d]=",i,j); scanf("%f", &a[i]); } } /* Reading Intial Guess Vector */ printf("Enter Initial Guess Vector: \n"); for(i=1;i<=n;i++) printf("x[%d]=",i); scanf("%f",&x[i]); } /* Initializing Lambda_Old */ lambda_old = 1; /* Multiplication */ up: for(i=1;i<=n;i++) temp = 0.0; for(j=1;j<n; j++) { temp temp + a[i]*x0]; } x_new[i] = temp; } /* Replacing */ for(i=1;i<=n;i++) { x[i] = x_new[i]; } lambda_new = fabs(x[1]); /* Finding Largest */ for(i=2;i<n;i++) { if(fabs(x[i])>lambda_new) lambda_new fabs(x[i]); } /* Normalization */ for(i=1;i<=n;i++) x[i]=x[i]/lambda_new; } /* Display */ printf("\n\nITERATION-%d: \n", step); printf("Eigen Value = %f\n", lambda_new); printf("Eigen Vector:\n"); for(i=1;i<=n;i++) printf("%f\t", x[i]); } /* Checking Accuracy */ if(fabs(lambda_new-lambda_old)>error) { lambda_old-lambda_new; step++; goto up; } } return(0);

#include<stdio.h> #include<conio.h> #include<math.h> #include<process.h> #include<string.h> float fun(float); void main() { float result =1; float a,b; float h,sum; int i,j; int n; printf("\n\n Enter the range: "); printf("\n\n Lower limit a:- "); scanf("%f",&a); printf("\n\n upper limit b:- "); scanf("%f",&b); printf("\n\n Enter number of subintervals :- "); scanf("%d",&n); h = (b-a)/n; sum=0; sum-fun(a)+fun(b); for(i=1;i<n;i++) { sum+= 2*fun(a+i); } result sum*h/2; printf("\n\n\n Value of the integral is %6.4f\t", result); printf("\n\n\n Press Enter to Exit"); getch(); } float fun(float x) { float temp; temp = 1/(1*(x*x)); return temp; }