1 Roots of quadratic equation
File: quadratic.c
#include <stdio.h>
#include <math.h>
void main(){
float a,b,c,d,r1,r2;
printf("\nEnter the coefficients a b c : ");
scanf("%f%f%f" ,&a,&b,&c);
if (a==0) {
printf("\nNot a quadratic equation");
}
d=b*b-4*a*c;
r1=-b/(2*a);
if (d==0){
printf("\nThe roots are real and equal");
printf("\nThe root =%f \n",r1);
}
else if(d>0){
r2=sqrt(d)/(2*a);
printf("\nThe roots are real and distinct");
printf("\nFirst root=%f \nSecond root=%f \n",r1+r2,r1-r2);
}else {
r2= sqrt(fabs(d))/(2*a);
printf("\nThe roots are imaginary");
printf("\nFirst root =%f +i%f", r1,r2);
printf("\nSecond root=%f -i%f \n", r1,r2);
}
}
The program is compiled and executed in terminal as follows:
$ gcc quadratic.c -o quadratic -lm
$ ./quadratic
Enter the coefficients a b c : 1 2 2
The roots are imaginary
First root =-1.000000 +i1.000000
Second root=-1.000000 -i1.000000
2 Rank List
File: ranklist.c
#include <stdio.h>
#include <string.h>
#define NMAX 20
1
#define NAMELENGTH 10
struct studentType{
char name[NAMELENGTH];
int marks;
};
void main(){
int n,i,j,marks;
char name[NAMELENGTH];
struct studentType students[NMAX];
printf("\nEnter number of students (< %d): ",NMAX);
scanf("%d",&n);
for(i=0;i<n;i++){
printf("\nEnter name and marks of sudent#%d: ",i+1);
scanf("%s %d",students[i].name,&students[i].marks);
}
for(i=0;i<n-1;i++)
for(j=0;j<n-1-i;j++){
if(students[j].marks < students[j+1].marks){
strcpy(name,students[j].name);
strcpy(students[j].name,students[j+1].name);
strcpy(students[j+1].name,name);
marks=students[j].marks;
students[j].marks=students[j+1].marks;
students[j+1].marks=marks;
}
}
printf("\nRank List");
printf("\n%7s %-11s %5s","Rank","Name","Marks");
for(i=0;i<n;i++){
printf("\n%7d %-11s %5d",i+1,students[i].name,students[i].marks);
}
printf("\n");
}
Compilation and execution
$ gcc ranklist.c -o ranklist
$ ./ranklist
Enter number of students (< 20): 4
Enter name and marks of sudent#1: alpha 34
Enter name and marks of sudent#2: beta 89
Enter name and marks of sudent#3: gamma 78
2
Enter name and marks of sudent#4: epsilon 96
Rank List
Rank Name Marks
1 epsilon 96
2 beta 89
3 gamma 78
4 alpha 34
3 Print values sin(x), tan(x) and exp(−x)
File name: funcval.c
compilation command: gcc funcval.c -o funcval -lm
exection: ./funcval
#include<stdio.h>
#include<math.h>
#define PI 3.141593
void main() {
float x,xr;
printf("\n enter the value of x in degree :");
scanf("%f" ,&x );
xr= x*PI/180;
printf("\n sine(%f) =%f \t tan(%f) =%f", x,sin(xr),x,tan(xr));
printf("\n e to the power of -%f is %f",x,exp(-x));
}
4 Ploting Functions
To plot a function f (x), we will create a text file with x-values in the first
column and the corresponding f (x) in the second column. In the examples
we discuss in this section, we assume the filename to be ’data.dat’. There are
many softwares available to produce graphs from the data files. Here, we use
the program Gnuplot.
4.1 Gnuplot
This program can be invoked from the terminal using the command gnuplot.
This will start the program with the prompt gnuplot>. The terminal on gnuplot
start-up looks like the following
$ gnuplot
G N U P L O T
Version 4.6 patchlevel 6 last modified September 2014
3
Build System: Linux x86_64
Copyright (C) 1986-1993, 1998, 2004, 2007-2014
Thomas Williams, Colin Kelley and many others
gnuplot home: http://www.gnuplot.info
faq, bugs, etc: type "help FAQ"
immediate help: type "help" (plot window: hit ’h’)
Terminal type set to ’x11’
gnuplot>
Now we discuss some gnuplot-commands that are used in creating two-dimensional
graphs.
command to plot a file:
plot <filename> <xcolumn>:<ycolumn> using <linestyle>
command to set labels of X- and Y- axes:
set xlabel <xlabelname>
set ylabel <ylabelname>
Here, terms within angular brackets < > should be replaced by appropriate
values. An example is provided below.
Example
• Set X-axis label to ‘x’
• Set Y-axis label to ‘f(x)’
• Plot the ‘data.dat’ with x-data from colum 1 and f (x) data from
column 2. Set the legend of the plot-curve to be ‘function f(x)’
The following set of gnuplot commands achieve this.
gnuplot> set xlabel ’x’
gnuplot> set ylabel ’f(x)’
gnuplot> plot ’data.dat’ using 1:2 with linespoints
title ’function f(x)’
The graph will appear on a separate window.
We can also save this graph to an image file. For example, to save the current
graph to a PNG file named ’plot.png’ we use the following set of commands.
gnuplot> set terminal png
Terminal type set to ’png’
Options are ’nocrop font "/usr/share/fonts/truetype/liberation
/LiberationSans-Regular.ttf,12" fontscale 1.0 size 640,480 ’
gnuplot> set output ’sin.png’
gnuplot> replot
We can exit gnuplot using the command exit.
4
gnuplot> exit
4.2 Plot sin(x)
The C program will create a file named ’sin.dat’ with x in the first column and
sin(x) in the second column. Gnuplot can generate the graph from this file.
File: plotsin.c
#include<stdio.h>
#include<math.h>
#define PI 3.141593
void main()
{
FILE *fp;
float xstart,xend,xstep,xr,x;
printf("\n enter initial and final values of x and step in degrees :
");
scanf("%f %f %f",&xstart, &xend, &xstep);
fp=fopen("sin.dat","w");
fprintf(fp,"#x \t sin(x) ");
for(x=xstart;x<=xend;x+=xstep){
xr=x*PI/180;
fprintf(fp,"\n %f %f",x,sin(xr));
}
fclose(fp);
}
Compilation and execution of the program is done the following way.
$ gcc plotsin.c -o plotsin -lm
$ ./plotsin
enter initial and final values of x and step in degrees :
0 360 5
The contents of the file can be viewed using a texteditor, for example gedit or
scratch.
First few lines of the file: sin.dat
#x sin(x)
0.000000 0.000000
5.000000 0.087156
10.000000 0.173648
15.000000 0.258819
20.000000 0.342020
25.000000 0.422618
30.000000 0.500000
35.000000 0.573576
.
.
.
5
The graph can be plotted using the following commands in gnuplot
$ gnuplot
G N U P L O T
Version 4.6 patchlevel 6 last modified September 2014
Build System: Linux x86_64
Copyright (C) 1986-1993, 1998, 2004, 2007-2014
Thomas Williams, Colin Kelley and many others
gnuplot home: http://www.gnuplot.info
faq, bugs, etc: type "help FAQ"
immediate help: type "help" (plot window: hit ’h’)
Terminal type set to ’x11’
gnuplot> set xlabel ’x(degrees)’
gnuplot> set ylabel ’sin(x)’
gnuplot> plot ’sin.dat’ using 1:2 with linespoints
title ’sine function’
The graph looks like
The graph can be saved into a PNG file named ‘sin.png’ using the following
set of commands
gnuplot> set terminal png
Terminal type set to ’png’
Options are ’nocrop font "/usr/share/fonts/truetype/liberation
6
/LiberationSans-Regular.ttf,12" fontscale 1.0 size 640,480 ’
gnuplot> set output ’sin.png’
gnuplot> replot
4.3 Plot tan(x)
File: plottan.c
#include<stdio.h>
#include<math.h>
#define PI 3.141593
void main()
{
FILE *fp;
float xstart,xend,xstep,xr,x;
printf("\n enter initial and final values of x and step in degrees :
");
scanf("%f %f %f",&xstart, &xend, &xstep);
fp=fopen("tan.dat","w");
fprintf(fp,"#x \t tan(x) ");
for(x=xstart;x<=xend;x+=xstep){
xr=x*PI/180;
fprintf(fp,"\n %f %f",x,tan(xr));
}
fclose(fp);
}
Compilation and execution
$ gcc plottan.c -o plottan -lm
$ ./plottan
enter initial and final values of x and step in degrees :
-85 85 5
This creates the file ‘tan.dat’ from which Gnuplot creates the graph
gnuplot> set xlabel ’x’
gnuplot> set ylabel ’tan(x)’
gnuplot> plot ’tan.dat’ u 1:2 with linespoints
title ’tan function’
7
4.4 Plot exp(−x)
File: plotexp.c
#include<stdio.h>
#include<math.h>
void main()
{
FILE *fp;
float x,xstart,xend,xstep;
printf("\n enter initial and final values of x and step : ");
scanf("%f %f %f",&xstart, &xend, &xstep);
fp=fopen("exp.dat","w");
fprintf(fp,"#x \t exp(-x) ");
for(x=xstart;x<=xend;x+=xstep){
fprintf(fp,"\n %f %f",x,exp(-x));
}
fclose(fp);
}
Compilation and execution
$ gcc plotexp.c -o plotexp -lm
$ ./plotexp
enter initial and final values of x and step : 0 10 0.1
8
This creates the file ‘exp.dat’ from which Gnuplot creates the graph
gnuplot> set xlabel ’x’
gnuplot> set ylabel ’exp(-x)’
gnuplot> plot ’exp.dat’ u 1:2 with linespoints
title ’function exp(-x)’
5 Product of two matrices
File: matrix.c
#include<stdio.h>
#include<stdlib.h>
#define ROWS 10
#define COLS 10
void getmatrix(int m, int n, int a[ROWS][COLS]){
int i,j;
for(i=0;i<m;i++)
for(j=0;j<n;j++)
scanf("%d",&a[i][j]);
}
void printmatrix(int m, int n, int a[ROWS][COLS]){
int i,j;
9
for(i=0;i<m;i++){
for(j=0;j<n;j++)
printf("%d ",a[i][j]);
printf("\n");
}
}
void main(){
int a[ROWS][COLS],b[ROWS][COLS],c[ROWS][COLS];
int m,n,p,q,i,j,k;
printf("\n Enter the order of first matrix\n");
scanf("%d%d",&m,&n);
printf("\n Enter the order of the second matrix \n");
scanf("%d%d",&p,&q);
if(n!=p){
printf("\n The orders do not match,cannot multiply\n");
exit(0);
}
printf("\n Enter the elements in first matrix\n");
getmatrix(m,n,a);
printf("\n Enter the elements in second matrix\n");
getmatrix(p,q,b);
for(i=0;i<m;i++){
for(k=0;k<q;k++){
c[i][k]=0;
for(j=0;j<n;j++)
c[i][k]=c[i][k]+a[i][j]*b[j][k];
}
}
printf("\n First matrix \n");
printmatrix(m,n,a);
printf("\n second matrix \n");
printmatrix(p,q,b);
printf("\n The product of the above matrices:\n");
printmatrix(m,q,c);
}
Compilation and execution
$ gcc matrix.c -o matrix
$ ./matrix
Enter the order of first matrix
2 2
Enter the order of the second matrix
10
2 2
Enter the elements in first matrix
3 4
7 8
Enter the elements in second matrix
10 2
1 7
First matrix
3 4
7 8
second matrix
10 2
1 7
The product of the above matrices:
34 34
78 70
6 Scalar and Vector Products
File name: vecproducts.c
compilation command: gcc vecproducts.c -o vecproducts
exection: ./vecproducts
#include <stdio.h>
void main(){
float a1,a2,a3,b1,b2,b3,c1,c2,c3,scalarproduct=0;
printf("\nEnter (3D)vector a:");
scanf("%f %f %f",&a1,&a2,&a3);
printf("\nEnter (3D)vector b:");
scanf("%f %f %f",&b1,&b2,&b3);
scalarproduct=a1*b1+a2*b2+a3*b3;
c1=a2*b3-a3*b2;
c2=a3*b1-a1*b3;
c3=a1*b2-a2*b1;
printf("\nVector a: [%f \t %f \t %f]\n",a1,a2,a3);
printf("\nVector b: [%f \t %f \t %f]\n",b1,b2,b3);
printf("\nScalrproduct= %f\n",scalarproduct);
printf("\nVector product: [%f \t %f \t %f]\n",c1,c2,c3);
}
The compilation and execution in the terminal looks like the following
11
$ gcc vecproducts.c -o vecproducts
$ ./vecproducts
Enter (3D)vector a:1 3 4
Enter (3D)vector b:2 6 -3
Vector a: [1.000000 3.000000 4.000000]
Vector b: [2.000000 6.000000 -3.000000]
Scalrproduct= 8.000000
Vector product: [-33.000000 11.000000 0.000000]
7 Numerical Integration: Trapezoidal Rule
File: trapezoidal.c
#include <stdio.h>
double fn(double x){
return 1/(1+x*x);
}
double trapz(double f(double), double a, double b, int n){
double h, sumt,x;
h=(b-a)/n;
sumt=(f(a)+f(b))/2;
for(x=a+h;x<b;x+=h){
sumt=sumt+f(x);
}
return (sumt*h);
}
void main(){
int n;
double a,b;
printf("\nEnter number of intervals( >2),initial and final points:
");
scanf("%d %lf %lf",&n,&a,&b);
if((n%2)) n++;
printf("\nTrapezoidal rule: n=%d, integral=%lf\n",n,trapz(fn,a,b,n));
}
$ gcc trapezoidal.c -o trapezoidal
$ ./trapezoidal
Enter number of intervals( >2),initial and final points: 20 0 1
12
Trapezoidal rule: n=20, integral=0.785294
8 Numerical Integration: Simpson’s Rule
File: simpson.c
#include <stdio.h>
#include <math.h>
double fn(double x){
return sin(x);
}
double simpson(double f(double), double a, double b, int n){
double h,x,sums;
h=(b-a)/n;
sums=(f(a)+f(b))/2+2*f(b-h);
for(x=a+h;x<b-h;x+=2*h){
sums=sums+2*f(x)+f(x+h);
}
return (sums*h*2/3);
}
void main(){
int n;
double a,b;
printf("\nEnter number of intervals(even, >2),initial and final
points: ");
scanf("%d %lf %lf",&n,&a,&b);
if((n%2)) n++;
printf("\nSimpson rule: n=%d, Integral=%lf\n",n,simpson(fn,a,b,n));
}
Compilation and Execution:
$ gcc simpson.c -o simpson -lm
$ ./simpson
Enter number of intervals(even, >2),initial and final points:
20 0 3.14
Simpson rule: n=20, Integral=2.000006
9 Ordinary Differential Equations: Runge-Kutta
Method
File: rk2.c
13
#include <stdio.h>
float f(float y, float t){
return -y;
}
void main(){
float y,t,yout,h,tstart,tend,k1,k2;
FILE *fp;
y=10;
h=0.1;
tstart=0;
tend=10;
fp=fopen("rk2.dat","w");
fprintf(fp,"\n%f %f",tstart,y);
for(t=tstart;t<tend;t+=h){
k1=f(y,t);
k2=f(y+h*k1,t+h);
yout=y+(k1+k2)*h/2;
y=yout;
fprintf(fp,"\n%f %f",t+h,y);
}
fclose(fp);
}
Compilation and Execution:
$ gcc rk2.c -o rk2
$ ./rk2
Plotting the solution y(t) using Gnuplot:
gnuplot> set xlabel ’t’
gnuplot> set ylabel ’y(t)’
gnuplot> plot ’rk2.dat’ u 1:2 with linespoints title
’solution of dy/dt=-y, y(0)=10’
14
15
