Higher Engineering Mathematics - B. S. Grewal Companion Text

Scilab Textbook Companion for Higher Engineering Mathematics by B. S.
Grewal1
Created by Karan Arora and Kush Garg B.Tech. (pursuing) Civil Engineering Indian Institute of Technology Roorkee College Teacher Self Cross-Checked by Santosh Kumar, IIT Bombay August 10, 2013
by a grant from the National Mission on Education through ICT, http://spoken-tutorial.org/NMEICT-Intro. This Textbook Companion and Scilab codes written in it can be downloaded from the Textbook Companion Project section at the website http://scilab.in
1 Funded
Book Description
Title: Higher Engineering Mathematics Author: B. S. Grewal Publisher: Khanna Publishers, New Delhi Edition: 40 Year: 2007 ISBN: 8174091955
Scilab numbering policy used in this document and the relation to the above book. Exa Example (Solved example) Eqn Equation (Particular equation of the above book) AP Appendix to Example(Scilab Code that is an Appednix to a particular Example of the above book) For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means a scilab code whose theory is explained in Section 2.3 of the book.
Contents
List of Scilab Codes 1 Solution of equation and curve tting 2 Determinants and Matrices 4 Dierentiation and Applications 5 Partial Dierentiation And Its Applications 6 Integration and its Applications 9 Innite Series 10 Fourier Series 13 Linear Dierential Equations 21 Laplace Transform 22 Integral Transform 23 Statistical Methods 24 Numerical Methods 26 Dierence Equations and Z Transform 27 Numerical Solution of Ordinary Dierential Equations 5 15 25 40 57 61 69 74 85 94 108 111 124 134 142
28 Numerical Solution of Partial Dierential Equations 34 Probability and Distributions 35 Sampling and Inference
161 171 189
List of Scilab Codes

Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa 1.1 1.2 1.3 1.6 1.7 1.11 1.12 1.13 1.14 1.15 1.16 1.17 1.18 1.19 1.20 1.21 1.22 1.23 1.24 1.25 2.1 2.2 2.3 2.4 5.8 2.16 2.17 2.18 nding the roots of quadratic equations . . . . . . . nding the roots of equation containing one variable nding the roots of equation containing one variable nding the roots of equation containing one variable nding the roots of equation containing one variable forming an equation with known roots . . . . . . . . forming an equation under restricted conditions . . . nding the roots of equation containing one variable nding the roots of equation containing one variable nding the roots of equation containing one variable nding the roots of equation containing one variable nding the roots of equation containing one variable Finding the roots of equation containing one variable Finding the roots of equation containing one variable Finding the roots of equation containing one variable Finding the roots of equation containing one variable Finding the roots of equation containing one variable Finding the solution of equation by drawing graphs . Finding the solution of equation by drawing graphs . Finding the solution of equation by drawing graphs . Calculating Determinant . . . . . . . . . . . . . . . Calculating Determinant . . . . . . . . . . . . . . . Calculating Determinant . . . . . . . . . . . . . . . Calculating Determinant . . . . . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . product of two matrices . . . . . . . . . . . . . . . . Product of two matrices . . . . . . . . . . . . . . . . Product and inverse of matrices . . . . . . . . . . . . 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 15 16 16 16 17 17 18 18 19 19 19 20 20 20 21 21 21 22 23 25 25 26 26 26 27 27 27
Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa
2.19 2.20 2.23 2.24.1 2.24.2 2.25 2.26 2.28 2.31 2.32 2.34.1 2.34.2 2.38 2.39 2.42 2.43 2.44 2.45 2.46 2.47 2.48 2.49 2.50 2.51 2.52 2.53 2.54 4.4.1 4.5 4.6 4.7 4.8 4.9 4.10 4.11 4.16 4.17 4.18
Solving equation of matrices . . . . . . . . . . . . Nth power of a given matrix . . . . . . . . . . . . Inverse of matrix . . . . . . . . . . . . . . . . . . . Rank of a matrix . . . . . . . . . . . . . . . . . . . Rank of a matrix . . . . . . . . . . . . . . . . . . . Inverse of matrix . . . . . . . . . . . . . . . . . . . eigen values vectors rank of matrix . . . . . . . . . Inverse of a matrix . . . . . . . . . . . . . . . . . . Solving equation using matrices . . . . . . . . . . . Solving equation using matrices . . . . . . . . . . . predicting nature of equation using rank of matrix predicting nature of equation using rank of matrix Inverse of a matrix . . . . . . . . . . . . . . . . . . Transpose and product of matrices . . . . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and characteristic equation . . . . . . eigen values and characteristic equation . . . . . . eigen values and characteristic equation . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and vectors of given matrix . . . . . . eigen values and vectors of given matrix . . . . . . Hermitian matrix . . . . . . . . . . . . . . . . . . . tranpose and inverse of complex matrix . . . . . . Unitary matrix . . . . . . . . . . . . . . . . . . . . nding nth derivative . . . . . . . . . . . . . . . . nding nth derivative . . . . . . . . . . . . . . . . nding nth derivative . . . . . . . . . . . . . . . . nding nth derivative . . . . . . . . . . . . . . . . proving the given dierential equation . . . . . . . proving the given dierential equation . . . . . . . proving the given dierential equation . . . . . . . verify roles theorem . . . . . . . . . . . . . . . . . expansion using maclaurins series . . . . . . . . . . expanding function as fourier series of sine term . . expansion using maclaurins series . . . . . . . . . . 6
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
28 28 28 29 29 29 29 30 30 30 31 31 32 32 32 33 33 34 35 35 36 36 37 37 37 38 38 40 40 41 42 42 43 44 45 46 46 47
Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa
4.19 4.20 4.21 4.22 4.32 4.34 4.35 4.37 4.39 4.41 4.43 4.46 4.47 4.48 4.52 4.61 5.5 5.14 5.25.1 5.25.2 5.25.3 5.26 5.30 6.1.1 6.1.2 6.2.1 6.2.2 4.2.3 6.2.3 6.4.1 4.4.2 6.5 6.6.1 6.7.1 6.8 6.9.1 6.9.2 6.10
expansion using maclaurins series . . . . . . . . . . . . expansion using taylors series . . . . . . . . . . . . . . taylor series . . . . . . . . . . . . . . . . . . . . . . . . evaluating limit . . . . . . . . . . . . . . . . . . . . . . tangent to curve . . . . . . . . . . . . . . . . . . . . . nding equation of normal . . . . . . . . . . . . . . . . nding angle of intersection of curve . . . . . . . . . . prove given tangent statement . . . . . . . . . . . . . nding angle of intersection of curve . . . . . . . . . . nding pedal equation of parabola . . . . . . . . . . . nding radius of curvature of cycloid . . . . . . . . . . radius of curvature of cardoid . . . . . . . . . . . . . . cordinates of centre of curvature . . . . . . . . . . . . proof statement cycloid . . . . . . . . . . . . . . . . . maxima and minima . . . . . . . . . . . . . . . . . . . nding the asymptotes of curve . . . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . Partial derivative of given function . . . . . . . . . . . indenite integral . . . . . . . . . . . . . . . . . . . . indenite integral . . . . . . . . . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . Denite Integration of a function . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . reducing indenite integral to simpler form . . . . . . Indenite Integration of a function . . . . . . . . . . . Getting the manual input of a variable and integration Denite Integration of a function . . . . . . . . . . . . Denite Integration of a function . . . . . . . . . . . . denite integral . . . . . . . . . . . . . . . . . . . . . . 7
47 48 48 49 50 50 51 51 52 53 53 54 54 55 55 55 57 57 58 58 59 59 60 61 61 61 62 62 62 63 63 63 64 64 65 65 65 65
Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa
6.12 6.13 6.14 6.15 6.16 6.24 9.1 9.1.3 9.2.1 9.2.2 9.5.1 9.5.2 9.7.1 9.7.3 9.8.1 9.8.2 9.10.1 9.10.2 9.11.1 9.11.2 10.1 10.2 10.3 10.4 10.5
Exa 10.6 Exa 10.7 Exa 10.8 Exa Exa Exa Exa Exa Exa 10.9 10.10 10.11 10.12 10.13 10.14
Denite Integration of a function . . . . . . . . . . . . sum of innite series . . . . . . . . . . . . . . . . . . . nding the limit of the function . . . . . . . . . . . . . Denite Integration of a function . . . . . . . . . . . . Denite Integration of a function . . . . . . . . . . . . Calculating the area under two curves . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the sum of series upto innity . . . . . . . . . . to check for the type of series . . . . . . . . . . . . . . to check the type of innite series . . . . . . . . . . . . to check the type of innite series . . . . . . . . . . . . to check the type of innite series . . . . . . . . . . . . to check the type of innite series . . . . . . . . . . . . to nd the sum of series upto innity . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . to nd the limit at innity . . . . . . . . . . . . . . . nding fourier series of given function . . . . . . . . . nding fourier series of given function . . . . . . . . . nding fourier series of given function . . . . . . . . . nding fourier series of given function . . . . . . . . . nding fourier series of given function in interval minus pi to pi . . . . . . . . . . . . . . . . . . . . . . . . . . nding fourier series of given function in interval minus l to l . . . . . . . . . . . . . . . . . . . . . . . . . . . . nding fourier series of given function in interval minus pi to pi . . . . . . . . . . . . . . . . . . . . . . . . . . nding fourier series of given function in interval minus pi to pi . . . . . . . . . . . . . . . . . . . . . . . . . . nding half range sine series of given function . . . . . nding half range cosine series of given function . . . . expanding function as fourier series of sine term . . . . nding fourier series of given function . . . . . . . . . nding complex form of fourier series . . . . . . . . . . practical harmonic analysis . . . . . . . . . . . . . . . 8
66 66 66 67 67 67 69 69 69 70 70 70 71 71 71 72 72 72 72 73 74 74 75 75 76 77 77 78 78 79 80 80 81 81
Exa 10.15 practical harmonic analysis . . . . Exa 10.16 practical harmonic analysis . . . . Exa 10.17 practical harmonic analysis . . . . Exa 13.1 solvinf linear dierential equation . Exa 13.2 solving linear dierential equation Exa 13.3 solving linear dierential equation Exa 13.4 solving linear dierential equation Exa 13.5 nding particular integral . . . . . Exa 13.6 nding particular integral . . . . . Exa 13.7 nding particular integral . . . . . Exa 13.8 nding particular integral . . . . . Exa 13.9 nding particular integral . . . . . Exa 13.10 nding particular integral . . . . . Exa 13.11 solving the given linear equation . Exa 13.12 solving the given linear equation . Exa 13.13 solving the given linear equation . Exa 13.14 solving the given linear equation . Exa 21.1.1 nding laplace transform . . . . . Exa 21.1.2 nding laplace transform . . . . . Exa 21.1.3 nding laplace transform . . . . . Exa 21.2.1 nding laplace transform . . . . . Exa 21.2.2 nding laplace transform . . . . . Exa 21.2.3 nding laplace transform . . . . . Exa 21.4.1 nding laplace transform . . . . . Exa 21.4.2 nding laplace transform . . . . . Exa 21.5 nding laplace transform . . . . . Exa 21.7 nding laplace transform . . . . . Exa 21.8.1 nding laplace transform . . . . . Exa 21.8.2 nding laplace transform . . . . . Exa 21.8.3 nding laplace transform . . . . . Exa 21.8.4 nding laplace transform . . . . . Exa 21.9.1 nding laplace transform . . . . . Exa 21.9.2 nding laplace transform . . . . . Exa 21.10.1nding laplace transform . . . . . Exa 21.10.3nding laplace transform . . . . . Exa 21.11.1nding inverse laplace transform . Exa 21.11.2nding inverse laplace transform . Exa 21.12.1nding inverse laplace transform . 9
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
82 82 83 85 85 86 86 87 87 88 88 89 89 90 90 91 92 94 94 94 95 95 95 96 96 96 97 97 97 98 98 98 99 99 99 100 100 100
Exa 21.12.3nding inverse laplace transform . . . . . . . . . . . . Exa 21.13.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.13.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.14.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.14.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.15.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.15.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.16.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.16.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.16.3nding inverse laplace transform . . . . . . . . . . . . Exa 21.17.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.17.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.19.1nding inverse laplace transform . . . . . . . . . . . . Exa 21.19.2nding inverse laplace transform . . . . . . . . . . . . Exa 21.28.1nding laplace transform . . . . . . . . . . . . . . . . Exa 21.28.2nding laplace transform . . . . . . . . . . . . . . . . Exa 21.34 nding laplace transform . . . . . . . . . . . . . . . . Exa 22.1 nding fourier sine integral . . . . . . . . . . . . . . . Exa 22.2 nding fourier transform . . . . . . . . . . . . . . . . . Exa 22.3 nding fourier transform . . . . . . . . . . . . . . . . . Exa 22.4 nding fourier sine transform . . . . . . . . . . . . . . Exa 22.5 nding fourier cosine transform . . . . . . . . . . . . . Exa 22.6 nding fourier sine transform . . . . . . . . . . . . . . Exa 23.1 Calculating cumulative frequencies of given using iterations on matrices . . . . . . . . . . . . . . . . . . . . . Exa 23.2 Calculating mean of of statistical data performing iterations matrices . . . . . . . . . . . . . . . . . . . . . . Exa 23.3 Analysis of statistical data performing iterations on matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . Exa 23.4 Analysis of statistical data . . . . . . . . . . . . . . . Exa 23.5 Finding the missing frequency of given statistical data using given constants . . . . . . . . . . . . . . . . . . Exa 23.6 Calculating average speed . . . . . . . . . . . . . . . . Exa 23.7 Calculating mean and standard deviation performing iterations on matrices . . . . . . . . . . . . . . . . . . . Exa 23.8 Calculating mean and standard deviation performing iterations on matrices . . . . . . . . . . . . . . . . . . .
101 101 101 102 102 103 103 103 104 104 104 105 105 106 106 106 107 108 108 109 109 109 110 111 112 113 114 114 115 115 117
10
Exa 23.9
Analysis of statistical data performing iterations on matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . Exa 23.10 Calculating mean and standard deviation of dierent statistical data when put together . . . . . . . . . . . Exa 23.12 Calculating median and quartiles of given statistical data performing iterations on matrices . . . . . . . . . . . . Exa 23.13 Calculating coecient of correlation . . . . . . . . . . Exa 24.1 nding the roots of equation . . . . . . . . . . . . . . Exa 24.3 nding the roots of equation by the method of false statement . . . . . . . . . . . . . . . . . . . . . . . . . Exa 24.4 nding rea roots of equation by regula falsi method . . Exa 24.5 real roots of equation by newtons method . . . . . . . Exa 24.6 real roots of equation by newtons method . . . . . . . Exa 24.7 evaluating square root by newtons iterative method . . Exa 24.10 solving equations by guass elimination method . . . . Exa 24.12 solving equations by guass elimination method . . . . Exa 24.13 solving equations by guass elimination method . . . . Exa 26.2 nding dierence equation . . . . . . . . . . . . . . . . Exa 26.3 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.4 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.6 rming bonacci dierence equation . . . . . . . . . . Exa 26.7 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.8 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.10 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.11 solving dierence equation . . . . . . . . . . . . . . . . Exa 26.12 solving simultanious dierence equation . . . . . . . . Exa 26.15.2Z transform . . . . . . . . . . . . . . . . . . . . . . . . Exa 26.16 evaluating u2 and u3 . . . . . . . . . . . . . . . . . . . Exa 27.1 solving ODE with picards method . . . . . . . . . . . Exa 27.2 solving ODE with picards method . . . . . . . . . . . Exa 27.5 solving ODE using Eulers method . . . . . . . . . . . Exa 27.6 solving ODE using Eulers method . . . . . . . . . . . Exa 27.7 solving ODE using Modied Eulers method . . . . . . Exa 27.8 solving ODE using Modied Eulers method . . . . . . Exa 27.9 solving ODE using Modied Eulers method . . . . . . Exa 27.10 solving ODE using runge method . . . . . . . . . . . . Exa 27.11 solving ODE using runge kutta method . . . . . . . . Exa 27.12 solving ODE using runge kutta method . . . . . . . . 11
118 120 120 121 124 125 125 126 127 128 128 130 132 134 135 135 136 136 137 138 138 139 140 140 142 142 143 144 144 145 146 147 148 148
Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa Exa
solving ODE using runge kutta method . . . . . . . . solving ODE using milnes method . . . . . . . . . . . solving ODE using runge kutta and milnes method . . solving ODE using adamsbashforth method . . . . . . solving ODE using runge kutta and adams method . . solving simultanious ODE using picards method . . . solving ssecond ODE using runge kutta method . . . . solving ODE using milnes method . . . . . . . . . . . classication of partial dierential equation . . . . . . solving elliptical equation . . . . . . . . . . . . . . . . evaluating function satisfying laplace equation . . . . . solution of poissons equation . . . . . . . . . . . . . . solving parabolic equation . . . . . . . . . . . . . . . . solving heat equation . . . . . . . . . . . . . . . . . . solving wave equation . . . . . . . . . . . . . . . . . . solving wave equation . . . . . . . . . . . . . . . . . . Calculating probability . . . . . . . . . . . . . . . . . Calculating the number of permutations . . . . . . . . Number of permutations . . . . . . . . . . . . . . . . . Calculating the number of committees . . . . . . . . . Finding the number of committees . . . . . . . . . . . Finding the number of committees . . . . . . . . . . . Finding the probability of getting a four in a single throw of a die . . . . . . . . . . . . . . . . . . . . . . . Exa 34.4.2 Finding the probability of getting an even number in a single throw of a die . . . . . . . . . . . . . . . . . . . Exa 34.5 Finding the probability of 53 sundays in a leap year . Exa 34.6 probability of getting a number divisible by 4 under given conditions . . . . . . . . . . . . . . . . . . . . . Exa 34.7 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.8 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.9.1 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.9.2 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.9.3 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.13 probability of drawing an ace or spade from pack of 52 cards . . . . . . . . . . . . . . . . . . . . . . . . . . . Exa 34.14.1Finding the probability . . . . . . . . . . . . . . . . . Exa 34.15.1Finding the probability . . . . . . . . . . . . . . . . . 12
27.13 27.14 27.15 27.16 27.17 27.18 27.19 27.20 28.1 28.2 28.3 28.4 28.5 28.6 28.7 28.8 34.1 34.2.1 34.2.2 34.3.1 34.3.2 34.3.3 34.4.1
149 150 152 154 155 157 159 160 161 161 162 163 165 166 168 169 171 171 171 172 172 172 173 173 173 174 174 175 175 176 176 177 177 177
Exa 34.15.2Finding the probability . . . . . . . . . . . . . . . . . Exa 34.15.3Finding the probability . . . . . . . . . . . . . . . . . Exa 34.16 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.17 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.18 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.19.1Finding the probability . . . . . . . . . . . . . . . . . Exa 34.19.2Finding the probability . . . . . . . . . . . . . . . . . Exa 34.19.3Finding the probability . . . . . . . . . . . . . . . . . Exa 34.20 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.22 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.23 Finding the probability . . . . . . . . . . . . . . . . . Exa 34.25 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.26 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.27 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.28 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.29 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.30 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.31 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.33 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.34 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.35 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.38 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.39 nding the probability . . . . . . . . . . . . . . . . . . Exa 34.40 nding the probability . . . . . . . . . . . . . . . . . . Exa 35.1 calculating the SD of given sample . . . . . . . . . . . Exa 35.2 Calculating SD of sample . . . . . . . . . . . . . . . . Exa 35.3 Analysis of sample . . . . . . . . . . . . . . . . . . . . Exa 35.4 Analysis of sample . . . . . . . . . . . . . . . . . . . . Exa 35.5 Checking whether real dierence will be hidden . . . . Exa 35.6 Checking whether given sample can be regarded as a random sample . . . . . . . . . . . . . . . . . . . . . . Exa 35.9 Checking whethet samples can be regarded as taken from the same population . . . . . . . . . . . . . . . . Exa 35.10 calculating SE of dierence of mean hieghts . . . . . . Exa 35.12 Mean and standard deviation of a given sample . . . . Exa 35.13 Mean and standard deviation of a given sample . . . . Exa 34.15 Standard deviation of a sample . . . . . . . . . . . . .
178 178 178 179 179 179 180 180 181 181 181 182 182 183 183 183 184 185 185 186 186 187 187 188 189 189 190 191 191 192 192 193 193 194 195
13
List of Figures
1.1 1.2 1.3 6.1 Finding the solution of equation by drawing graphs . . . . . Finding the solution of equation by drawing graphs . . . . . Finding the solution of equation by drawing graphs . . . . . Calculating the area under two curves . . . . . . . . . . . . . 22 23 24 68
14
Chapter 1 Solution of equation and curve tting
Scilab code Exa 1.1 nding the roots of quadratic equations

1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =2*( x ^3) + x ^2 -13* x +6 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
Scilab code Exa 1.2 nding the roots of equation containing one variable
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =3*( x ^3) -4*( x ^2) + x +88 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
15
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^3 -7*( x ^2) +36 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^4 -2*( x ^3) -21*( x ^2) +22* x +40 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =2*( x ^4) -15*( x ^3) +35*( x ^2) -30* x +8 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
16
Scilab code Exa 1.11 forming an equation with known roots

1 2 3 4 5 6 7 8 9 10 11 12 13 14
clear clc x = poly ([0] , x ) ; x1 = poly ([0] , x1 ) ; x2 = poly ([0] , x2 ) ; x3 = poly ([0] , x3 ) ; p = x ^3 -3*( x ^2) +1 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) roots ( p ) disp ( l e t ) x1 =0.6527036 x2 = -0.5320889 x3 =2.8793852 disp ( s o t h e e q u a t i o n whose r o o t s a r e c u b e o f t h e r o o t s o f a b o v e e q u a t i o n i s ( xx1 3 ) ( xx2 3 ) ( x x3 3 ) =0 => ) 15 p1 =( x - x1 ^3) *( x - x2 ^3) *( x - x3 ^3)
Scilab code Exa 1.12 forming an equation under restricted conditions

1 2 3 4 5 6 7 8 9 10 11 12 13
clear clc x = poly ([0] , x ) ; x1 = poly ([0] , x1 ) ; x2 = poly ([0] , x2 ) ; x3 = poly ([0] , x3 ) ; x4 = poly ([0] , x4 ) ; x5 = poly ([0] , x5 ) ; x6 = poly ([0] , x6 ) ; p = x ^3 -6*( x ^2) +5* x +8 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) roots ( p ) disp ( l e t ) 17
14 x1 = -0.7784571 15 x2 =2.2891685 16 x3 =4.4892886 17 disp ( now , s i n c e we want e q u a t i o n whose sum o f
18 19 20 21 22 23 24 25
r o o t s i s 0 . sum o f r o o t s o f a b o v e e q u a t i o n i s 6 , s o we w i l l d e c r e a s e ) disp ( v a l u e o f e a c h r o o t by 2 i . e . x4=x1 2 ) x4 = x1 -2 disp ( x5=x2 2 ) x5 = x2 -2 disp ( x6=x3 2 ) x6 = x3 -2 disp ( hence , t h e r e q u i r e d e q u a t i o n i s ( xx4 ) ( xx5 ) ( xx6 ) =0 > ) p1 =( x - x4 ) *( x - x5 ) *( x - x6 )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =6*( x ^5) -41*( x ^4) +97*( x ^3) -97*( x ^2) +41* x -6 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =6*( x ^6) -25*( x ^5) +31*( x ^4) -31*( x ^2) +25* x -6 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
18
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^3 -3*( x ^2) +12* x +16 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p =28*( x ^3) -9*( x ^2) +1 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^3+ x ^2 -16* x +20 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
19
Scilab code Exa 1.18 Finding the roots of equation containing one variable
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^3 -3*( x ^2) +3 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^4 -12*( x ^3) +41*( x ^2) -18* x -72 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^4 -2*( x ^3) -5*( x ^2) +10* x -3 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
20
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^4 -8*( x ^2) -24* x +7 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
1 clear 2 clc 3 x = poly ([0] , x ) ; 4 p = x ^4 -6*( x ^3) -3*( x ^2) +22* x -6 5 disp ( t h e r o o t s o f a b o v e e q u a t i o n a r e ) 6 roots ( p )
Scilab code Exa 1.23 Finding the solution of equation by drawing graphs
1 2 3 4 5 6 7 8 9 10
clear clc xset ( window ,1) xtitle ( My Graph , X a x i s , Y a x i s ) x = linspace (1 ,3 ,30) y1 =3 - x y2 = %e ^( x -1) plot (x , y1 , o ) plot (x , y2 , + ) legend ( 3 x , %e ( x 1) ) 21
Figure 1.1: Finding the solution of equation by drawing graphs

11
disp ( from t h e graph , i t i s c l e a r t h a t t h e p o i n t o f i n t e r s e c t i o n i s n e a r l y x =1.43 )
1 2 3 4 5 6 7 8 9 10 11
clear clc xset ( window ,2) xtitle ( My Graph , X a x i s , Y a x i s ) x = linspace (1 ,3 ,30) y1 = x y2 = sin ( x ) + %pi /2 plot (x , y1 , o ) plot (x , y2 , + ) legend ( x , s i n ( x )+%pi /2 ) disp ( from t h e graph , i t i s c l e a r t h a t t h e p o i n t o f i n t e r s e c t i o n i s n e a r l y x =2.3 )
22
1 2 3 4 5 6 7 8 9 10 11
clear clc xset ( window ,3) xtitle ( My Graph , X a x i s , Y a x i s ) x = linspace (0 ,3 ,30) y1 = - sec ( x ) y2 = cosh ( x ) plot (x , y1 , o ) plot (x , y2 , + ) legend ( s e c ( x ) , c o s h ( x ) ) disp ( from t h e graph , i t i s c l e a r t h a t t h e p o i n t o f i n t e r s e c t i o n i s n e a r l y x =2.3 )
23
24
Chapter 2 Determinants and Matrices
Scilab code Exa 2.1 Calculating Determinant

1 clc 2 syms a ; 3 syms h ; 4 syms g ; 5 syms b ; 6 syms f ; 7 syms c ; 8 A =[ a h g ; h b f ; g f c ] 9 det ( A )

1 clear 2 clc 3 a =[0 1 2 3;1 0 3 0;2 3 0 1;3 0 1 2] 4 disp ( d e t e r m i n a n t o f a i s ) 5 det ( a )
25

1 clc 2 syms a ; 3 syms b ; 4 syms c ; 5 A =[ a a ^2 a ^3 -1; b b ^2 b ^3 -1; c c ^2 c ^3 -1] 6 det ( A )

1 clear 2 clc 3 a =[21 17 7 10;24 22 6 10;6 8 2 3;6 7 1 2] 4 disp ( d e t e r m i n a n t o f a i s ) 5 det ( a )
Scilab code Exa 5.8 Partial derivative of given function

1 2 3 4 5 6 7 8 9 10
clc syms x y u=x^y a = diff (u , y ) b = diff (a , x ) c = diff (b , x ) d = diff (u , x ) e = diff (d , y ) f = diff (e , x ) disp ( c l e a r l y , c=f ) 26
Scilab code Exa 2.16 product of two matrices

1 2 3 4 5 6 7 8
clear clc A =[0 1 2;1 2 3;2 3 4] B =[1 -2; -1 0;2 -1] disp ( AB= ) A*B disp ( BA= ) B*A
Scilab code Exa 2.17 Product of two matrices

1 2 3 4 5 6 7 8 9
clear clc A =[1 3 0; -1 2 1;0 0 2] B =[2 3 4;1 2 3; -1 1 2] disp ( AB= ) A*B disp ( BA= ) B*A disp ( c l e a r l y AB i s n o t e q u a l t o BA )
Scilab code Exa 2.18 Product and inverse of matrices

1 clear 2 clc 3 A =[3 2 2;1 3 1;5 3 4] 4 C =[3 4 2;1 6 1;5 6 4]
27
5 disp ( AB=C >B=i n v (A) C ) 6 B = inv ( A ) * C
Scilab code Exa 2.19 Solving equation of matrices

1 2 3 4 5 6
clear clc A =[1 3 2;2 0 -1;1 2 3] I = eye (3 ,3) disp ( A3 4 A2 3A+11 I= ) A ^3 -4* A *A -3* A +11* I
Scilab code Exa 2.20 Nth power of a given matrix

1 2 3 4 5
clc A =[11 -25;4 -9] n = input ( E n t e r t h e v a l u e o f n ) ; d i s p ( calculating A ^ n ) ; A n
Scilab code Exa 2.23 Inverse of matrix

1 clear 2 clc 3 A =[1 1 3;1 3 -3; -2 -4 -4] 4 disp ( i n v e r s e o f A i s ) 5 inv ( A )
28
Scilab code Exa 2.24.1 Rank of a matrix

1 clear 2 clc 3 A =[1 2 3;1 4 2;2 6 5] 4 disp ( Rank o f A i s ) 5 rank ( A )
Scilab code Exa 2.24.2 Rank of a matrix

1 clear 2 clc 3 A =[0 1 -3 -1;1 0 1 1;3 1 0 2;1 1 -2 0] 4 disp ( Rank o f A i s ) 5 rank ( A )
Scilab code Exa 2.25 Inverse of matrix

1 clear 2 clc 3 A =[1 1 3;1 3 -3; -2 -4 -4] 4 disp ( i n v e r s e o f A i s ) 5 inv ( A )
Scilab code Exa 2.26 eigen values vectors rank of matrix

1 clear 2 clc 3 A =[2 3 -1 -1;1 -1 -2 -4;3 1 3 -2;6 3 0 -7] 4 [ R P ]= spec ( A )
29
5 6
disp ( r a n k o f A ) rank ( A )
Scilab code Exa 2.28 Inverse of a matrix

1 clear 2 clc 3 A =[1 1 1;4 3 -1;3 5 3] 4 disp ( i n v e r s e o f A = ) 5 inv ( A )
Scilab code Exa 2.31 Solving equation using matrices

1 clear 2 clc 3 disp ( t h e e q u a t i o n s can be r e w r i t t e n a s AX=B where 4 5 6 7 8 9 10 11
X=[ x1 ; x2 ; x3 ; x4 ] and ) A =[1 -1 1 1;1 1 -1 1;1 1 1 -1;1 1 1 1] B =[2; -4;4;0] disp ( d e t e r m i n a n t o f A= ) det ( A ) disp ( i n v e r s e o f A = ) inv ( A ) disp ( X= ) inv ( A ) * B
Scilab code Exa 2.32 Solving equation using matrices

1 clear 2 clc
30
3 4 5 6 7 8
disp ( t h e e q u a t i o n s can be r e w r i t t e n a s AX=B where X=[ x ; y ; z ] and ) A =[5 3 7;3 26 2;7 2 10] B =[4;9;5] disp ( d e t e r m i n a n t o f A= ) det ( A ) disp ( S i n c e d e t (A) =0 , hence , t h i s s y s t e m o f e q u a t i o n w i l l have i n f i n i t e s o l u t i o n s . . hence , t h e s y s t e m i s c o n s i s t e n t )
Scilab code Exa 2.34.1 predicting nature of equation using rank of matrix
1 2 3 4 5 6
clc A =[1 2 3;3 4 4;7 10 12] disp ( r a n k o f A i s ) p = rank ( A ) if p ==3 then disp ( e q u a t i o n s have o n l y a t r i v i a l s o l u t i o n : x=y=z =0 ) 7 else 8 disp ( e q u a t i o n s have i n f i n i t e no . o f s o l u t i o n s . ) 9 end
Scilab code Exa 2.34.2 predicting nature of equation using rank of matrix
1 2 3 4 5 6
clc A =[4 2 1 3;6 3 4 7;2 1 0 1] disp ( r a n k o f A i s ) p = rank ( A ) if p ==4 then disp ( e q u a t i o n s have o n l y a t r i v i a l =0 ) 7 else 31
s o l u t i o n : x=y=z
8 9
disp ( e q u a t i o n s have i n f i n i t e no . o f s o l u t i o n s . ) end
Scilab code Exa 2.38 Inverse of a matrix

1 clear ; 2 clc ; 3 disp ( t h e g i v e n e q u a t i o n s can be w r i t t e n a s Y=AX 4 5 6 7 8 9
where ) A =[2 1 1;1 1 2;1 0 -2] disp ( d e t e r m i n a n t o f A i s ) det ( A ) disp ( s i n c e , i t s non s i n g u l a r , h e n c e t r a n s f o r m a t i o n i s r e g u l a r ) disp ( i n v e r s e o f A i s ) inv ( A )
Scilab code Exa 2.39 Transpose and product of matrices

1 2 3 4 5 6 7 8
clear clc A =[ -2/3 1/3 2/3;2/3 2/3 1/3;1/3 -2/3 2/3] disp ( A t r a n s p o s e i s e q u a l t o ) A disp ( A ( t r a n s p o s e o f A)= ) A *A disp ( hence , A i s o r t h o g o n a l )
Scilab code Exa 2.42 eigen values and vectors of given matrix
32
1 clear 2 clc 3 A =[5 4;1 2] 4 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f
5 6 7 8 9
transformation and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) [ R P ]= spec ( A ) disp ( R i s n o r m a l i s e d . l e t U r e p r e s e n t s u n n o r m a l i s e d v e r s i o n of r ) U (: ,1) = R (: ,1) * sqrt (17) ; U (: ,2) = R (: ,2) * sqrt (2) disp ( two e i g e n v e c t o r s a r e t h e two c o l u m n s o f U )
1 clear 2 clc 3 A =[1 1 3;1 5 1;3 1 1] 4 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f
5 6 7 8 9 10
transformation and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) [ R P ]= spec ( A ) disp ( R i s n o r m a l i s e d . l e t U r e p r e s e n t s u n n o r m a l i s e d v e r s i o n of r ) U (: ,1) = R (: ,1) * sqrt (2) ; U (: ,2) = R (: ,2) * sqrt (3) ; U (: ,3) = R (: ,3) * sqrt (6) disp ( t h r e e e i g e n v e c t o r s a r e t h e t h r e e c o l u m n s o f U )
1
clear 33
2 clc 3 A =[3 1 4;0 2 6;0 0 5] 4 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f
5 6 7 8 9 10
transformation and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) [ R P ]= spec ( A ) disp ( R i s n o r m a l i s e d . l e t U r e p r e s e n t s u n n o r m a l i s e d v e r s i o n of r ) U (: ,1) = R (: ,1) * sqrt (1) ; U (: ,2) = R (: ,2) * sqrt (2) ; U (: ,3) = R (: ,3) * sqrt (14) disp ( t h r e e e i g e n v e c t o r s a r e t h e t h r e e c o l u m n s o f U )
Scilab code Exa 2.45 eigen values and characteristic equation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
clear clc x = poly ([0] , x ) A =[1 4;2 3] I = eye (2 ,2) disp ( e i g e n v a l u e s o f A a r e ) spec ( A ) disp ( l e t ) a = -1; b =5; disp ( hence , t h e c h a r a c t e r i s t i c e q u a t i o n i s ( xa ) ( xb ) ) p =( x - a ) *( x - b ) disp ( A2 4 A5 I= ) A ^2 -4* A -5* I disp ( i n v e r s e o f A= ) inv ( A )
34

1 2 3 4 5 6 7 8 9 10 11
clear clc x = poly ([0] , x ) A =[1 1 3;1 3 -3; -2 -4 -4] disp ( e i g e n v a l u e s o f A a r e ) spec ( A ) disp ( l e t ) a =4.2568381; b =0.4032794; c = -4.6601175; disp ( hence , t h e c h a r a c t e r i s t i c e q u a t i o n i s ( xa ) ( xb ) ( xc ) ) 12 p =( x - a ) *( x - b ) *( x - c ) 13 disp ( i n v e r s e o f A= ) 14 inv ( A )

1 2 3 4 5 6 7 8 9 10 11
clear clc x = poly ([0] , x ) A =[2 1 1;0 1 0;1 1 2] I = eye (3 ,3) disp ( e i g e n v a l u e s o f A a r e ) spec ( A ) disp ( l e t ) a =1; b =1; c =3;
35
disp ( hence , t h e c h a r a c t e r i s t i c e q u a t i o n i s ( xa ) ( xb ) ( xc ) ) 13 p =( x - a ) *( x - b ) *( x - c ) 14 disp ( A8 5 A7+7 A6 3 A5+A4 5 A3+8 A2 2 A+I = ) 15 A ^8 -5* A ^7+7* A ^6 -3* A ^5+ A ^4 -5* A ^3+8* A ^2 -2* A + I
12
1 clear 2 clc 3 A =[ -1 2 -2;1 2 1; -1 -1 0] 4 disp ( R i s m a t r i x o f t r a n s f o r m a t i o n and D i s a
d i a g o n a l matrix ) 5 [ R D ]= spec ( A )
1 clear 2 clc 3 A =[1 1 3;1 5 1;3 1 1] 4 disp ( R i s m a t r i x o f t r a n s f o r m a t i o n and D i s a
d i a g o n a l matrix )
5 [ R D ]= spec ( A ) 6 disp ( R i s n o r m a l i s e d , l e t P d e n o t e s u n n o r m a l i s e d 7 8 9 10 11
v e r s i o n o f R . Then ) P (: ,1) = R (: ,1) * sqrt (2) ; P (: ,2) = R (: ,2) * sqrt (3) ; P (: ,3) = R (: ,3) * sqrt (6) disp ( A4= ) A ^4
36
1 clear 2 clc 3 disp ( 3 x 2+5 y 2+3 z 2 2 y z +2 z x 2 x y ) 4 disp ( The m a t r i x o f t h e g i v e n q u a d r a t i c form i s ) 5 A =[3 -1 1; -1 5 -1;1 -1 3] 6 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f t r a n s f o r m a t i o n
and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) 7 [ R P ]= spec ( A ) 8 disp ( so , c a n o n i c a l form i s 2 x 2+3 y 2+6 z 2 )
1 clear 2 clc 3 disp ( 2 x1 x2 +2 x1 x3 2 x2 x3 ) 4 disp ( The m a t r i x o f t h e g i v e n q u a d r a t i c form i s ) 5 A =[0 1 1;1 0 -1;1 -1 0] 6 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f t r a n s f o r m a t i o n
and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) 7 [ R P ]= spec ( A ) 8 disp ( so , c a n o n i c a l form i s 2 x2+y2+ z 2 )
Scilab code Exa 2.52 Hermitian matrix

1
clear 37
2 3 4 5 6 7 8
clc A =[2+ %i 3 -1+3* %i ; -5 %i 4 -2* %i ] disp ( A= ) A disp ( AA= ) A *A disp ( c l e a r l y ,AA i s h e r m i t i a n m a t r i x )
Scilab code Exa 2.53 tranpose and inverse of complex matrix

1 clear 2 clc 3 A =[(1/2) *(1+ %i ) (1/2) *( -1+ %i ) ;(1/2) *(1+ %i ) (1/2) *(1 4 5 6 7 8 9 10 11
%i ) ] disp ( A= ) A disp ( AA= ) A *A disp ( AA= ) A * A disp ( i n v e r s e o f A i s ) inv ( A )
Scilab code Exa 2.54 Unitary matrix

1 2 3 4 5 6 7 8
clear clc A =[0 1+2* %i ; -1+2* %i 0] I = eye (2 ,2) disp ( I A= ) I-A disp ( i n v e r s e o f ( I+A)= ) inv ( I + A ) 38
9 disp ( ( ( I A) ( i n v e r s e ( I+A) ) ) ( ( I A) ( i n v e r s e ( I+A) ) )= ) 10 ((( I - A ) *( inv ( I + A ) ) ) ) *(( I - A ) *( inv ( I + A ) ) ) 11 disp ( ( ( I A) ( i n v e r s e ( I+A) ) ) ( ( I A) ( i n v e r s e ( I+A) ) ) = ) 12 (( I - A ) *( inv ( I + A ) ) ) *((( I - A ) *( inv ( I + A ) ) ) ) 13 disp ( c l e a r l y , t h e p r o d u c t i s an i d e n t i t y m a t r i x .
hence , i t
i s a un ita ry matrix )
39
Chapter 4 Dierentiation and Applications
Scilab code Exa 4.4.1 nding nth derivative

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
// q u e s 4 . 1 // c l e a r // cd SCI // cd ( . . ) // cd ( . . ) // e x e c s y m b o l i c . s c e clc disp ( we have t o f i n d yn f o r F=c o s x c o s 2 x c o s 3 x ) ; syms x F = cos ( x ) * cos (2* x ) * cos (3* x ) ; n = input ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; d i s p ( calculating yn ) ; yn= d i f f ( F , x , n ) d i s p ( the expression for yn is ) ; d i s p ( yn ) ;
40
Scilab code Exa 4.5 nding nth derivative

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
// q u e s 4 . 1 // c l e a r // cd SCI // cd ( . . ) // cd ( . . ) // e x e c s y m b o l i c . s c e clc disp ( we have t o f i n d yn f o r F=c o s x c o s 2 x c o s 3 x ) ; syms x F = x /(( x -1) *(2* x +3) ) ; n = input ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n : ) ; d i s p ( calculating yn ) ; yn= d i f f ( F , x , n ) d i s p ( the expression for yn is ) ; d i s p ( yn ) ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
// q u e s 4 . 1 // c l e a r // cd SCI // cd ( . . ) // cd ( . . ) // e x e c s y m b o l i c . s c e clc disp ( we have t o f i n d yn f o r F=c o s x c o s 2 x c o s 3 x ) ; syms x a F = x /( x ^2+ a ^2) ; n = input ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n : ) ; d i s p ( calculating yn ) ; yn= d i f f ( F , x , n ) d i s p ( the expression for yn is ) ; d i s p ( yn ) ; 41

1 // q u e s 4 . 1 2 // c l e a r 3 // cd SCI 4 // cd ( . . ) 5 // cd ( . . ) 6 // e x e c s y m b o l i c . s c e 7 clc 8 disp ( we have t o f i n d yn f o r F=c o s x c o s 2 x c o s 3 x ) ; 9 syms x a 10 F = %e ^( x ) *(2* x +3) ^3; 11 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n : ) ; 12 disp ( c a l c u l a t i n g yn ) ; 13 yn = diff (F ,x , n ) 14 disp ( t h e e x p r e s s i o n f o r yn i s ) ; 15 disp ( yn ) ;
Scilab code Exa 4.8 proving the given dierential equation

1 // q u e s 4 . 1 2 // c l e a r 3 // cd SCI 4 // cd ( . . ) 5 // cd ( . . ) 6 // e x e c s y m b o l i c . s c e 7 clc 8 disp ( y=( s i n 1) x ) s i g n i n v e r s e x ) ; 9 syms x 10 y =( asin ( x ) ) ^2; 11 disp ( we have t o p r o v e (1 x 2 ) y ( n+2) (2n+1) xy ( n+1)n
2 yn ) ; 42
12 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; 13 disp ( c a l c u l a t i n g yn f o r v a r i o u s v a l u e s o f n ) ; 14 for n =1:4 15 16 F =(1 - x ^2) * diff (y ,x , n +2) -(2* n +1) * x * diff (y ,x , n +1) -( n 17 18 19 20 21 22 end 23 disp ( Hence p r o v e d ) ;
^2+ a ^2) * diff (y ,x , n ) ; disp ( n ) ; disp ( t h e e x p r e s s i o n f o r yn i s disp ( F ) ; disp ( Which i s e q u a l t o 0 ) ;
);

1 // q u e s 4 . 1 2 // c l e a r 3 // cd SCI 4 // cd ( . . ) 5 // cd ( . . ) 6 // e x e c s y m b o l i c . s c e 7 clc 8 disp ( y=e ( a ( s i n 1) x ) ) s i g n i n v e r s e x ) ; 9 syms x a 10 y = %e ^( a *( asin ( x ) ) ) ; 11 disp ( we have t o p r o v e (1 x 2 ) y ( n+2) (2n+1) xy ( n+1) (
n2+a 2 ) yn ) ; 12 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; 13 disp ( c a l c u l a t i n g yn f o r v a r i o u s v a l u e s o f n ) ; 14 for n =1:4

15 16 17
// yn= d i f f ( F , x , n ) F =(1 - x ^2) * diff (y ,x , n +2) -(2* n +1) * x * diff (y ,x , n +1) -( n ^2+ a ^2) * diff (y ,x , n ) ; 43
18 disp ( n ) ; 19 disp ( t h e e x p r e s s i o n f o r yn i s 20 disp ( F ) ; 21 disp ( Which i s e q u a l t o 0 ) ; 22 23 end 24 disp ( Hence p r o v e d ) ;
);

1 clc 2 disp ( y ( 1 /m)+y (1/m) =2x ); 3 disp ( OR y ( 2 /m) 2xy ( 1 /m) +1 ) ; 4 disp ( OR y =[ x+(x 2 1) ] m and y =[ x (x 2 1) ] m ) ; 5 6 syms x m 7 disp ( For y =[ x+(x 2 1) ] m ) ; 8 y =( x +( x ^2 -1) ) ^ m 9 disp ( we have t o p r o v e ( x 2 1) y ( n+2) +(2n+1) xy ( n+1)+(
n2 m 2 ) yn ) ; 10 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; 11 disp ( c a l c u l a t i n g yn f o r v a r i o u s v a l u e s o f n ) ; 12 for n =1:4

13 14 15
16 17 18 19 20 21 end 22 disp ( For y =[x (x 2 1) ] m ) ; 23 y =( x -( x ^2 -1) ) ^ m
// yn= d i f f ( F , x , n ) F =( x ^2 -1) * diff (y ,x , n +2) +(2* n +1) * x * diff (y ,x , n +1) +( n ^2 - m ^2) * diff (y ,x , n ) ; disp ( n ) ; disp ( t h e e x p r e s s i o n f o r yn i s ) ; disp ( F ) ; disp ( Which i s e q u a l t o 0 ) ;
44
disp ( we have t o p r o v e ( x 2 1) y ( n+2) +(2n+1) xy ( n+1)+( n2 m 2 ) yn ) ; 25 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; 26 disp ( c a l c u l a t i n g yn f o r v a r i o u s v a l u e s o f n ) ; 27 for n =1:4

24 28 29 30
31 32 33 34 35 36 end 37 disp ( Hence p r o v e d ) ;
// yn= d i f f ( F , x , n ) F =( x ^2 -1) * diff (y ,x , n +2) +(2* n +1) * x * diff (y ,x , n +1) +( n ^2 - m ^2) * diff (y ,x , n ) ; disp ( n ) ; disp ( t h e e x p r e s s i o n f o r yn i s ) ; disp ( F ) ; disp ( Which i s e q u a l t o 0 ) ;
Scilab code Exa 4.11 verify roles theorem

1 clc 2 disp ( f o r
3 4 5 6 7 8 y1 = diff (y , x ) ; 9 disp ( y1 ) ; 10 disp ( p u t t i n g t h i s 11
r o l e s t h e o r e m F9x ) s h o u l d be d i f f e r e n t i a b l e i n ( a , b ) and f ( a )=f ( b ) ) ; disp ( Here f ( x )=s i n ( x ) / e x ) ; disp ( ) ; syms x y = sin ( x ) / %e ^ x ;
t o z e r o we g e t t a n ( x ) =1 i e x=p i /4
); disp ( v a l u e p i /2 l i e s b/w 0 and p i . Hence r o l e s theorem i s v e r i f i e d );
45
Scilab code Exa 4.16 expansion using maclaurins series

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// q u e s 1 6 disp ( M a c l a u r i n s s e r i e s ) ; disp ( f ( x )=f ( 0 )+x f 1 ( 0 )+x 2 / 2 ! f 2 ( 0 )+x 3 / 3 ! f 3 ( 0 ) + . . . . . . ); syms x a // f u n c t i o n y=f ( a ) y = tan ( a ) ; // e n d f u n c t i o n n = input ( e n t e r t h e number o f e x p r e s s i o n i n s e r i e s : ); a =1; t = eval ( y ) ; a =0; for i =2: n y1 = diff (y , a ,i -1) ; t = t + x ^( i -1) * eval ( y1 ) / factorial (i -1) ; end disp ( t )
Scilab code Exa 4.17 expanding function as fourier series of sine term
1 2 3
// q u e s 1 6 disp ( M a c l a u r i n s s e r i e s ) ; disp ( f ( x )=f ( 0 )+x f 1 ( 0 )+x 2 / 2 ! f 2 ( 0 )+x 3 / 3 ! f 3 ( 0 ) + . . . . . . ); 4 syms x a

5 6 7
y = %e ^( sin ( a ) ) ; n = input ( e n t e r t h e number o f e x p r e s s i o n i n s e r i s : ); 46
8 9 10 11 12 13 14 15
a =0; t = eval ( y ) ; a =0; for i =2: n y1 = diff (y , a ,i -1) ; t = t + x ^( i -1) * eval ( y1 ) / factorial (i -1) ; end disp ( t )

1 2 3
// q u e s 1 8 disp ( M a c l a u r i n s s e r i e s ) ; disp ( f ( x )=f ( 0 )+x f 1 ( 0 )+x 2 / 2 ! f 2 ( 0 )+x 3 / 3 ! f 3 ( 0 ) + . . . . . . ); 4 syms x a

5 6 7
8 9 10 11 12 13 14 15
y = log (1+( sin ( a ) ) ^2) ; n = input ( e n t e r t h e number o f d i f f e r e n t i a t i o n involved in m a c l a u r i n s s e r i e s : ); a =0; t = eval ( y ) ; a =0; for i =2: n y1 = diff (y , a ,i -1) ; t = t + x ^( i -1) * eval ( y1 ) / factorial (i -1) ; end disp ( t )

1 2
// q u e s 1 9 disp ( M a c l a u r i n s s e r i e s ) ; 47
disp ( f ( x )=f ( 0 )+x f 1 ( 0 )+x 2 / 2 ! f 2 ( 0 )+x 3 / 3 ! f 3 ( 0 ) + . . . . . . ); 4 syms x a b

5 6 7
y = %e ^( a * asin ( b ) ) ; n = input ( e n t e r t h e number o f e x p r e s s i o n i n s e r i s : ); 8 b =0; 9 t = eval ( y ) ;

10 11 for i =2: n 12 y1 = diff (y , b ,i -1) ; 13 t = t + x ^( i -1) * eval ( y1 ) / factorial (i -1) ; 14 end 15 disp ( t )
Scilab code Exa 4.20 expansion using taylors series

1 2 3 4 5 6
// q u e s 2 0 disp ( Advantage o f s c i l a b i s t h a t we can c a l c u l a t e l o g 1 . 1 d i r e c t l y without u s i n g Taylor s e r i e s ); disp ( Use o f t a y l o r s e r i e s a r e g i v e n i n s u b s e q u e n t examples ); y = log (1.1) ; disp ( l o g ( 1 . 1 ) = ) ; disp ( log (1.1) ) ;
Scilab code Exa 4.21 taylor series

1 2 3
// q u e s 2 1 disp ( T a y l o r s e r i e s ) ; disp ( f ( x+h )=f ( x )+h f 1 ( x )+h 2 / 2 ! f 2 ( x )+h 3 / 3 ! f 3 ( x ) + . . . . . . ); 48
4 5 6 7 8
disp ( To f i n f t h e t a y l o r e x p a n s i o n o f tan 1( x+h ) ) syms x h y = atan ( x ) ; n = input ( e n t e r t h e number o f e x p r e s s i o n i n s e r i s : );
9 10 t = y ; 11 12 for i =2: n 13 y1 = diff (y , x ,i -1) ; 14 t = t + h ^( i -1) *( y1 ) / factorial (i -1) ; 15 end 16 disp ( t )
Scilab code Exa 4.22 evaluating limit

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
// q u e s 2 2 disp ( Here we n e e d t o f i n d f i n d t h e l i m i t o f f ( x ) a t x=0 ) syms x y =( x * %e ^x - log (1+ x ) ) / x ^2; // d i s p ( The l i m i t a t x=0 i s : ) ; // l = l i m i t ( y , x , 0 ) ; // d i s p ( l ) f =1; while f ==1 yn = x * %e ^x - log (1+ x ) ; yd = x ^2; yn1 = diff ( yn , x ,1) ; yd1 = diff ( yd , x ,1) ; x =0; a = eval ( yn1 ) ; b = eval ( yd1 ) ; if a == b then 49
18 yn = yn1 ; 19 yd = yd1 ; 20 else 21 f =0; 22 23 end 24 end 25 h = a / b ; 26 disp ( h ) ;
Scilab code Exa 4.32 tangent to curve

1 2 3 4 5 6 7 8 9 10 11
// q u e s 32 disp ( E q u a t i o n o f t a n g e n t ) ; syms x a y ; f =( a ^(2/3) -x ^(2/3) ) ^(3/2) ; s = diff (f , x ) ; Y1 = s *( - x ) + y ; X1 = - y / s * x ; g =x -( Y1 - s *( X1 - x ) ) ; disp ( E q u a t i o n i s g=0 where g i s ) ; disp ( g ) ;
Scilab code Exa 4.34 nding equation of normal

1 2 3 4 5 6 7
// q u e s 3 4 disp ( E q u a t i o n o f t a n g e n t ) ; syms x a t y xo = a *( cos ( t ) + t * sin ( t ) ) ; yo = a *( sin ( t ) -t * cos ( t ) ) ; s = diff ( xo , t ) / diff ( yo , t ) ; y = yo + s *( x - xo ) ; 50
8 9
disp ( y= ) ; disp ( y ) ;
Scilab code Exa 4.35 nding angle of intersection of curve

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// q u e s 3 5 disp ( The two g i v e n c u r v e s a r e x=4y and y2=4x which i n t e r s e c t s a t ( 0 , 0 ) and ( 4 , 4 ) ) ; disp ( for (4 ,4) ) ; x =4; syms x y1=x 2 / 4 ; y2 =2 x ( 1 / 2 ) ; m1= d i f f ( y1 , x , 1 ) ; m2= d i f f ( y2 , x , 1 ) ; x =4; m1=e v a l (m1) ; m2=e v a l (m2) ; d i s p ( A n g l e b e t w e e n them i s ( r a d i a n s ) : ) ; t=a t a n ( ( m1m2) /(1+m1 m2) ) ; disp ( t ) ;
Scilab code Exa 4.37 prove given tangent statement

1 2 3 4 5 6 7 8
// q u e s 3 7 syms a t x = a *( cos ( t ) + log ( tan ( t /2) ) ) ; y = a * sin ( t ) ; s = diff (x ,t ,1) / diff (y ,t ,1) ; disp ( l e n g t h o f t a n g e n t ) ; l = y *(1+ s ) ^(0.5) ; disp ( l ) ; 51
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
disp ( c h e c k i n g f o r i t s d e p e n d e n c y on t ) f =1 t =0; k = eval ( l ) ; for i =1:10 t=i; if ( eval ( l ) ~= k ) f =0; end end if ( f ==1) disp ( v e r i f i e d and e q u a l t o a ) ; disp ( s u b t a n g e n t ) ; m=y/s; disp ( m ) ;
Scilab code Exa 4.39 nding angle of intersection of curve

1 // q u e s 3 9 2 clc 3 disp ( A n g l e o f 4 disp ( p o i n t o f
i n t e r s e c t i o n ); i n t e r s e c t i o n o f r= s i n t +c o s t and r =2 s i n t i s t=p i /4 ) ; 5 disp ( t a n u=dQ/ d r r ) ; 6 syms Q ;
r1 =2* sin ( Q ) ; r2 = sin ( Q ) + cos ( Q ) ; u = atan ( r1 * diff ( r2 ,Q ,1) ) ; Q = %pi /4; u = eval ( u ) ; disp ( The a n g l e a t p o i n t o f i n t e r s e c t i o n i n r a d i a n s i s : ); 14 disp ( u ) ; 52
7 8 9 10 11 12 13
Scilab code Exa 4.41 nding pedal equation of parabola

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// q u e s 4 1 clc disp ( t a n u=dQ/ d r r ) ; syms Q a ; r =2* a /(1 - cos ( Q ) ) ; u = atan ( r / diff ( r2 ,Q ,1) ) ; u = eval ( u ) ; p = r * sin ( u ) ; syms r ; Q = acos (1 -2* a / r ) ; // c o s (Q) =12 a / r ; p = eval ( p ) ; disp ( p ) ;
Scilab code Exa 4.43 nding radius of curvature of cycloid

1 2 3 4 5 6 7 8 9 10
// q u e s 4 3 syms a t x = a *( t + sin ( t ) ) ; y = a *(1 - cos ( t ) ) ; s2 = diff (y ,t ,2) / diff (x ,t ,2) ; s1 = diff (y ,t ,1) / diff (x ,t ,1) ; r =(1+ s1 ^2) ^(3/2) / s2 ; disp ( The r a d i u s o f c u r v a t u r e i s : disp ( r ) ; 53
);
Scilab code Exa 4.46 radius of curvature of cardoid

1 2 3 4 5 6 7 8 9 10 11
// q u e s 4 6 disp ( r a d i u s o f c u r v a t u r e ) ; syms a t r = a *(1 - cos ( t ) ) ; r1 = diff (r ,t ,1) ; l =( r ^2+ r1 ^2) ^(3/2) /( r ^2+2* r1 ^2 - r * r1 ) ; syms r ; t = acos (1 - r / a ) ; l = eval ( l ) ; disp ( l ) ; disp ( Which i s p r o p o r t i o n a l t o r 0 . 5 ) ;
Scilab code Exa 4.47 cordinates of centre of curvature

1 2 3 4 5 6 7 8 9 10 11 12
// q u s 4 7 disp ( The c e n t r e o f c u r v a t u r e ) ; syms x a y y =2*( a * x ) ^0.5; y1 = diff (y ,x ,1) ; y2 = diff (y ,x ,2) ; xx =x - y1 *(1+ y1 ) ^2/ y2 ; yy = y +(1+ y1 ^2) / y2 ; disp ( t h e c o o r d i n a t e s x , y a r e r e s p : ) ; disp ( xx ) ; disp ( yy ) ;
54
Scilab code Exa 4.48 proof statement cycloid

1 2 3 4 5 6 7 8 9 10 11 12 13 14
// q u e s 4 8 disp ( c e n t r e o f c u r v a t u r e o f g i v e n c y c l o i d syms a t x = a *( t - sin ( t ) ) ; y = a *(1 - cos ( t ) ) ; y1 = diff (y ,t ,1) ; y2 = diff (y ,t ,2) ; xx =x - y1 *(1+ y1 ) ^2/ y2 ; yy = y +(1+ y1 ^2) / y2 ;
);
disp ( t h e c o o r d i n a t e s x , y a r e r e s p : ) ; disp ( xx ) ; disp ( yy ) ; disp ( which a n o t h e r p a r a m e t r i c e q u a t i o n o f c y c l o i d );
Scilab code Exa 4.52 maxima and minima

1 2 3 4 5 6 7 8 9
// e r r o r // q u e s 5 2 disp ( To f i n d t h e maxima and minima o f g i v e n f u n c t i o n put f 1 ( x ) =0 ) ; syms x // x=p o l y ( 0 , x ) ; f =3* x ^4 -2* x ^3 -6* x ^2+6* x +1; k = diff (f , x ) ; x = poly (0 , x ) ; k = eval ( k ) ;
Scilab code Exa 4.61 nding the asymptotes of curve 55
1 // q u e s 61 2 clc 3 disp ( t o f i n d t h e a s s y m p t o t e o f g i v e n c u r v e ) ; 4 syms x y 5 f = x ^2* y ^2 - x ^2* y - x * y ^2+ x + y +1; 6 // a=d e g r e e s ( f , x ) ; 7 f1 = coeffs (f ,x ,2) ; 8 disp ( a s s y m p t o t e s p a r a l l e l t o x x i s i s g i v e n by f 1 =0
where f 1 i s : ) ;
9 disp ( factor ( f1 ) ) ; 10 f2 = coeffs (f ,y ,2) ; 11 disp ( a s s y m p t o t e s p a r a l l e l
t o y a x i s i s g i v e n by f 2
=0 and f 2 i s : ) ; 12 disp ( factor ( f2 ) ) ;
56
Chapter 5 Partial Dierentiation And Its Applications

1 2 3 4 5 6 7
clc syms x y z v =( x ^2+ y ^2+ z ^2) ^( -1/2) a = diff (v ,x ,2) b = diff (v ,y ,2) c = diff (v ,z ,2) a+b+c

1 2 3 4 5 6
clc syms x y u = asin (( x + y ) /( x ^0.5+ y ^0.5) ) a = diff (u , x ) b = diff (u , y ) c = diff (a , x ) 57
7 8 9 10 11 12
d = diff (b , y ) e = diff (b , x ) x*a+y*b (1/2) * tan ( u ) ( x ^2) * c +2* x * y * e +( y ^2) * d ( - sin ( u ) * cos (2* u ) ) /(4*( cos ( u ) ) ^3)
Scilab code Exa 5.25.1 Partial derivative of given function

1 2 3 4 5 6 7 8 9 10
clc syms r l x = r * cos ( l ) y = r * sin ( l ) a = diff (x , r ) b = diff (x , l ) c = diff (y , r ) d = diff (y , l ) A =[ a b ; c d ] det ( A )

1 2 3 4 5 6 7 8 9 10 11
clc syms r l z x = r * cos ( l ) y = r * sin ( l ) m=z a = diff (x , r ) b = diff (x , l ) c = diff (x , z ) d = diff (y , r ) e = diff (y , l ) f = diff (y , z ) 58
12 13 14 15 16
g = diff (m , r ) h = diff (m , l ) i = diff (m , z ) A =[ a b c ; d e f ; g h i ] det ( A )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
clc syms r l m x = r * cos ( l ) * sin ( m ) y = r * sin ( l ) * sin ( m ) z = r * cos ( m ) a = diff (x , r ) b = diff (x , m ) c = diff (x , l ) d = diff (y , r ) e = diff (y , m ) f = diff (y , l ) g = diff (z , r ) h = diff (z , m ) i = diff (z , l ) A =[ a b c ; d e f ; g h i ] det ( A )

1 2 3 4 5 6
clc syms x1 x2 x3 y1 =( x2 * x3 ) / x1 y2 =( x3 * x1 ) / x2 y3 =( x1 * x2 ) / x3 a = diff ( y1 , x1 ) 59
7 8 9 10 11 12 13 14 15 16
b = diff ( y1 , x2 ) c = diff ( y1 , x3 ) d = diff ( y2 , x1 ) e = diff ( y2 , x2 ) f = diff ( y2 , x3 ) g = diff ( y3 , x1 ) h = diff ( y3 , x2 ) i = diff ( y3 , x3 ) A =[ a b c ; d e f ; g h i ] det ( A )

1 2 3 4 5 6 7 8 9 10
clc syms x y u = x *(1 - y ^2) ^0.5+ y *(1 - x ^2) ^0.5 v = asin ( x ) + asin ( y ) a = diff (u , x ) b = diff (u , y ) c = diff (v , x ) d = diff (v , y ) A =[ a b ; c d ] det ( A )
60
Chapter 6 Integration and its Applications
Scilab code Exa 6.1.1 indenite integral

1 // q u e s 1 2 disp ( I n d e f i n i t e i n t e g r a l ) ; 3 syms x 4 f = integ (( sin ( x ) ) ^4 , x ) ; 5 disp ( f ) ;
Scilab code Exa 6.1.2 indenite integral

1 // q u e s 1 2 disp ( I n d e f i n i t e i n t e g r a l ) ; 3 syms x 4 f = integ (( cos ( x ) ) ^7 , x ) ; 5 disp ( f ) ;
Scilab code Exa 6.2.1 denite integral
61
1 // q u e s 1 2 disp ( d e f i n i t e i n t e g r a l ) ; 3 syms x 4 f = integ (( cos ( x ) ) ^6 ,x ,0 , %pi /2) ; 5 disp ( float ( f ) ) ;
Scilab code Exa 6.2.2 Denite Integration of a function

1 // no o u t p u t 2 // q u e s 1 3 clc 4 disp ( d e f i n i t e i n t e g r a l ) ; 5 syms x a 6 g = x ^7/( a ^2 - x ^2) ^1/2 7 f = integ (g ,x ,0 , a ) ; 8 disp ( float ( f ) ) ;

1 // e r r o r no o u t p u t 2 // q u e s 4 3 clc 4 disp ( d e f i n i t e i n t e g r a l ) ; 5 syms x a 6 g = x ^3*(2* a *x - x ^2) ^(1/2) ; 7 f = integ (g ,x ,0 ,2* a ) ; 8 disp ( f ) ;
62
1 // no o u t p u t 2 // q u e s 1 3 clc 4 disp ( d e f i n i t e i n t e g r a l ) ; 5 syms x a n 6 g =1/( a ^2+ x ^2) ^ n ; 7 f = integ (g ,x ,0 , %inf ) ; 8 disp ( f ) ;

1 // q u e s 4 2 clc 3 disp ( d e f i n i t e i n t e g r a l ) ; 4 syms x 5 g =( sin (6* x ) ) ^3*( cos (3* x ) ) ^7; 6 f = integ (g ,x ,0 , %pi /6) ; 7 disp ( float ( f ) ) ;

1 // q u e s 4 2 clc 3 disp ( d e f i n i t e i n t e g r a l ) ; 4 syms x 5 g = x ^4*(1 - x ^2) ^(3/2) ; 6 f = integ (g ,x ,0 ,1) ; 7 disp ( float ( f ) ) ;
Scilab code Exa 6.5 denite integral 63
1 2 3 4 5 6 7 8 9 10 11 12 13 14
// e r r o r no i n t e r n a l e r r o r // q u e s 5 clc disp ( d e f i n i t e i n t e g r a l ) ; syms x m n n = input ( E n t e r n : ) ; m = input ( E n t e r m : ) ; g =( cos ( x ) ) ^ m * cos ( n * x ) ; f = integ (g ,x ,0 , %pi /2) ; disp ( float ( f ) ) ; g2 =( cos ( x ) ) ^( m -1) * cos (( n -1) * x ) ; f2 = m /( m + n ) * integ ( g2 ,x ,0 , %pi /2) ; disp ( float ( f2 ) ) ; disp ( Equal ) ;
Scilab code Exa 6.6.1 reducing indenite integral to simpler form

1 2 3 4 5 6 7 8 9
// q u e s 6 clc disp ( d e f i n i t e i n t e g r a l ) ; syms x a n = input ( E n t e r n : ) ; g = exp ( a * x ) *( sin ( x ) ) ^ n ; f = integ (g , x ) ; disp ( f ) ;
Scilab code Exa 6.7.1 Indenite Integration of a function

1 clc 2 syms x 3 disp ( integ ( tan ( x ) ^5 , x ) )
64
Scilab code Exa 6.8 Getting the manual input of a variable and integration
1 2 3 4 5 6
clc n = input ( E n t e r t h e v a l u e o f n ) ; p= i n t e g r a t e ( ( tan ( x ) ) ^( n -1) , x ,0 , %pi /4) q = integrate ( ( t a n ( x ) ) ( n+1) , x ,0 , %pi /4) disp ( n ( p+q )= ) disp ( n *( p + q ) )

1 clear 2 clc 3 integrate ( s e c ( x ) 4 , x ,0 , %pi /4)

1 clear 2 clc 3 integrate ( 1/ s i n ( x ) 3 , x , %pi /3 , %pi /2)
Scilab code Exa 6.10 denite integral

1 2 // q u e s 8 3 clc
65
4 syms x 5 g = x * sin ( x ) ^6* cos ( x ) ^4; 6 f = integ (g ,x ,0 , %pi ) ; 7 disp ( float ( f ) ) ;
Scilab code Exa 6.12 Denite Integration of a function

1 clear 2 clc 3 integrate ( s i n ( x ) 0 . 5 / ( s i n ( x ) 0.5+ c o s ( x ) 0 . 5 ) , x
,0 , %pi /2)
Scilab code Exa 6.13 sum of innite series

1 2 // q u e s 1 3 3 clc 4 syms x 5 disp ( The summation i s
e q u i v a l e n t to i n t e g r a t i o n of 1/(1+ x 2 ) from 0 t o 1 ) ; 6 g =1/(1+ x ^2) ; 7 f = integ (g ,x ,0 ,1) ; 8 disp ( float ( f ) ) ;
Scilab code Exa 6.14 nding the limit of the function

1 // q u e s 1 4 2 clc 3 syms x
66
disp ( The summation i s e q u i v a l e n t t o i n t e g r a t i o n o f l o g (1+ x ) from 0 t o 1 ) ; 5 g = log (1+ x ) ; 6 f = integ (g ,x ,0 ,1) ; 7 disp ( float ( f ) ) ;
4

1 clear 2 clc 3 integrate ( x s i n ( x ) 8 c o s ( x ) 4 , x ,0 , %pi )

1 clear 2 clc 3 integrate ( l o g ( s i n ( x ) ) , x ,0 , %pi /2)
Scilab code Exa 6.24 Calculating the area under two curves
1 2 3 4 5 6 7 8 9 10
clear clc xset ( window ,1) xtitle ( My Graph , X a x i s , Y a x i s ) x = linspace ( -5 ,10 ,70) y1 =( x +8) /2 y2 = x ^2/8 plot (x , y1 , o ) plot (x , y2 , + ) legend ( ( x +8) /2 , x 2 / 8 ) 67
Figure 6.1: Calculating the area under two curves

11
disp ( from t h e graph , i t i s c l e a r t h a t t h e p o i n t s o f i n t e r s e c t i o n a r e x=4 and x =8. ) 12 disp ( So , o u r r e g i o n o f i n t e g r a t i o n i s from x=4 t o x =8 ) 13 integrate ( ( x +8)/2 x 2 / 8 , x , -4 ,8)
68
Chapter 9 Innite Series
Scilab code Exa 9.1 to nd the limit at innity

1 clc 2 syms n ; 3 f =((1/ n ) ^2 -2*(1/ n ) ) /(3*(1/ n ) ^2+(1/ n ) ) 4 disp ( limit (f ,n ,0) ) ;
Scilab code Exa 9.1.3 to nd the limit at innity

1 clc 2 syms n ; 3 f =3+( -1) ^ n 4 limit (f ,n , %inf )
Scilab code Exa 9.2.1 to nd the sum of series upto innity

1 clc 2 syms n
69
3 disp ( 1+2+3+4+5+6+7+....+n + . . . . . = ) 4 p =1/ n *(1/ n +1) /2 5 disp ( limit (p ,n ,0) ) ;
Scilab code Exa 9.2.2 to check for the type of series

1 clc 2 disp ( 5 4 1+5 4 1+5 4 1+5 4 1+.........=0 ,5 ,1
a c c o r d i n g t o t h e no . o f t e r m s . ) 3 disp ( c l e a r l y , i n t h i s c a s e sum d o e s n t t e n d t o a u n i q u e l i m i t . hence , s e r i e s i s o s c i l l a t o r y . )
Scilab code Exa 9.5.1 to check the type of innite series

1 clc 2 syms n ; 3 v =1/((1/ n ) ^2) 4 u =(2/ n -1) /(1/ n *(1/ n +1) *(1/ n +2) ) 5 disp ( limit ( u /v ,n ,0) ) ; 6 disp ( b o t h u and v c o n v e r g e and d i v e r g e
together ,
hence u i s convergent )

1 clc 2 syms n ; 3 v=n 4 u =((1/ n ) ^2) /((3/ n +1) *(3/ n +4) *(3/ n +7) ) 5 disp ( limit ( u /v ,n ,0) ) ; 6 disp ( b o t h u and v c o n v e r g e and d i v e r g e
together ,
hence u i s d i v e r g e n t ) 70

1 clc 2 syms n 3 disp ( u =(( n+1) 0 . 5 1 ) / ( ( n+2) 3 1)=> ) 4 // put n=1/n 5 u =((1+1/(1/ n ) ) -(1/ n ) ^( -0.5) ) /(((1/ n ) ^5/2) *((1+2/(1/ n 6 7 8 9
) ) ^3 -(1/ n ) ^( -3) ) ) v =(1/ n ) ^( -5/2) disp ( limit ( u /v ,n ,0) ) ; // d i s p ( = 1 ) disp ( s i n c e , v i s c o n v e r g e n t , s o u i s a l s o conzavergent . )

1 clc 2 syms n 3 disp ( integ (1/( n * log ( n ) ) ,n ,2 , %inf ) ) ;
Scilab code Exa 9.8.1 to nd the sum of series upto innity

1 clc 2 syms x n ; 3 // put n=1/n 4 u =( x ^(2*(1/ n ) -2) ) /(((1/ n ) +1) *(1/ n ) ^0.5) 5 v =( x ^(2*(1/ n ) ) ) /((1/ n +2) *(1/ n +1) ^0.5) 6 disp ( limit ( u /v ,n ,0) ) ;
71

1 clc 2 syms x n ; 3 // put n=1/n 4 u =((2^(1/ n ) -2) *( x ^(1/ n -1) ) ) /(2^(1/ n ) +1) 5 v =((2^((1/ n ) +1) -2) *( x ^(1/ n ) ) ) /(2^(1/ n +1) +1) 6 disp ( limit ( u /v ,n ,0) ) ;

1 clc 2 syms x n ; 3 u =1/(1+ x ^( - n ) ) ; 4 v =1/(1+ x ^( -n -1) ) ; 5 disp ( limit ( u /v ,n ,0) ) ;

1 clc 2 syms a b n ; 3 l =( b +1/ n ) /( a +1/ n ) 4 disp ( limit (l ,n ,0) ) ;
72
1 clc 2 syms x n ; 3 disp ( u = ( ( 4 . 7 . . . . ( 3 n+1) ) x n ) / ( 1 . 2 . . . . . n ) ) 4 disp ( v = ( ( 4 . 7 . . . . ( 3 n+4) x ( n+1) ) / ( 1 . 2 . . . . . ( n+1) ) ) 5 disp ( l =u / v=> ) 6 l =(1+ n ) /((3+4* n ) * x ) 7 disp ( limit (l ,n ,0) )

1 clc 2 syms x n ; 3 u =((( factorial ( n ) ) ^2) * x ^(2* n ) ) / factorial (2* n ) 4 v =((( factorial ( n +1) ) ^2) * x ^(2*( n +1) ) ) / factorial (2*( n
+1) ) 5 limit ( u /v ,n , %inf )
73
Chapter 10 Fourier Series
Scilab code Exa 10.1 nding fourier series of given function

1 // q u e s 1 2 clc 3 disp ( f i n d i n g t h e 4 5 6 7 8 9 10 11 12 13
fourier
s e r i e s of given function )
; syms x ao =1/ %pi * integ ( exp ( -1* x ) ,x ,0 ,2* %pi ) ; s = ao /2; n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f s i n o r c o s terms i n the ex p an si o n : ); for i =1: n ai =1/ %pi * integ ( exp ( - x ) * cos ( i * x ) ,x ,0 ,2* %pi ) ; bi =1/ %pi * integ ( exp ( - x ) * sin ( i * x ) ,x ,0 ,2* %pi ) ; s = s + float ( ai ) * cos ( i * x ) + float ( bi ) * sin ( i * x ) ; end disp ( float ( s ) ) ;
74
1 2 3 4 5 6 7 8
// e r r o r // q u e s 2 disp ( To f i n d t h e f o u r i e r t r a n s f o r m o f g i v e n f u n c t i o n ); syms x s F = integ ( exp ( %i * s * x ) ,x , -1 ,1) ; disp ( F ) ; // p r o d u c e s e r r o r > F1 = integ ( sin ( x ) /x ,x ,0 , %inf ) ;

1 // q u e s 3 2 clc 3 disp ( f i n d i n g t h e
fourier
;
4 syms x 5 ao =1/ %pi *( integ ( -1* %pi * x ^0 ,x , - %pi ,0) + integ (x ,x ,0 , %pi
));
6 s = ao /2; 7 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f 8 9 10 11 12 13
s i n or c o s terms i n the ex p an si o n : ); for i =1: n ai =1/ %pi *( integ ( -1* %pi * cos ( i * x ) ,x , -1* %pi ,0) + integ ( x * cos ( i * x ) ,x ,0 , %pi ) ) ; bi =1/ %pi *( integ ( -1* %pi * x ^0* sin ( i * x ) ,x , -1* %pi ,0) + integ ( x * sin ( i * x ) ,x ,0 , %pi ) ) ; s = s + float ( ai ) * cos ( i * x ) + float ( bi ) * sin ( i * x ) ; end disp ( float ( s ) ) ;
Scilab code Exa 10.4 nding fourier series of given function 75
fourier
;
4 syms x l 5 ao =1/ l * integ ( exp ( -1* x ) ,x , -l , l ) ; 6 s = ao /2 7 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f
s i n or
cos terms in the expansion :
);
8 for i =1: n 9 ai =1/ l * integ ( exp ( - x ) * cos ( i * %pi * x / l ) ,x , -l , l ) ; 10 bi =1/ l * integ ( exp ( - x ) * sin ( i * %pi * x / l ) ,x , -l , l ) ; 11 s = s + float ( ai ) * cos ( i * %pi * x / l ) + float ( bi ) * sin ( i * %pi * x
/l);
12 end 13 disp ( float ( s ) ) ;
Scilab code Exa 10.5 nding fourier series of given function in interval minus pi to pi
fourier
;
4 syms x l 5 s =0; 6 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f
sin
terms in the expansion : 7 for i =1: n
);
8 9 bi =2/ %pi * integ ( x * sin ( i * x ) ,x ,0 , %pi ) ; 10 s = s + float ( bi ) * sin ( i * x ) ; 11 end 12 disp ( float ( s ) ) ;
76
Scilab code Exa 10.6 nding fourier series of given function in interval minus l to l
1 // e r r o r no o u t p u t 2 // q u e s 6 3 clc 4 disp ( f i n d i n g t h e 5 6 7 8 9 10 11 12 13 14
f o u r i e r s e r i e s of given function ) ; syms x l ao =2/ l * integ ( x ^2 ,x ,0 , l ) ; s = float ( ao ) /2; n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f s i n o r c o s terms i n the ex p an si o n : ); for i =1: n ai =2/ l * integ ( x ^2* cos ( i * %pi * x / l ) ,x ,0 , l ) ; // b i =1/ l i n t e g ( exp ( x ) s i n ( i x ) , x , l , l ) ; s = s + float ( ai ) * cos ( i * %pi * x / l ) ; end disp ( float ( s ) ) ;
fourier
;
4 syms x 5 ao =2/ %pi *( integ ( cos ( x ) ,x ,0 , %pi /2) + integ ( - cos ( x ) ,x ,
%pi /2 , %pi ) ) ; 6 s = ao /2;
77
7 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f
s i n or
);
8 for i =1: n 9 ai =2/ %pi *( integ ( cos ( x ) * cos ( i * x ) ,x ,0 , %pi /2) + integ ( 10 11 12 13
cos ( x ) * cos ( i * x ) ,x , %pi /2 , %pi ) ) ; // b i =1/%pi ( i n t e g ( 1 %pi x 0 s i n ( i x ) , x , 1 %pi , 0 ) + i n t e g ( x s i n ( i x ) , x , 0 , %pi ) ) ; s = s + float ( ai ) * cos ( i * x ) ; end disp ( float ( s ) ) ;
1 // q u e s 8 2 clc 3 disp ( f i n d i n g t h e 4 5 6 7 8 9 10 11 12 13
fourier
; syms x ao =2/ %pi *( integ ((1 -2* x / %pi ) ,x ,0 , %pi ) ) ; s = ao /2; n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f s i n o r c o s terms i n the ex p an si o n : ); for i =1: n ai =2/ %pi *( integ ((1 -2* x / %pi ) * cos ( i * x ) ,x ,0 , %pi ) ) ; // b i =1/%pi ( i n t e g ( 1 %pi x 0 s i n ( i x ) , x , 1 %pi , 0 ) + i n t e g ( x s i n ( i x ) , x , 0 , %pi ) ) ; s = s + float ( ai ) * cos ( i * x ) ; end disp ( float ( s ) ) ;
Scilab code Exa 10.9 nding half range sine series of given function 78
fourier
;
4 syms x l 5 6 s =0; 7 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f
s i n or
);
8 for i =1: n 9 // a i =1/ l i n t e g ( exp ( x ) c o s ( i %pi x / l ) , x , l , l ) ; 10 bi = integ ( x * sin ( i * %pi * x /2) ,x ,0 ,2) ; 11 s = s + float ( bi ) * sin ( i * %pi * x /2) ; 12 end 13 disp ( float ( s ) ) ;
Scilab code Exa 10.10 nding half range cosine series of given function
1 // q u e s 1 0 2 clc 3 disp ( f i n d i n g t h e 4 5 6 7 8 9 10 11 12 13
f o u r i e r s e r i e s of given function ) ; syms x ao =2/2*( integ (x ,x ,0 ,2) ) ; s = ao /2; n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f s i n o r c o s terms i n the ex p an si o n : ); for i =1: n ai =2/2*( integ ( x * cos ( i * %pi * x /2) ,x ,0 ,2) ) ; // b i =1/%pi ( i n t e g ( 1 %pi x 0 s i n ( i x ) , x , 1 %pi , 0 ) + i n t e g ( x s i n ( i x ) , x , 0 , %pi ) ) ; s = s + float ( ai ) * cos ( i * %pi * x /2) ; end disp ( float ( s ) ) ;
79
Scilab code Exa 10.11 expanding function as fourier series of sine term
1 // q u e s 3 2 clc 3 disp ( f i n d i n g t h e 4 5 6 7 8 9 10 11 12
f o u r i e r s e r i e s of given function ) ; syms x ao =0; s = ao ; n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f s i n o r c o s terms i n the ex p an si o n : ); for i =1: n bi =2/1*( integ ((1/4 - x ) * sin ( i * %pi * x ) ,x ,0 ,1/2) + integ (( x -3/4) * sin ( i * %pi * x ) ,x ,1/2 ,1) ) ; s = s + float ( bi ) * sin ( i * %pi * x ) ; end disp ( float ( s ) ) ;

fourier
;
4 syms x 5 ao =1/ %pi * integ ( x ^2 ,x , - %pi , %pi ) ; 6 s = ao /2; 7 n = input ( e n t e r t h e no o f t e r m s u p t o e a c h o f
s i n or
);
8 for i =1: n 9 ai =1/ %pi * integ (( x ^2) * cos ( i * x ) ,x , - %pi , %pi ) ; 10 bi =1/ %pi * integ (( x ^2) * sin ( i * x ) ,x , - %pi , %pi ) ;
80
11 s = s + float ( ai ) * cos ( i * x ) + float ( bi ) * sin ( i * x ) ; 12 end 13 disp ( float ( s ) ) ;
Scilab code Exa 10.13 nding complex form of fourier series

1 // q u e s 1 3 2 clc 3 disp ( The c o m p l e x form o f
s e r i e s i s summation o f f ( n , x ) where n v a r i e s from % i n f t o % i n f and f ( n , x ) i s g i v e n by : ) ; 4 syms n x 5 cn =1/2* integ ( exp ( - x ) * exp ( - %i * %pi * n * x ) ,x , -1 ,1) ; 6 fnx = float ( cn ) * exp ( %i * n * %pi * x ) ;
7 8
disp ( float ( fnx ) ) ;
Scilab code Exa 10.14 practical harmonic analysis

1 2 3 4 5 6 7 8 9 10 11 12
// q u e s 1 5 // yo = [ 1 . 8 0 1 . 1 0 0 . 3 0 0 . 1 6 1 . 5 0 1 . 3 0 2 . 1 6 1 . 2 5 1 . 3 0 1.52 1.76 2.00] // x0 =[0 %pi /6 %pi /3 %pi /2 2 %pi /3 5 %pi /6 %pi 7 %pi /6 4 %pi /3 3 %pi /2 5 %pi /3 11 %pi / 6 ] disp ( P r a c t i c a l h a r m o n i c a n a l y s i s ) ; syms x xo = input ( I n p u t xo m a t r i x : ) ; yo = input ( I n p u t yo m a t r i x : ) ; ao =2* sum ( yo ) / length ( xo ) ; s = ao /2; n = input ( No o f s i n o r c o s term i n e x p a n s i o n : ) ; for i =1: n an =2* sum ( yo .* cos ( i * xo ) ) / length ( yo ) ; 81
13 14 15 16 17
bn =2* sum ( yo .* sin ( i * xo ) ) / length ( yo ) ; s = s + float ( an ) * cos ( i * x ) + float ( bn ) * sin ( i * x ) ; end disp ( s ) ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
// e r r o r // q u e s 1 5 , 1 6 , 1 7 // yo = [ 1 . 9 8 1 . 3 0 1 . 0 5 1 . 3 0 0.88 .25 1 . 9 8 ] // x0 =[0 1/6 1/3 1/2 2/3 5/6 1 ] disp ( P r a c t i c a l h a r m o n i c a n a l y s i s ) ; syms x T xo = input ( I n p u t xo m a t r i x ( i n f a c t o r o f T) : ) ; yo = input ( I n p u t yo m a t r i x : ) ; ao =2* sum ( yo ) / length ( xo ) ; s = ao /2; n = input ( No o f s i n o r c o s term i n e x p a n s i o n : ) ; i =1 an =2*( yo .* cos ( i * xo *2* %pi ) ) / length ( yo ) ; bn =2*( yo .* sin ( i * xo *2* %pi ) ) / length ( yo ) ; s = s + float ( an ) * cos ( i * x *2* %pi / T ) + float ( bn ) * sin ( i * x *2* %pi / T ) ; disp ( s ) ; disp ( D i r e c t c u r r e n t : ) ; i = sqrt ( an ^2+ bn ^2) ;

1 2
// e r r o r // q u e s 1 5 , 1 6 , 1 7 82
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
// yo = [ 1 . 9 8 1 . 3 0 1 . 0 5 1 . 3 0 0.88 .25 1 . 9 8 ] // x0 =[0 1/6 1/3 1/2 2/3 5/6 1 ] disp ( P r a c t i c a l h a r m o n i c a n a l y s i s ) ; syms x T xo = input ( I n p u t xo m a t r i x ( i n f a c t o r o f T) : ) ; yo = input ( I n p u t yo m a t r i x : ) ; ao =2* sum ( yo ) / length ( xo ) ; s = ao /2; n = input ( No o f s i n o r c o s term i n e x p a n s i o n : ) ; i =1 an =2*( yo .* cos ( i * xo *2* %pi ) ) / length ( yo ) ; bn =2*( yo .* sin ( i * xo *2* %pi ) ) / length ( yo ) ; s = s + float ( an ) * cos ( i * x *2* %pi / T ) + float ( bn ) * sin ( i * x *2* %pi / T ) ; disp ( s ) ; disp ( D i r e c t c u r r e n t : ) ; i = sqrt ( an ^2+ bn ^2) ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14
// e r r o r // q u e s 1 5 , 1 6 , 1 7 // yo = [ 1 . 9 8 1 . 3 0 1 . 0 5 1 . 3 0 0.88 .25 1 . 9 8 ] // x0 =[0 1/6 1/3 1/2 2/3 5/6 1 ] disp ( P r a c t i c a l h a r m o n i c a n a l y s i s ) ; syms x T xo = input ( I n p u t xo m a t r i x ( i n f a c t o r o f T) : ) ; yo = input ( I n p u t yo m a t r i x : ) ; ao =2* sum ( yo ) / length ( xo ) ; s = ao /2; n = input ( No o f s i n o r c o s term i n e x p a n s i o n : ) ; i =1 an =2*( yo .* cos ( i * xo *2* %pi ) ) / length ( yo ) ; bn =2*( yo .* sin ( i * xo *2* %pi ) ) / length ( yo ) ; 83
15 16 17 18 19
s = s + float ( an ) * cos ( i * x *2* %pi / T ) + float ( bn ) * sin ( i * x *2* %pi / T ) ; disp ( s ) ; disp ( D i r e c t c u r r e n t : ) ; i = sqrt ( an ^2+ bn ^2) ;
84
Chapter 13 Linear Dierential Equations
Scilab code Exa 13.1 solvinf linear dierential equation

1 // q u e s 1 2 clc 3 disp ( s o l u t i o n 4 5 6 7 8 9 10 11 12 13 14 15 16
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; syms c1 c2 x m = poly (0 , m ) ; f = m ^2+ m -2; r = roots ( f ) ; disp ( r ) ; y =0; // f o r i =1: l e n g t h ( r ) // syms c ( i ) // y=y+c ( i ) exp ( r ( i ) x ) ; // end y = c1 * exp ( r (1) * x ) + c2 * exp ( r (2) * x ) ; disp ( y= ) ; disp ( y ) ;
Scilab code Exa 13.2 solving linear dierential equation 85
1 // q u e s 2 2 clc 3 disp ( s o l u t i o n 4 5 6 7 8 9 10 11 12
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; syms c1 c2 x ; m = poly (0 , m ) ; f = m ^2+6* m +9; r = roots ( f ) ; disp ( r ) ; disp ( r o o t s a r e e q u a l s o s o l u t i o n i s g i v e n by : ) ; disp ( y= ) ; y =( c1 + x * c2 ) * exp ( r (1) * x ) ; disp ( y ) ;
Scilab code Exa 13.3 solving linear dierential equation

1 // q u e s 4 2 clc 3 disp ( s o l u t i o n 4 5 6 7 8 9 10 11
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; syms c1 c2 c3 x m = poly (0 , m ) ; f = m ^3+ m ^2+4* m +4; r = roots ( f ) ; disp ( r ) ; y = c1 * exp ( r (1) * x ) + c2 * exp ( r (2) * x ) + c3 * exp ( r (3) * x ) ; disp ( y= ) ; disp ( real ( y ) ) ;
Scilab code Exa 13.4 solving linear dierential equation

1 // q u e s 4 2 clc
86
3 4 5 6 7 8 9 10 11
disp ( s o l u t i o n o f t h e g i v e n l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; m = poly (0 , m ) ; syms c1 c2 c3 c4 x f = m ^4+4; r = roots ( f ) ; disp ( r ) ; y = c1 * exp ( r (1) * x ) + c2 * exp ( r (2) * x ) + c3 * exp ( r (3) * x ) + c4 * exp ( r (4) * x ) ; disp ( y= ) ; disp ( real ( y ) ) ;
Scilab code Exa 13.5 nding particular integral

1 // q u e s 5 2 clc 3 disp ( s o l u t i o n 4 5 6 7 8 9
of the given l i n e a r e q u a t i o n i s g i v e n by : ) ; m = poly (0 , m ) ; f = m ^2+5* m +6; // f o r p a r t i c u l a r s o l u t i o n a=1 y = exp ( x ) / horner (f ,1) ; disp ( y ) ; disp ( y ) ;
differential

1 // q u e s 6 2 clc 3 disp ( s o l u t i o n
of the given l i n e a r e q u a t i o n i s g i v e n by : ) ; 4 m = poly (0 , m ) ; 5 f =( m +2) *( m -1) ^2; 87
differential
6 r = roots ( f ) ; 7 disp ( r ) ; 8 disp ( y=1/ f (D) [ exp ( 2x )+exp ( x ) exp ( x ) ) ; 9 disp ( u s i n g 1/ f (D) exp ( ax )=x / f 1 (D) exp ( ax ) i f 10 11 12 13 14 15
f (m) =0
); y1 = x * exp ( -2* x ) /9; y2 = exp ( - x ) /4; y3 = x ^2* exp ( x ) /6; y = y1 + y2 + y3 ; disp ( y= ) ; disp ( y ) ;

1 // q u e s 7 2 clc 3 disp ( s o l u t i o n 4 5 6
7 8 9 10
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; m = poly (0 , m ) ; f = m ^3+1; disp ( U s i n g t h e i d e n t i t y 1/ f (D 2 ) s i n ( ax+b ) [ o r c o s ( ax+b ) ]=1/ f ( a 2 ) s i n ( ax+b ) [ o r c o s ( ax+b ) ] t h i s e q u a t i o n can be r e d u c e d t o ) ; disp ( y =(4D+1) /65 c o s ( 2 x 1) ) ; y =( cos (2* x -1) +4* diff ( cos (2* x -1) ,x ) ) /65; disp ( y= ) ; disp ( y ) ;

1 // q u e s 8 2 clc
88
3 4 5 6 7 8
9 10 11 12
disp ( s o l u t i o n o f t h e g i v e n l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; m = poly (0 , m ) ; f = m ^3+4* m ; disp ( u s i n g 1/ f (D) exp ( ax )=x / f 1 (D) exp ( ax ) i f f (m) =0 ); disp ( y=x 1 / ( 3D2+4) s i n 2 x ) ; disp ( U s i n g t h e i d e n t i t y 1/ f (D 2 ) s i n ( ax+b ) [ o r c o s ( ax+b ) ]=1/ f ( a 2 ) s i n ( ax+b ) [ o r c o s ( ax+b ) ] t h i s e q u a t i o n can be r e d u c e d t o ) ; disp ( y=x /8 s i n 2 x ) ; disp ( y= ) ; y = - x * sin (2* x ) /8; disp ( y ) ;

1 // q u e s 9 2 clc 3 disp ( s o l u t i o n
of the given l i n e a r e q u a t i o n i s g i v e n by : ) ; 4 m = poly (0 , m ) ;

5 6
differential
7 8 9 10 11
disp ( y =1/(D(D+1) ) [ x2+2x +4] can be w r i t t e n a s (1 D+ D 2 ) /D [ x2+2x +4] which i s c o m b i n a t i o n o f d i f f e r e n t i a t i o n and i n t e g r a t i o n ) ; g = x ^2+2* x +4; f =g - diff (g , x ) + diff (g ,x ,2) ; y = integ (f , x ) ; disp ( y= ) ; disp ( y ) ;
Scilab code Exa 13.10 nding particular integral 89
1 // e r r o r 2 clc 3 disp ( s o l u t i o n
of the given l i n e a r e q u a t i o n i s g i v e n by : ) ;
differential
Scilab code Exa 13.11 solving the given linear equation

1 // q u e s 1 1 2 clc 3 disp ( s o l u t i o n 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; disp ( CF + PI ) ; syms c1 c2 x m = poly (0 , m ) ; f =( m -2) ^2; r = roots ( f ) ; disp ( r ) ; disp ( CF i s g i v e n by ) ; cf =( c1 + c2 * x ) * exp ( r (1) * x ) ; disp ( cf ) ; disp ( ) ; disp ( PI =8 { 1/(D 2) 2 [ exp ( 2 x ) ] + { 1 / (D 2) 2 [ s i n ( 2 x ) ] + { 1 / (D 2) 2 [ x 2 ] } ) ; disp ( u s i n g i d e n t i t i e s i t r e d u c e s t o : ) ; pi =4* x ^2* exp (2* x ) + cos (2* x ) +4* x +3; disp ( pi ) ; y = cf + pi ; disp ( The s o l u t i o n i s : y= ) ; disp ( y ) ;

1
// q u e s 1 2 90
2 clc 3 4 disp ( s o l u t i o n 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; disp ( CF + PI ) ; syms c1 c2 x m = poly (0 , m ) ; f =( m ^2 -4) ; r = roots ( f ) ; disp ( r ) ; disp ( CF i s g i v e n by ) ; cf = c1 * exp ( r (1) * x ) + c2 * exp ( r (2) * x ) ; disp ( cf ) ; disp ( ) ; disp ( PI =8 { 1/(D2 4) [ x s i n h ( x ) ] ) ; disp ( u s i n g i d e n t i t i e s i t r e d u c e s t o : ) ; pi = - x /6*( exp ( x ) - exp ( - x ) ) -2/18*( exp ( x ) + exp ( - x ) ) ; disp ( pi ) ; y = cf + pi ; disp ( The s o l u t i o n i s : y= ) ; disp ( y ) ;

1 // q u e s 1 2 2 clc 3 4 disp ( s o l u t i o n 5 6 7 8 9 10
of the given l i n e a r e q u a t i o n i s g i v e n by : ) ; disp ( CF + PI ) ; syms c1 c2 x m = poly (0 , m ) ; f =( m ^2 -1) ; r = roots ( f ) ; disp ( r ) ; 91
differential
11 disp ( CF i s g i v e n by ) ; 12 cf = c1 * exp ( r (1) * x ) + c2 * exp ( r (2) * x ) ; 13 disp ( cf ) ; 14 disp ( ) ; 15 disp ( PI = { 1/(D2 1) [ x s i n ( 3 x )+c o s ( x ) ] ) ; 16 disp ( u s i n g i d e n t i t i e s i t r e d u c e s t o : ) ; 17 pi = -1/10*( x * sin (3* x ) +3/5* cos (3* x ) ) - cos ( x ) /2; 18 disp ( pi ) ; 19 y = cf + pi ; 20 disp ( The s o l u t i o n i s : y= ) ; 21 disp ( y ) ;

1 // q u e s 1 4 2 clc 3 4 disp ( s o l u t i o n 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
of the given l i n e a r d i f f e r e n t i a l e q u a t i o n i s g i v e n by : ) ; disp ( CF + PI ) ; syms c1 c2 c3 c4 x m = poly (0 , m ) ; f =( m ^4+2* m ^2+1) ; r = roots ( f ) ; disp ( r ) ; disp ( CF i s g i v e n by ) ; cf = real (( c1 + c2 * x ) * exp ( r (1) * x ) +( c3 + c4 * x ) * exp ( r (3) * x ) ) ; disp ( cf ) ; disp ( ) ; disp ( PI = { 1/(D4+2 D+1) [ x 2 c o s ( x ) ] ) ; disp ( u s i n g i d e n t i t i e s i t r e d u c e s t o : ) ; pi = -1/48*(( x ^4 -9* x ^2) * cos ( x ) -4* x ^3* sin ( x ) ) ; disp ( pi ) ; y = cf + pi ; 92
20 21
disp ( The s o l u t i o n i s : y= ) ; disp ( y ) ;
93
Chapter 21 Laplace Transform
Scilab code Exa 21.1.1 nding laplace transform

1 2 3 4
// q u e s 1 ( i ) disp ( To f i n d t h e l a p l a c e o f g i v e n f u n c t i o n i n t ) ; syms t s disp ( laplace ( sin (2* t ) * sin (3* t ) ,t , s ) ) ;

1 2 3 4
// q u e s 1 ( i i ) disp ( To f i n d t h e l a p l a c e o f g i v e n f u n c t i o n i n t ) ; syms t s disp ( laplace (( cos ( t ) ) ^2 ,t , s ) ) ;

1 2
// q u e s 1 ( i i ) disp ( To f i n d t h e l a p l a c e o f g i v e n f u n c t i o n i n t ) ; 94
3 4
syms t s disp ( laplace (( sin ( t ) ) ^3 ,t , s ) ) ;

1 // q u e s 1 ( i i ) 2 disp ( To f i n d t h e l a p l a c e o f g i v e n f u n c t i o n 3 syms t s 4 f = exp ( -3* t ) *(2* cos (5* t ) -3* sin (5* t ) ) ; 5 disp ( laplace (f ,t , s ) ) ;
i n t );

1 // q u e s 1 ( i i ) 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s 5 f = exp (3* t ) *( sin ( t ) ) ^2; 6 disp ( laplace (f ,t , s ) ) ;
o f g i v e n f u n c t i o n i n t );

1 // q u e s 1 ( i i ) 2 clc 3 disp ( To f i n d t h e l a p l a c e o f g i v e n 4 syms t s 5 f = exp (4* t ) *( cos ( t ) * sin (2* t ) ) ; 6 disp ( laplace (f ,t , s ) ) ;
f u n c t i o n i n t );
95

1 // q u e s 1 ( i i ) 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = t * sin ( a * t ) ; 6 disp ( laplace (f ,t , s ) ) ;

1 // q u e s 4 ( i i ) 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = t * cos ( a * t ) ; 6 l = laplace (f ,t , s ) ; 7 disp ( l ) ;
Scilab code Exa 21.5 nding laplace transform

1 // e r r o r 2 // q u e s 5 3 clc 4 syms t s u 5 f = integ ( exp ( - s * t ) * t /u ,t ,0 , u ) + integ ( exp ( - s * t ) ,t ,u ,
%inf ) ; 6 disp ( f ) ;
96
Scilab code Exa 21.7 nding laplace transform

1 // q u e s 7 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = sin ( a * t ) / t ; 6 disp ( laplace (f ,t , s ) ) ;

1 // q u e s 7 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = t * cos ( a * t ) ; 6 disp ( laplace (f ,t , s ) ) ;

1 // q u e s 7 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = t ^2* sin ( a * t ) ; 6 disp ( laplace (f ,t , s ) ) ;
97

1 // q u e s 7 2 clc 3 disp ( To f i n d t h e 4 syms t s a 5 f = exp ( -3* t ) * t ^3; 6 l = laplace (f ,t , s ) 7 disp ( l ) ;
l a p l a c e o f g i v e n f u n c t i o n i n t );

1 // q u e s 7 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a 5 f = exp ( - t ) * t * sin (3* t ) ; 6 l = laplace (f ,t , s ) 7 disp ( l ) ;

1 // e r r o r 2 // q u e s 7 3 clc 4 disp ( To f i n d t h e 5 syms t s a 6 f =(1 - exp ( t ) ) / t ; 7 8 l = laplace (f ,t , s ) 9 disp ( l ) ;
l a p l a c e o f g i v e n f u n c t i o n i n t );
98

1 // q u e s 9 2 clc 3 disp ( To f i n d t h e l a p l a c e 4 syms t s a b 5 f =( cos ( a * t ) - cos ( b * t ) ) / t ; 6 7 l = laplace (f ,t , s ) 8 disp ( l ) ;

1 // q u e s 1 0 ( i ) 2 clc 3 disp ( To f i n d t h e t h e g i v e n
i n t e g r a l f i n d the l a p l a c e o f t s i n ( t ) and put s =2 ) ; 4 syms t s m 5 f = sin ( t ) * t ;

6 7 l = laplace (f ,t , s ) 8 s =2 9 10 disp ( eval ( l ) ) ;

1 2
// e r r o r // q u e s 1 0 99
3 clc 4 disp ( To f i n d t h e l a p l a c e o f g i v e n 5 syms t s a b 6 f = integ ( exp ( t ) * sin ( t ) /t ,t ,0 , t ) ; 7 8 l = laplace (f ,t , s ) 9 disp ( l ) ;
f u n c t i o n i n t );
Scilab code Exa 21.11.1 nding inverse laplace transform

1 2 3 4 5 6
// q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t f =( s ^2 -3* s +4) / s ^3; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6
// q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t f =( s +2) /(2* s ^2 -4* s +13) ) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;
100
1 2 3 4 5 6
// q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t f =((2* s ^2 -6* s +5) /( s ^3 -6* s ^2+11* s -6) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6
// q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t f =(4* s +5) /(( s -1) ^2*( s +2) ) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6
// q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t f =(5* s +3) /(( s -1) *( s ^2+2* s +5) ) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;
101
1 2 3 4 5 6 7 8
// e r r o r no o u t p u t // q u e s 1 1 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f = s /( s ^4+4* a ^4) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6 7
// q u e s 1 4 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f = s ^2/( s -2) ^3; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6 7
// q u e s 1 4 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f =( s +3) /(( s ^2 -4* s +13) ) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;
102

1 2 3 4 5 6 7
// no o u t p // q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f =1/( s *( s ^2+ a ^2) ) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;

1 2 3 4 5 6 7
// q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f =1/( s *( s + a ) ^3) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;
Scilab code Exa 21.16.1 nding inverse laplace transform // no o u t p // q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a
1 2 3
103
5 f = s /(( s ^2+ a ^2) ^2) ; 6 il = ilaplace (f ,s , t ) ; 7 disp ( il ) ;

1 2 3 4 5 6 7
// no o u t p u t // q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); syms s t a f = s ^2/(( s ^2+ a ^2) ^2) ; il = ilaplace (f ,s , t ) ; disp ( il ) ;
Scilab code Exa 21.16.3 nding inverse laplace transform // no o u t p u t // q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a
1 2 3 5 6 f =1/(( s ^2+ a ^2) ^2) ; 7 il = ilaplace (f ,s , t ) ; 8 disp ( il ) ;

1
// no o u t p u t 104
// q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a
2 3 5 6 f =( s +2) /( s ^2*( s +1) *( s -2) ) ; 7 il = ilaplace (f ,s , t ) ; 8 disp ( il ) ;
Scilab code Exa 21.17.2 nding inverse laplace transform // no o u t p u t // q u e s 1 5 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a
1 2 3 5 6 f =( s +2) /( s ^2+4* s +5) ^2; 7 il = ilaplace (f ,s , t ) ; 8 disp ( il ) ;
Scilab code Exa 21.19.1 nding inverse laplace transform // e r r o r no o u t p u t // q u e s 1 8 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a
1 2 3 5 6 f = s /( s ^2+ a ^2) ^2; 7 il = ilaplace (f ,s , t ) ; 8 disp ( il ) ;
105
Scilab code Exa 21.19.2 nding inverse laplace transform // e r r o r no o u t p u t // q u e s 1 8 disp ( To f i n d t h e i n v e r s e l a p l a c e t r a n s f o r m o f t h e f u n c t i o n ); 4 syms s t a b
1 2 3 5 6 f = s ^2/(( s ^2+ a ^2) *( s ^2+ b ^2) ) ; 7 il = ilaplace (f ,s , t ) ; 8 disp ( il ) ;

1 // q u e s 2 8 2 syms s t 3 f = integ ( exp ( - s * t ) *( t -1) ,t ,1 ,2) + integ ( exp ( - s * t ) *(3 - t )
,t ,2 ,3) ; 4 disp ( L a p l a c e o f g i v e n f u n c t i o n i s ) ; 5 disp ( f ) ;

1 // q u e s 2 8 2 syms s t 3 f = integ ( exp ( - s * t ) * exp ( - t ) ,t ,0 ,2) ; 4 disp ( L a p l a c e o f g i v e n f u n c t i o n i s ) ; 5 disp ( f ) ;
106
Scilab code Exa 21.34 nding laplace transform // e r r o r no o u t p u t // q u e s 3 4 disp ( t o f i n d t h e l a p l a c e t r a n s f o r m o f p e r i o d i c f u n c t i o n ); 4 syms w t s 5 f =1/(1 - exp ( -2* %pi * s / w ) ) * integ ( exp ( -1* s * t ) * sin ( w * t ) ,t ,0 , %pi / w ) ; 6 disp ( f )
1 2 3
107
Chapter 22 Integral Transform
Scilab code Exa 22.1 nding fourier sine integral

1 // e r r o r 2 // q u e s 1 3 disp ( To f i n d t h e f o u r i e r s i n e i n t e g r a l ) ; 4 syms x t u 5 fs =2/ %pi * integ ( sin ( u * x ) ,u ,0 , %inf ) *( integ ( x ^0* sin ( u * t
) ,t ,0 , %inf ) ) ; 6 disp ( fs ) ;
Scilab code Exa 22.2 nding fourier transform

1 2 3 4 5 6 7 8
// e r r o r // q u e s 2 disp ( To f i n d t h e f o u r i e r t r a n s f o r m o f g i v e n f u n c t i o n ); syms x s F = integ ( exp ( %i * s * x ) ,x , -1 ,1) ; disp ( F ) ; // p r o d u c e s e r r o r > F1 = integ ( sin ( x ) /x ,x ,0 , %inf ) ; 108
Scilab code Exa 22.3 nding fourier transform

1 2 3 4 5 6 7 8
// e r r o r // q u e s 3 disp ( To f i n d t h e f o u r i e r t r a n s f o r m o f g i v e n f u n c t i o n ); syms x s F = integ ( exp ( %i * s * x ) *(1 - x ^2) ,x , -1 ,1) ; disp ( F ) ; // p r o d u c e s e r r o r > F1 = integ (( x * cos ( x ) - sin ( x ) ) / x ^3* cos ( x /2) ,x ,0 , %inf ) ;
Scilab code Exa 22.4 nding fourier sine transform

1 // e r r o r 2 // q u e s 1 3 disp ( To f i n d t h e f o u r i e r s i n e t r a n s f o r m ) ; 4 syms x s m 5 // f u n c t t i o n f s =f ( x ) 6 fs = integ ( sin ( s * x ) * exp ( - x ) ,x ,0 , %inf ) ; 7 disp ( fs ) ; 8 // i n t e g p r o d u c e s e r r o r 9 f = integ ( x * sin ( m * x ) /(1+ x ^2) ,x ,0 , %inf ) ; 10 disp ( f ) ;
Scilab code Exa 22.5 nding fourier cosine transform

1 2
// q u e s 5 syms x s 109
3 disp ( F o u r i e r c o s i n e t r a n s f o r m ) ; 4 f = integ ( x * cos ( s * x ) ,x ,0 ,1) + integ ((2 - x ) * cos ( s * x ) ,x
,1 ,2) ; 5 disp ( f )
Scilab code Exa 22.6 nding fourier sine transform

1 // q u e s 6 2 syms x s a 3 disp ( F o u r i e r c o s i n e t r a n s f o r m ) ; 4 f = integ ( exp ( - a * x ) / x * sin ( s * x ) ,x ,0 , %inf ) ; 5 disp ( f )
110
Chapter 23 Statistical Methods
Scilab code Exa 23.1 Calculating cumulative frequencies of given using iterations on matrices
1 clear 2 clc 3 disp ( t h e
4 5 6 7 8 9 10 11 12 13 14 15
f i r s t row o f A d e n o t e s t h e no . o f s t u d e n t s f a l l i n g i n t h e marks g r o u p s t a r t i n g from (5 10) . . . t i l l (40 45) ) A (1 ,:) =[5 6 15 10 5 4 2 2]; disp ( t h e s e c o n d row d e n o t e s c u m u l a t i v e f r e q u e n c y ( l e s s than ) ) A (2 ,1) =5; for i =2:8 A (2 , i ) = A (2 ,i -1) + A (1 , i ) ; end disp ( t h e t h i r d row d e n o t e s c u m u l a t i v e f r e q u e n c y ( more t h a n ) ) A (3 ,1) =49; for i =2:8 A (3 , i ) = A (3 ,i -1) -A (1 ,i -1) ; end disp ( A )
111
Scilab code Exa 23.2 Calculating mean of of statistical data performing iterations matrices
1 clc 2 disp ( t h e
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
f i r s t row o f A r e p r e s e n t s t h e mid v a l u e s o f weekly e a r n i n g s having i n t e r v a l o f 2 in each c l a s s =x ) A (1 ,:) =[11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41] disp ( t h e s e c o n d row d e n o t e s t h e no . o f e m p l o y e e s o r i n o t h e r words f r e q u e n c y=f ) A (2 ,:) =[3 6 10 15 24 42 75 90 79 55 36 26 19 13 9 7] disp ( t h i r d row d e n o t e s f x ) for i =1:16 A (3 , i ) = A (1 , i ) * A (2 , i ) ; end disp ( f o u r t h row d e n o t e s u=(x 25) /2 ) for i =1:16 A (4 , i ) =( A (1 , i ) -25) /2 end disp ( f i f t h row d e n o t e s f x ) for i =1:16 A (5 , i ) = A (4 , i ) * A (2 , i ) ; end A b =0; disp ( sum o f a l l e l e m e n t s o f t h i r d row= ) for i =1:16 b += A (3 , i ) end disp ( b ) f =0; disp ( sum o f a l l e l e m e n t s o f s e c o n d row= ) for i =1:16 112
28 29 30 31 32 33 34 35 36 37 38 39
f += A (2 , i ) end disp ( f ) disp ( mean= ) b/f d =0; disp ( sum o f a l l e l e m e n t s o f f i f t h row= ) for i =1:16 d += A (5 , i ) end disp ( mean by s t e p d e v i a t i o n method= ) 25+(2* d / f )
Scilab code Exa 23.3 Analysis of statistical data performing iterations on matrices
4 5 6 7 8 9 10 11 12 13
f i r s t row o f A d e n o t e s t h e no . o f s t u d e n t s f a l l i n g i n t h e marks g r o u p s t a r t i n g from (5 10) . . . t i l l (40 45) ) A (1 ,:) =[5 6 15 10 5 4 2 2]; disp ( t h e s e c o n d row d e n o t e s c u m u l a t i v e f r e q u e n c y ( l e s s than ) ) A (2 ,:) =[5 11 26 36 41 45 47 49] disp ( t h e t h i r d row d e n o t e s c u m u l a t i v e f r e q u e n c y ( more t h a n ) ) A (3 ,:) =[49 44 38 23 13 8 4 2] disp ( median f a l l s i n t h e c l a s s (15 20) = l +(( n/2 c ) h ) / f= ) 15+((49/2 -11) *5) /15 disp ( l o w e r q u a r t i l e a l s o f a l l s i n t h e c l a s s (15 20) =) Q1 =15+((49/4 -11) *5) /15 disp ( u p p e r q u a r t i l e a l s o f a l l s i n t h e c l a s s (25 30) 113
=) 14 Q3 =25+((3*49/4 -36) *5) /5 15 disp ( s e m i i n t e r q u a r t i l e r a n g e= ) 16 ( Q3 - Q1 ) /2
Scilab code Exa 23.4 Analysis of statistical data

4 5 6 7 8 9 10 11 12 13 14 15
f i r s t row o f A d e n o t e s t h e r o l l no . o f s t u d e n t s form 1 t o 10 and t h a t o f B d e n o t e s form 11 t o 20 ) A (1 ,:) =[1 2 3 4 5 6 7 8 9 10]; B (1 ,:) =[11 12 13 14 15 16 17 18 19 20]; disp ( t h e s e c o n d row o f A annd B d e n o t e s t h e c o r r e s p o n d i n g marks i n p h y s i c s ) A (2 ,:) =[53 54 52 32 30 60 47 46 35 28]; B (2 ,:) =[25 42 33 48 72 51 45 33 65 29]; disp ( t h e t h i r d row d e n o t e s t h e c o r r e s p o n d i n g marks in chemistry ) A (3 ,:) =[58 55 25 32 26 85 44 80 33 72]; B (3 ,:) =[10 42 15 46 50 64 39 38 30 36]; disp ( median marks i n p h y s i c s =a r i t h m e t i c mean o f 10 thand 11 t h s t u d e n t = ) (28+25) /2 disp ( median marks i n c h e m i s t r y =a r i t h m e t i c mean o f 10 thand 11 t h s t u d e n t = ) (72+10) /2
Scilab code Exa 23.5 Finding the missing frequency of given statistical data using given constants
1
clear 114
2 3 4 5 6 7 8 9 10 11
clc disp ( l e t t h e m i s s s i n g f r e q u e n c i e s be f 1 a n d f 2 ) disp ( sum o f g i v e n f r e q u e n c i e s =12+30+65+25+18= ) c =12+30+65+25+18 disp ( so , f 1+f 2 =229 c= ) 229 - c disp ( median =46=40+(114.5 (12+30+ f 1 ) ) 1 0 / 6 5 ) ) disp ( f 1 =33.5=34 ) f1 =34 f2 =45
Scilab code Exa 23.6 Calculating average speed

1 2 3 4 5 6 7 8
clear clc syms s ; disp ( l e t t h e e q i d i s t a n c e be s , t h e n ) t1 = s /30 t2 = s /40 t3 = s /50 disp ( a v e r a g e s p e e d= t o t a l d i s t a n c e / t o t a l t i m e t a k e n ) 9 3* s /( t1 + t2 + t3 )
Scilab code Exa 23.7 Calculating mean and standard deviation performing iterations on matrices
1 clear 2 clc 3 disp ( t h e f i r s t row d e n o t e s t h e s i z e o f i t e m ) 4 A (1 ,:) =[6 7 8 9 10 11 12]; 5 disp ( t h e s e c o n d row d e n o t e s t h e c o r r e s p o n d i n g
frequency ( f ) ) 115
6 A (2 ,:) =[3 6 9 13 8 5 4]; 7 disp ( t h e t h i r d row d e n o t e s t h e c o r r e s p o n d i n g 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
deviation (d) ) A (3 ,:) =[ -3 -2 -1 0 1 2 3]; disp ( t h e f o u r t h row d e n o t e s t h e c o r r e s p o n d i n g f d ) for i =1:7 A (4 , i ) = A (2 , i ) * A (3 , i ) ; end disp ( t h e f i f t h row d e n o t e s t h e c o r r e s p o n d i n g f d 2 ) for i =1:7 A (5 , i ) = A (2 , i ) *( A (3 , i ) ^2) ; end A b =0; disp ( sum o f f o u r t h row e l e m e n t s= ) for i =1:7 b = b + A (4 , i ) ; end disp ( b ) c =0 disp ( sum o f f i f t h row e l e m e n t s= ) for i =1:7 c = c + A (5 , i ) ; end disp ( c ) d =0; disp ( sum o f a l l f r e q u e n c i e s = ) for i =1:7 d = d + A (2 , i ) ; end disp ( d ) disp ( mean=9+b /d= ) 9+ b / d disp ( s t a n d a r d d e v i a t i o n =( c / d ) 0 . 5 ) ( c / d ) ^0.5
116
Scilab code Exa 23.8 Calculating mean and standard deviation performing iterations on matrices
1 clc 2 disp ( t h e
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
f i r s t row o f A r e p r e s e n t s t h e mid v a l u e s o f wage c l a s s e s h a v i n g i n t e r v a l o f 8 i n e a c h c l a s s =x ) A (1 ,:) =[8.5 16.5 24.5 32.5 40.5 48.5 56.5 64.5 72.5] disp ( t h e s e c o n d row d e n o t e s t h e no . o f men o r i n o t h e r words f r e q u e n c y=f ) A (2 ,:) =[2 24 21 18 5 3 5 8 2] disp ( t h i r d row d e n o t e s f x ) for i =1:9 A (3 , i ) = A (1 , i ) * A (2 , i ) ; end disp ( f o u r t h row d e n o t e s d=(x 3 2 . 5 ) /8 ) for i =1:9 A (4 , i ) =( A (1 , i ) -32.5) /8 end disp ( f i f t h row d e n o t e s f d ) for i =1:9 A (5 , i ) = A (4 , i ) * A (2 , i ) ; end disp ( s i x t h row d e n o t e s f ( d 2 ) ) for i =1:9 A (6 , i ) = A (4 , i ) ^2* A (2 , i ) ; end A b =0; disp ( sum o f a l l e l e m e n t s o f s i x t h row= ) for i =1:9 b += A (6 , i ) end disp ( b ) 117
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
f =0; disp ( sum o f a l l e l e m e n t s o f s e c o n d row= ) for i =1:9 f += A (2 , i ) end disp ( f ) disp ( mean= ) b/f d =0; disp ( sum o f a l l e l e m e n t s o f f i f t h row= ) for i =1:9 d += A (5 , i ) end disp ( mean wage= ) 32.5+(8* d / f ) disp ( s t a n d a r d d e v i a t i o n= ) 8*( b /f -( d / f ) ^2)
Scilab code Exa 23.9 Analysis of statistical data performing iterations on matrices
1 clear 2 clc 3 disp ( t h e 4 5 6 7 8 9 10 11
f i r s t row o f A d e n o t e s t h e s c o r e s o f A and t h a t o f B d e n o t e s t h a t o f B ) A (1 ,:) =[12 115 6 73 7 19 119 36 84 29]; B (1 ,:) =[47 12 16 42 4 51 37 48 13 0]; disp ( t h e s e c o n d row o f A annd B d e n o t e s t h e corresponding deviation ) for i =1:10 A (2 , i ) = A (1 , i ) -51; B (2 , i ) = B (1 , i ) -51; end disp ( t h e t h i r d row o f A and B d e n o t e s t h e corresponding deviation square ) 118
12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49
for i =1:10 A (3 , i ) = A (2 , i ) ^2; B (3 , i ) = B (2 , i ) ^2; end A B b =0; disp ( sum o f s e c o n d row e l e m e n t s o f for i =1:10 b = b + A (2 , i ) ; end disp ( b ) c =0; disp ( sum o f s e c o n d row e l e m e n t s o f for i =1:10 c = c + B (2 , i ) ; end disp ( c ) d =0; disp ( sum o f t h i r d row e l e m e n t s o f for i =1:10 d = d + A (3 , i ) ; end disp ( d ) e =0; disp ( sum o f s e c o n d row e l e m e n t s o f for i =1:10 e = e + B (3 , i ) ; end disp ( e ) disp ( a r i t h m e t i c mean o f A= ) f =51+ b /10 disp ( s t a n d a r d d e v i a t i o n o f A= ) g =( d /10 -( b /10) ^2) ^0.5 disp ( a r i t h m e t i c mean o f B= ) h =51+ c /10 disp ( s t a n d a r d d e v i a t i o n o f A= ) i =( e /10 -( c /10) ^2) ^0.5 119
A=b= )
B=c= )
A=d= )
B=e= )
50 51 52 53
disp ( c o e f f i c i e n t o f v a r i a t i o n o f A= ) ( g / f ) *100 disp ( c o e f f i c i e n t o f v a r i a t i o n o f B= ) ( i / h ) *100
Scilab code Exa 23.10 Calculating mean and standard deviation of different statistical data when put together
1 clear 2 clc 3 disp ( i f m i s t h e mean o f e n t i r e data , t h e n ) 4 m =(50*113+60*120+90*115) /(50+60+90) 5 disp ( i f s i s t h e s t a n d a r d d e v i a t i o n o f e n t i r e data ,
then ) 6 s =(((50*6^2) +(60*7^2) +(90*8^2) +(50*3^2) +(60*4^2) +(90*1^2) ) /200) ^0.5
Scilab code Exa 23.12 Calculating median and quartiles of given statistical data performing iterations on matrices
4 5 6 7 8 9 10
f i r s t row o f A d e n o t e s t h e no . o f p e r s o n s f a l l i n g i n t h e w e i g h t g r o u p s t a r t i n g from (70 80) . . . t i l l (140 150) ) A (1 ,:) =[12 18 35 42 50 45 20 8]; disp ( t h e s e c o n d row d e n o t e s c u m u l a t i v e f r e q u e n c y ) A (2 ,1) =12; for i =2:8 A (2 , i ) = A (2 ,i -1) + A (1 , i ) ; end disp ( median f a l l s i n t h e c l a s s (110 120) = l +(( n/2 c ) h ) / f= ) 120
11 Q2 =110+(8*10) /50 12 disp ( l o w e r q u a r t i l e 13 14 15 16 17
a l s o f a l l s in the c l a s s (90 100)= ) Q1 =90+(57.5 -30) *10/35 disp ( u p p e r q u a r t i l e a l s o f a l l s i n t h e c l a s s (120 130)= ) Q3 =120+(172.5 -157) *10/45 disp ( q u a r t i l e c o e f f i c i e n t o f s k e w n e s s= ) ( Q1 + Q3 -2* Q2 ) /( Q3 - Q1 )
Scilab code Exa 23.13 Calculating coecient of correlation

1 clear 2 clc 3 disp ( t h e 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
f i r s t row o f A d e n o t e s t h e c o r r e s p o n d i n g I .R. of students ) A (1 ,:) =[105 104 102 101 100 99 98 96 93 92]; disp ( t h e s e c o n d row d e n o t e s t h e c o r r e s p o n d i n g d e v i a t i o n of I .R. ) for i =1:10 A (2 , i ) = A (1 , i ) -99; end disp ( t h e t h i r d row d e n o t e s t h e s q u a r e o f corresponding d e v i a t i o n of I .R. ) for i =1:10 A (3 , i ) = A (2 , i ) ^2; end disp ( t h e f o u r t h row d e n o t e s t h e c o r r e s p o n d i n g E . R . of students ) A (4 ,:) =[101 103 100 98 95 96 104 92 97 94]; disp ( t h e f i f t h row d e n o t e s t h e c o r r e s p o n d i n g d e v i a t i o n of E.R. ) for i =1:10 A (5 , i ) = A (4 , i ) -98; end 121
19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
disp ( t h e s i x t h row d e n o t e s t h e s q u a r e o f corresponding d e v i a t i o n of E.R. ) for i =1:10 A (6 , i ) = A (5 , i ) ^2; end disp ( t h e s e v e n t h row d e n o t e s t h e p r o d u c t o f t h e two corresponding deviations ) for i =1:10 A (7 , i ) = A (2 , i ) * A (5 , i ) ; end A a =0; disp ( t h e sum o f e l e m e n t s o f f i r s t row=a ) for i =1:10 a = a + A (1 , i ) ; end a b =0; disp ( t h e sum o f e l e m e n t s o f s e c o n d row=b ) for i =1:10 b = b + A (2 , i ) ; end b c =0; disp ( t h e sum o f e l e m e n t s o f t h i r d row=c ) for i =1:10 c = c + A (3 , i ) ; end c d =0; disp ( t h e sum o f e l e m e n t s o f f o u r t h row=d ) for i =1:10 d = d + A (4 , i ) ; end d e =0; disp ( t h e sum o f e l e m e n t s o f f i f t h row=e ) for i =1:10 122
55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71
e = e + A (5 , i ) ; end e f =0; disp ( t h e sum o f e l e m e n t s o f s i x t h row=d ) for i =1:10 f = f + A (6 , i ) ; end f g =0; disp ( t h e sum o f e l e m e n t s o f s e v e n t h row=d ) for i =1:10 g = g + A (7 , i ) ; end g disp ( c o e f f i c i e n t o f c o r r e l a t i o n = ) g /( c * f ) ^0.5
123
Chapter 24 Numerical Methods
Scilab code Exa 24.1 nding the roots of equation

1 clc 2 clear 3 x = poly (0 , x ) ; 4 p = x ^3 -4* x -9 5 disp ( F i n d i n g r o o t s 6 7 8 9 10 11 12 13 14 15 16 17 18 19
o f t h i s e q u a t i o n by b i s e c t i o n method ) ; disp ( f ( 2 ) i s ve and f ( 3 ) i s +ve s o a r o o t l i e s b e t w e e n 2 and 3 ) ; l =2; m =3; function y = f ( x ) y = x ^3 -4* x -9; endfunction for i =1:4 k =1/2*( l + m ) ; if ( f ( k ) <0) l=k; else m=k; end end 124
20
disp ( k )
Scilab code Exa 24.3 nding the roots of equation by the method of false statement
1 2 3 4 5 6 7
// q u e s 2 disp ( f ( x )=xe x c o s ( x ) ) ; function y = f ( x ) y = x * %e ^( x ) - cos ( x ) ; endfunction
disp ( we a r e r e q u i r e d t o f i n d t h e r o o t s o f f ( x ) by t h e method o f f a l s e p o s i t i o n ) ; 8 disp ( f ( 0 )=ve and f ( 1 )=+ve s o s r o o t l i e b e t w e e n 0 and 1 ) ; 9 disp ( f i n d i n g t h e r o o t s by f a l s e p o s i t i o n method ) ;

10 11 12 13 14 15 16 17 18 19 20 21 22 23
l =0; m =1; for i =1:10 k =l -( m - l ) * f ( l ) /( f ( m ) -f ( l ) ) ; if ( f ( k ) <0) l=k; else m=k; end end // f p r i n t f ( The r o o t s o f t h e e q u a t i o n i s %g , k ) disp ( The r o o t o f t h e e q u a t i o n i s : ) ; disp ( k ) ;
Scilab code Exa 24.4 nding rea roots of equation by regula falsi method 125
1 2 3 4 5 6 7
// q u e s 2 disp ( f ( x )=x l o g ( x ) 1.2 ) ; function y = f ( x ) y = x * log10 ( x ) -1.2; endfunction
disp ( we a r e r e q u i r e d t o f i n d t h e r o o t s o f f ( x ) by t h e method o f f a l s e p o s i t i o n ) ; 8 disp ( f ( 2 )=ve and f ( 3 )=+ve s o s r o o t l i e b e t w e e n 2 and 3 ) ; 9 disp ( f i n d i n g t h e r o o t s by f a l s e p o s i t i o n method ) ; l =2; m =3; for i =1:3 k =l -( m - l ) * f ( l ) /( f ( m ) -f ( l ) ) ; if ( f ( k ) <0) l=k; else m=k; end end // f p r i n t f ( The r o o t s o f t h e e q u a t i o n i s %g , k ) disp ( The r o o t o f t h e e q u a t i o n i s : ) ; disp ( k ) ;
10 11 12 13 14 15 16 17 18 19 20 21 22 23
Scilab code Exa 24.5 real roots of equation by newtons method // q u e s 5 disp ( To f i n d t h e r o o t s o f f ( x ) =3x c o s ( x ) 1 by n e w t o n s method ) ; 3 disp ( f ( 0 )=ve and f ( 1 ) i s +ve s o a r o o t l i e s b e t w e e n 0 and 1 ) ; 4 l =0; 5 m =1;
1 2
126
6 function y = f ( x ) 7 y =3* x - cos ( x ) -1; 8 endfunction 9 x0 =0.6; 10 disp ( l e t u s t a k e x0 =0.6 a s t h e r o o t 11 12 13 14 15 16 17
i s c l o s e r to 1
); disp ( Root i s g i v e n by r=x0 f ( xn ) / d e r ( f ( xn ) ) ) ; disp ( a p p r o x i m a t e d r o o t i n e a c h s t e p s a r e ) ; for i =1:3 k = x0 - f ( x0 ) / derivative (f , x0 ) ; disp ( k ) ; x0 = k ; end
Scilab code Exa 24.6 real roots of equation by newtons method

1 // q u e s 7 2 clear 3 clc 4 disp ( To f i n d 5 6 7 8 9 10 11 12 13 14 15 16 17
s q u a r e r o o t o f 28 by n e w t o n s method l e t x=s q r t ( 2 8 ) i e x2 28=0 ) ; function y = f ( x ) y = x ^2 -28; endfunction disp ( To f i n d t h e r o o t s by n e w t o n s method ) ; disp ( f ( 5 )=ve and f ( 6 ) i s +ve s o a r o o t l i e s b e t w e e n 5 and 6 ) ; l =5; m =6; disp ( l e t u s t a k e x0 =5.5 ) ; disp ( Root i s g i v e n by r n=xn f ( xn ) / d e r ( f ( xn ) ) ) ; disp ( a p p r o x i m a t e d r o o t i n e a c h s t e p s a r e ) ; x0 =5.5; for i =1:4 k = x0 - f ( x0 ) / derivative (f , x0 ) ; 127
18 disp ( k ) ; 19 x0 = k ; 20 end
Scilab code Exa 24.7 evaluating square root by newtons iterative method
1 // q u e s 7 2 clear 3 clc 4 disp ( To f i n d 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
s q u a r e r o o t o f 28 by n e w t o n s method l e t x=s q r t ( 2 8 ) i e x2 28=0 ) ; function y = f ( x ) y = x ^2 -28; endfunction disp ( To f i n d t h e r o o t s by n e w t o n s method ) ; disp ( f ( 5 )=ve and f ( 6 ) i s +ve s o a r o o t l i e s b e t w e e n 5 and 6 ) ; l =5; m =6; disp ( l e t u s t a k e x0 =5.5 ) ; disp ( Root i s g i v e n by r n=xn f ( xn ) / d e r ( f ( xn ) ) ) ; disp ( a p p r o x i m a t e d r o o t i n e a c h s t e p s a r e ) ; x0 =5.5; for i =1:4 k = x0 - f ( x0 ) / derivative (f , x0 ) ; disp ( k ) ; x0 = k ; end
Scilab code Exa 24.10 solving equations by guass elimination method

1
// q u e s 10 , q u e s 11 128
// L i n e a r e q u a t i o n s y s t e m Ax=r by Gauss e l i m i n a t i o n method . 3 clc 4 clear

2 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
disp ( S o l u t i o n o f N e q u a t i o n [ A ] [ X] = [ r ] ) n = input ( E n t e r number o f E q u a t i o n s : ) ; A = input ( E n t e r M a t r i x [ A ] : ) ; r = input ( E n t e r M a t r i x [ r ] : ) ; D=A;d=r; // c r e a t e u p p e r t r i a n g u l a r m a t r i x s =0; for j =1: n -1 if A (j , j ) ==0 k=j; for k = k +1: n if A (k , j ) ==0 continue end break end B = A (j ,:) ; C = r ( j ) ; A (j ,:) = A (k ,:) ; r ( j ) = r ( k ) ; A (k ,:) = B ; r ( k ) = C ; end for i =1+ s :n -1 L = A ( i +1 , j ) / A (j , j ) ; A ( i +1 ,:) = A ( i +1 ,:) -L * A (j ,:) ; r ( i +1) = r ( i +1) -L * r ( j ) ; end s = s +1; end // S o l u t i o n o f e q u a t i o n s x ( n ) = r ( n ) / A (n , n ) ; for i =n -1: -1:1 sum =0; for j = i +1: n 129
39 sum = sum + A (i , j ) * x ( j ) ; 40 end 41 x ( i ) =(1/ A (i , i ) ) *( r ( i ) - sum ) ; 42 end 43 44 // h e c k i n g w i t h s c i l a b f u n c t i o n s 45 p = inv ( D ) * d ; 46 // Output 47 disp ( @
48 49 50 51 52
@) disp ( Output [ B ] [ x ] = [ b ] ) disp ( Upper r i a n g u l a r M a t r i x [ B ] = ) ; disp ( A ) disp ( M a t r i x [ b ] = ) ; disp ( r ) disp ( s o l u t i o n o f l i n e a r e q u a t i o n s : ) ; disp (x ) disp ( s o l v e w i t h m a t l a b f u n c t i o n s ( f o r c h e c k i n g ) : ) ; disp ( p )
Scilab code Exa 24.12 solving equations by guass elimination method // q u e s 10 , q u e s 11 // L i n e a r e q u a t i o n s y s t e m Ax=r by Gauss e l i m i n a t i o n method . 3 clc 4 clear
1 2 5 6 7 8 9 10 11 12 13 14
disp ( S o l u t i o n o f N e q u a t i o n [ A ] [ X] = [ r ] ) n = input ( E n t e r number o f E q u a t i o n s : ) ; A = input ( E n t e r M a t r i x [ A ] : ) ; r = input ( E n t e r M a t r i x [ r ] : ) ; D=A;d=r; // c r e a t e u p p e r t r i a n g u l a r m a t r i x s =0; for j =1: n -1 130
15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
if A (j , j ) ==0 k=j; for k = k +1: n if A (k , j ) ==0 continue end break end B = A (j ,:) ; C = r ( j ) ; A (j ,:) = A (k ,:) ; r ( j ) = r ( k ) ; A (k ,:) = B ; r ( k ) = C ; end for i =1+ s :n -1 L = A ( i +1 , j ) / A (j , j ) ; A ( i +1 ,:) = A ( i +1 ,:) -L * A (j ,:) ; r ( i +1) = r ( i +1) -L * r ( j ) ; end s = s +1; end // S o l u t i o n o f e q u a t i o n s x ( n ) = r ( n ) / A (n , n ) ; for i =n -1: -1:1 sum =0; for j = i +1: n sum = sum + A (i , j ) * x ( j ) ; end x ( i ) =(1/ A (i , i ) ) *( r ( i ) - sum ) ; end
// h e c k i n g w i t h s c i l a b f u n c t i o n s p = inv ( D ) * d ; // Output disp ( @ @) 48 disp ( Output [ B ] [ x ] = [ b ] ) 49 disp ( Upper r i a n g u l a r M a t r i x [ B ] = ) ; disp ( A ) 50 disp ( M a t r i x [ b ] = ) ; disp ( r ) 131
51 52
disp ( s o l u t i o n o f l i n e a r e q u a t i o n s : ) ; disp (x ) disp ( s o l v e w i t h m a t l a b f u n c t i o n s ( f o r c h e c k i n g ) : ) ; disp ( p )
Scilab code Exa 24.13 solving equations by guass elimination method // q u e s 10 , q u e s 11 // L i n e a r e q u a t i o n s y s t e m Ax=r by Gauss e l i m i n a t i o n method . 3 clc 4 clear
1 2 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
disp ( S o l u t i o n o f N e q u a t i o n [ A ] [ X] = [ r ] ) n = input ( E n t e r number o f E q u a t i o n s : ) ; A = input ( E n t e r M a t r i x [ A ] : ) ; r = input ( E n t e r M a t r i x [ r ] : ) ; D=A;d=r; // c r e a t e u p p e r t r i a n g u l a r m a t r i x s =0; for j =1: n -1 if A (j , j ) ==0 k=j; for k = k +1: n if A (k , j ) ==0 continue end break end B = A (j ,:) ; C = r ( j ) ; A (j ,:) = A (k ,:) ; r ( j ) = r ( k ) ; A (k ,:) = B ; r ( k ) = C ; end for i =1+ s :n -1 L = A ( i +1 , j ) / A (j , j ) ; 132
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
A ( i +1 ,:) = A ( i +1 ,:) -L * A (j ,:) ; r ( i +1) = r ( i +1) -L * r ( j ) ; end s = s +1; end // S o l u t i o n o f e q u a t i o n s x ( n ) = r ( n ) / A (n , n ) ; for i =n -1: -1:1 sum =0; for j = i +1: n sum = sum + A (i , j ) * x ( j ) ; end x ( i ) =(1/ A (i , i ) ) *( r ( i ) - sum ) ; end // h e c k i n g w i t h s c i l a b f u n c t i o n s p = inv ( D ) * d ; // Output disp ( @ @) disp ( Output [ B ] [ x ] = [ b ] ) disp ( Upper r i a n g u l a r M a t r i x [ B ] = ) ; disp ( A ) disp ( M a t r i x [ b ] = ) ; disp ( r ) disp ( s o l u t i o n o f l i n e a r e q u a t i o n s : ) ; disp (x ) disp ( s o l v e w i t h m a t l a b f u n c t i o n s ( f o r c h e c k i n g ) : ) ; disp ( p )
48 49 50 51 52
133
Chapter 26 Dierence Equations and Z Transform
Scilab code Exa 26.2 nding dierence equation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
// q u e s 2 syms n a b yn0 yn1 yn2 yn = a *2^ n + b *( -2) ^ n ; disp ( yn= ) ; disp ( yn ) ; n = n +1; yn = eval ( yn ) ; disp ( y ( n+1)=yn1= ) ; disp ( yn ) ; n = n +1; yn = eval ( yn ) ; disp ( y ( n+2)=yn2= ) ; disp ( yn ) ; disp ( E l i m i n a t i n g a b fropm t h e s e e q u a t i o n s we g e t : ); A =[ yn0 1 1; yn1 2 -2; yn2 4 4] y = det ( A ) ; disp ( The r e q u i r e d d i f f e r e n c e e q u a t i o n : ) ; disp ( y ) ; 134
19
disp ( =0 ) ;
Scilab code Exa 26.3 solving dierence equation

1 2 3 4 5 6 7 8 9 10 11
// q u e s 3 syms c1 c2 c3 disp ( C u m u l a t i v e f u n c t i o n i s g i v e n by E3 2 E2 5 E +6 =0 ) ; E = poly (0 , E ) ; f = E ^3 -2* E ^2 -5* E +6; r = roots ( f ) ; disp ( r ) ; disp ( There f o r t h e c o m p l e t e s o l u t i o n i s : ) ; un = c1 *( r (1) ) ^ n + c2 *( r (2) ) ^ n + c3 *( r (3) ) ^ n ; disp ( un= ) ; disp ( un ) ;

1 2 3 4 5 6 7 8 9 10 11
// q u e s 4 syms c1 c2 c3 n disp ( C u m u l a t i v e f u n c t i o n i s g i v e n by E2 2 E+1 =0 ) ; E = poly (0 , E ) ; f = E ^2 -2* E +1; r = roots ( f ) ; disp ( r ) ; disp ( There f o r t h e c o m p l e t e s o l u t i o n i s : ) ; un =( c1 + c2 * n ) *( r (1) ) ^ n ; disp ( un= ) ; disp ( un ) ;
135
Scilab code Exa 26.6 rming bonacci dierence equation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
// q u e s 6 syms c1 c2 c3 n disp ( For F i b o n a c c i S e r i e s yn2=yn1+yn0 ) ; disp ( s o C u m u l a t i v e f u n c t i o n i s g i v e n by E2 E1 =0 ) ; E = poly (0 , E ) ; f = E ^2 -E -1; r = roots ( f ) ; disp ( r ) ; disp ( There f o r t h e c o m p l e t e s o l u t i o n i s : ) ; un =( c1 ) *( r (1) ) ^ n + c2 *( r (2) ) ^ n ; disp ( un= ) ; disp ( un ) ; disp ( Now p u t t t i n g n =1 , y=0 and n=2 , y=1 we g e t ) ; disp ( c 1 =(5 s q r t ( 5 ) ) /10 c 2 =(5+ s q r t ( 5 ) ) /10 ) ; c1 =(5 - sqrt (5) ) /10; c2 =(5+ sqrt (5) ) /10; un = eval ( un ) ; disp ( un ) ;

1 2 3 4 5 6 7
// q u e s 4 syms c1 c2 c3 n disp ( C u m u l a t i v e f u n c t i o n i s g i v e n by E2 4 E+3 =0 ) ; E = poly (0 , E ) ; f = E ^2 -4* E +3; r = roots ( f ) ; disp ( r ) ; 136
8 9 10 11 12 13 14 15 16 17 18
disp ( There f o r t h e c o m p l e t e s o l u t i o n i s = c f + p i ) ; cf = c1 *( r (1) ) ^ n + c2 * r (2) ^ n ; disp ( CF= ) ; disp ( cf ) ; disp ( PI = 1 / (E2 4E+3) [ 5 n ] ) ; disp ( put E=5 ) ; disp ( We g e t PI=5n /8 ) ; pi =5^ n /8; un = cf + pi ; disp ( un= ) ; disp ( un ) ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
// q u e s 4 syms c1 c2 c3 n disp ( C u m u l a t i v e f u n c t i o n i s g i v e n by E2 4 E+4 =0 ) ; E = poly (0 , E ) ; f = E ^2 -4* E +4; r = roots ( f ) ; disp ( r ) ; disp ( There f o r t h e c o m p l e t e s o l u t i o n i s = c f + p i ) ; cf =( c1 + c2 * n ) * r (1) ^ n ; disp ( CF= ) ; disp ( cf ) ; disp ( PI = 1 / (E2 4E+4) [ 2 n ] ) ; disp ( We g e t PI=n ( n 1) / 2 2 ( n 2) ) ; pi = n *( n -1) / factorial (2) *2^( n -2) ; un = cf + pi ; disp ( un= ) ; disp ( un ) ;
137

1 // q u e s 1 0 2 clc 3 syms c1 c2 c3 n 4 disp ( C u m u l a t i v e f u n c t i o n
i s g i v e n by E2 4
=0 )
;
5 E = poly (0 , E ) ; 6 f = E ^2 -4; 7 r = roots ( f ) ; 8 disp ( r ) ; 9 disp ( There f o r t h e c o m p l e t e 10 11 12 13 14 15 16 17 18 19
solution i s = cf + pi ) ; cf =( c1 + c2 * n ) * r (1) ^ n ; disp ( CF= ) ; disp ( cf ) ; // p a r t i c u l a r i n t e g r a l c a l u l a t i o n m a n u a l l y disp ( PI = 1 / (E2 4) [ n2+n 1] ) ; disp ( We g e t PI=n 2/3 7/9 n 17/27 ) ; pi = - n ^2/3 -7/9* n -17/27; un = cf + pi ; disp ( un= ) ; disp ( un ) ;

1 // q u e s 1 1 2 clc 3 syms c1 c2 c3 n 4 disp ( C u m u l a t i v e f u n c t i o n
i s g i v e n by E2 2 E+1
=0 ) ; 5 E = poly (0 , E ) ; 138
6 f = E ^2+2* E -1; 7 r = roots ( f ) ; 8 disp ( r ) ; 9 disp ( There f o r t h e c o m p l e t e 10 11 12 13 14 15 16 17 18 19
solution i s = cf + pi ) ; cf =( c1 + c2 * n ) * r (1) ^ n ; disp ( CF= ) ; disp ( cf ) ; // p a r t i c u l a r i n t e g r a l c a l u l a t i o n m a n u a l l y disp ( PI = 1 / ( E 1) 2 [ n 2 2 n ] ) ; disp ( We g e t PI=2n ( n2 8 n+20 ) ; pi =2^ n *( n ^2 -8* n +20) ; un = cf + pi ; disp ( un= ) ; disp ( un ) ;
Scilab code Exa 26.12 solving simultanious dierence equation

1 // q u e s 1 2 2 clc 3 disp ( s i m p l i f i e d e q u a t i o n s a r e : ) ; 4 disp ( ( E 3) ux+vx=x . . . . . ( i ) 3 ux+(E 5) vx=4x . . . . . . ( i i 5 6 7 8 9 10 11 12 13 14
) ); disp ( S i m p l i f y i n g we g e t (E2 8E+12) ux=14x 4x ) ; syms c1 c2 c3 x disp ( C u m u l a t i v e f u n c t i o n i s g i v e n by E2 8 E+12 =0 ) ; E = poly (0 , E ) ; f = E ^2 -8* E +12; r = roots ( f ) ; disp ( r ) ; disp ( There f o r t h e c o m p l e t e s o l u t i o n i s = c f + p i ) ; cf = c1 * r (1) ^ x + c2 * r (2) ^ x ; disp ( CF= ) ; 139
15 disp ( cf ) ; 16 // p a r t i c u l a r i n t e g r a l c a l u l a t i o n m a n u a l l y 17 disp ( s o l v i n g f o r PI ) ; 18 disp ( We g e t PI= ) ; 19 pi = -4/5* x -19/25+4^ x /4; 20 ux = cf + pi ; 21 disp ( ux= ) ; 22 disp ( ux ) ; 23 disp ( P u t t i n g i n ( i ) we g e t vx= ) ; 24 vx = c1 *2^ x -3* c2 *6^ x -3/5* x -34/25 -4^ x /4; 25 disp ( vx ) ;
Scilab code Exa 26.15.2 Z transform

1 // q u e s 1 5 ( i i ) 2 syms n z 3 y = z ^( - n ) ; 4 f = symsum (y ,n ,0 , %inf ) ; 5 disp ( f ) ;
Scilab code Exa 26.16 evaluating u2 and u3

1 2 3 4 5 6 7 8 9 10 11
// q u e s 1 6 syms z // f =(2/ z 2+5/ z 3+14/ z 4 ) /(1 1/ z ) 4 f =(2/ z ^2+5/ z +14) /(1/ z -1) ^4 u0 = limit (f ,z ,0) ; u1 = limit (1/ z *( f - u0 ) ,z ,0) ; u2 = limit (1/ z ^2*( f - u0 - u1 * z ) ,z ,0) ; disp ( u2= ) ; disp ( u2 ) ; u3 = limit (1/ z ^3*( f - u0 - u1 *z - u2 * z ^2) ,z ,0) ; disp ( u3= ) ; 140
12
disp ( u3 ) ;
141
Chapter 27 Numerical Solution of Ordinary Dierential Equations
Scilab code Exa 27.1 solving ODE with picards method

1 2 3 4 5 6 7 8 9 10 11 12
// q u e s 1 syms x disp ( s o l u t i o n t h r o u g h p i c a r d s method ) ; n = input ( The no o f i t e r a t i o n s r e q u i r e d ) ; disp ( y ( 0 ) =1 and y ( x )=x+y ) ; yo =1; yn =1; for i = 1: n yn = yo + integ ( yn +x ,x ,0 , x ) ; end disp ( y= ) ; disp ( yn ) ;
Scilab code Exa 27.2 solving ODE with picards method

1
// e r r o r 142
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// q u e s 2 syms x disp ( s o l u t i o n t h r o u g h p i c a r d s method ) ; n = input ( The no o f i t e r a t i o n s r e q u i r e d ) ; disp ( y ( 0 ) =1 and y ( x )=x+y ) ; yo =1; y =1; for i = 1: n f =( y - x ) /( y + x ) ; y = yo + integ (f ,x ,0 , x ) ; end disp ( y= ) ; x =0.1; disp ( eval ( y ) ) ;
Scilab code Exa 27.5 solving ODE using Eulers method

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
// q u e s 5 clc disp ( S o l u t i o n u s i n g E u l e r s Method ) ; disp x y ; n = input ( I n p u t t h e number o f i t e r a t i o n : ) ; x =0; y =1; for i =1: n y1 = x + y ; y = y +0.1* y1 ; x = x +0.1; end disp ( The v a l u e o f y i s : ) ; disp ( y ) ;
143
Scilab code Exa 27.6 solving ODE using Eulers method

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
// q u e s 5 clc disp ( S o l u t i o n u s i n g E u l e r s Method ) ; disp x y ; n = input ( I n p u t t h e number o f i t e r a t i o n : ) ; x =0; y =1; for i =1: n y1 =( y - x ) /( y + x ) ; y = y +0.02* y1 ; x = x +0.1; disp ( y ) ; end disp ( The v a l u e o f y i s : ) ; disp ( y ) ;
Scilab code Exa 27.7 solving ODE using Modied Eulers method
1 2 3 4 5 6 7 8 9 10 11
// q u e s 7 clc disp ( S o l u t i o n u s i n g E u l e r s Method ) ; disp x y ; n = input ( I n p u t t h e number o f i t e r a t i o n : ) ; x =0.1; m =1; y =1; yn =1; y1 =1; k =1; 144
12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
for i =1: n yn = y ;
for i =1:4 m =( k + y1 ) /2; yn = y +0.1* m ; y1 =( yn + x ) ; disp ( yn ) ; end disp ( ) ; y = yn ; m = y1 ; yn = yn +0.1* m ; disp ( yn ) ; x = x +0.1; yn = y ; k=m; end disp ( The v a l u e o f y i s : ) ; disp ( y ) ;
1 2 3 4 5 6 7 8 9 10
// q u e s 7 clc disp ( S o l u t i o n u s i n g E u l e r s Method ) ; disp x y ; n = input ( I n p u t t h e number o f i t e r a t i o n : ) ; x =0.2; m =0.301; y =2; yn =2; y1 = log10 (2) ; 145
11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
k =0.301; for i =1: n yn = y ;
for i =1:4 m =( k + y1 ) /2; yn = y +0.2* m ; y1 = log10 ( yn + x ) ; disp ( yn ) ; end disp ( ) ; y = yn ; m = y1 ; yn = yn +0.2* m ; disp ( yn ) ; x = x +0.2; yn = y ; k=m; end disp ( The v a l u e o f y i s : ) ; disp ( y ) ;
1 2 3 4 5 6 7 8 9
// q u e s 7 clc disp ( S o l u t i o n u s i n g E u l e r s Method ) ; disp x y ; n = input ( I n p u t t h e number o f i t e r a t i o n : ) ; x =0.2; m =1; y =1; yn =1; 146
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
y1 =1; k =1; for i =1: n yn = y ;
for i =1:4 m =( k + y1 ) /2; yn = y +0.2* m ; y1 =( sqrt ( yn ) + x ) ; disp ( yn ) ; end disp ( ) ; y = yn ; m = y1 ; yn = yn +0.2* m ; disp ( yn ) ; x = x +0.2; yn = y ; k=m; end disp ( The v a l u e o f y i s : ) ; disp ( y ) ;
Scilab code Exa 27.10 solving ODE using runge method

1 // q u e s 1 0 2 disp ( Runges method ) ; 3 function y = f (x , y ) 4 y=x+y; 5 endfunction 6 7 x =0; 8 y =1;
147
9 10 11 12 13 14 15 16 17
h =0.2; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; kk = h * f ( x +h , y + k1 ) ; k3 = h * f ( x +h , y + kk ) ; k =1/6*( k1 +4* k2 + k3 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ;
Scilab code Exa 27.11 solving ODE using runge kutta method
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
// q u e s 1 1 disp ( Runga k u t t a method ) ; function y = f (x , y ) y=x+y; endfunction x =0; y =1; h =0.2; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ;
1
// q u e s 1 2 148
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
clc disp ( Runga k u t t a method ) ; function y = f (x , y ) y =( y ^2 - x ^2) /( x ^2+ y ^2) ; endfunction x =0; y =1; h =0.2; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ; disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; x =0.2; h =0.2; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ;
1 // q u e s 1 2 2 clc 3 disp ( Runga k u t t a method ) ;
149
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
function yy = f (x , y ) yy = x + y ^2; endfunction x =0; y =1; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ; disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; x =0.1; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y=y+k; disp ( y ) ;
Scilab code Exa 27.14 solving ODE using milnes method

1 // q u e s 1 4 2 clc 3 syms x 4 yo =0; 5 y =0;
150
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
h =0.2; f =x - y ^2; y = integ (f ,x ,0 , x ) ; y1 = eval ( yo + y ) ; disp ( y1= ) ; disp ( float ( y1 ) ) ; f =x - y ^2; y = integ (f ,x ,0 , x ) ; y2 = yo + y ; disp ( y2= ) ; disp ( float ( y2 ) ) ; // f u n c t i o n y=f ( x , y ) y =x - y ^2; // e n d f u n c t i o n y = integ (f ,x ,0 , x ) ; y3 = yo + y ; disp ( y3= ) ; disp ( float ( y3 ) ) ; disp ( d e t e r m i n i n g t h e i n i t i a l v a l u e s f o r m i l n e s method u s i n g y3 ) ; disp ( x =0.0 y0 =0.0 f 0 =0 ) ; disp ( x =0.2 y1= ) ; x =0.2; disp ( eval ( y1 ) ) ; y1 = eval ( y1 ) ; disp ( f 1= ) ; f1 = float ( eval (x - y1 ^2) ) ; disp ( f1 ) ; disp ( x =0.4 y2= ) ; x =0.4; disp ( float ( eval ( y2 ) ) ) ; disp ( f 2= ) ; f2 = float ( eval (x - y2 ^2) ) ; disp ( f2 ) ; y3= ) ;
26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 disp ( x =0.6 42 x =0.6;
151
43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
disp ( eval ( y3 ) ) ; disp ( f 3= ) ; f3 = float ( eval (x - y3 ^2) ) ; disp ( f3 ) ; // disp ( U s i n g p r e d i c t o r method t o f i n d y4 ) ; x =0.8; y4 = eval ( yo +4/3* h *(2* f1 - f2 +2* f3 ) ) ; disp ( y4= ) ; disp ( float ( y4 ) ) ; f4 = float ( eval (x - y ^2) ) ; disp ( f 4= ) ; disp ( f4 ) ; disp ( U s i n g p r e d i c t o r method t o f i n d y5 ) ; x =1.0; y5 = eval ( y1 +4/3* h *(2* f2 - f3 +2* f4 ) ) ; disp ( float ( y5 ) ) ; f5 = float ( eval (x - y ^2) ) ; disp ( f 5= ) ; disp ( f5 ) ; disp ( Hence y ( 1 )= ) ; disp ( float ( y5 ) ) ;
Scilab code Exa 27.15 solving ODE using runge kutta and milnes method
1 2 3 4 5 6 7 8 9 10
// q u e s 1 5 clc disp ( Runga k u t t a method ) ; function yy = f (x , y ) yy = x * y + y ^2; endfunction y0 =1; x =0; y =1; 152
11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; ka =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y1 = y + ka ; y = y + ka ; disp ( y ) ; // x = 0 . 1 ; // y1= f l o a t ( e v a l ( y ) ) ; disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; x =0.1; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; kb =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y2 = y + kb ; y = y + kb ; disp ( y ) ; // x = 0 . 2 ; // y2= f l o a t ( e v a l ( y ) ) ; disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; x =0.2; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; kc =1/6*( k1 +2* k2 +2* k3 + k4 ) ; 153
47 48 49 50 51 52 53 54 55 56 57 58 59 60 61
62 63 64 65 66 67 68 69 disp ( u s i n g c o r r e c t o r method : ) ; 70 y4 = y2 + h /3*( f2 +4* f3 + f4 ) ; 71 disp ( y4= ) ; 72 disp ( y4 ) ; 73 disp ( f 4= ) ; 74 f4 = f (0.4 , y4 ) ; 75 disp ( f4 ) ;
disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y3 = y + kc ; y = y + kc ; disp ( y ) ; // x = 0 . 3 ; // y3= f l o a t ( e v a l ( y ) ) ; f0 = f (0 , y0 ) ; f1 = f (0.1 , y1 ) ; f2 = f (0.2 , y2 ) ; f3 = f (0.3 , y3 ) ; disp ( y0 y1 y2 y3 a r e r e s p e c t i v e l y : ) ; disp ( y3 , y2 , y1 , y0 ) ; disp ( f 0 f 1 f 2 f 3 a r e r e s p e c t i v e l y : ) ; disp ( f3 ,f2 , f1 , f0 ) ; disp ( f i n d i n g y4 u s i n g p r e d i c t o r s m i l n e method x =0.4 ); h =0.1; y4 = y0 +4* h /3*(2* f1 - f2 +2* f3 ) ; disp ( y4= ) ; disp ( y4 ) ; disp ( f 4= ) ; f4 = f (0.4 , y4 ) ;
Scilab code Exa 27.16 solving ODE using adamsbashforth method

1 // q u e s 1 6 2 clc
154
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
function yy = f (x , y ) yy = x ^2*(1+ y ) ; endfunction y3 =1 y2 =1.233 y1 =1.548 y0 =1.979 f3 = f (1 , y3 ) f2 = f (1.1 , y2 ) f1 = f (1.2 , y1 ) f0 = f (1.3 , y0 ) disp ( u s i n g p r e d i c t o r method ) ; h =0.1 y11 = y0 + h /24*(55* f0 -59* f1 +37* f2 -9* f3 ) disp ( y11= ) ; disp ( y11 ) ; x =1.4; f11 = f (1.4 , y11 ) ; disp ( u s i n g c o r r e c t o r method ) ; y11 = y0 + h /24*(9* f11 +19* f0 -5* f1 + f2 ) ; disp ( y11= ) ; disp ( y11 ) ; f11 = f (1.4 , y11 ) ; disp ( f 1 1= ) ; disp ( f11 ) ;
Scilab code Exa 27.17 solving ODE using runge kutta and adams method
1 // q u e s 1 7 2 clc 3 disp ( Runga k u t t a method ) ; 4 5 function yy = f (x , y )
155
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
yy =x - y ^2; endfunction y0 =1; x =0; y =1; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; ka =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y1 = y + ka ; y = y + ka ; disp ( y ) ; // x = 0 . 1 ; // y1= f l o a t ( e v a l ( y ) ) ; disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; x =0.1; h =0.1; k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; kb =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y2 = y + kb ; y = y + kb ; disp ( y ) ; // x = 0 . 2 ; // y2= f l o a t ( e v a l ( y ) ) ;
disp ( t o f i n d y ( 0 . 4 ) put x =0.2 y=a b o v e v a l u e i e 1 . 1 9 6 h =0.2 ) ; 40 x =0.2; 41 h =0.1; 156
42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75
k1 = h * f (x , y ) ; k2 = h * f ( x +1/2* h , y +1/2* k1 ) ; k3 = h * f ( x +1/2* h , y +1/2* k2 ) ; k4 = h * f ( x +h , y + k3 ) ; kc =1/6*( k1 +2* k2 +2* k3 + k4 ) ; disp ( t h e r e q u i r e d a p p r o x i m a t e v a l u e i s : ) ; y3 = y + kc ; y = y + kc ; disp ( y ) ; // x = 0 . 3 ; // y3= f l o a t ( e v a l ( y ) ) ; f0 = f (0 , y0 ) ; f1 = f (0.1 , y1 ) ; f2 = f (0.2 , y2 ) ; f3 = f (0.3 , y3 ) ; disp ( y0 y1 y2 y3 a r e r e s p e c t i v e l y : ) ; disp ( y3 , y2 , y1 , y0 ) ; disp ( f 0 f 1 f 2 f 3 a r e r e s p e c t i v e l y : ) ; disp ( f3 ,f2 , f1 , f0 ) ; disp ( U s i n g adams method ) ; disp ( U s i n g t h e p r e d i c t o r ) ; h =0.1; y4 = y3 + h /24*(55* f3 -59* f2 +37* f1 -9* f0 ) ; x =0.4; f4 = f (0.4 , y4 ) ; disp ( y4= ) ; disp ( y4 ) ; disp ( u s i n g c o r r e c t o r method ) ; y4 = y3 + h /24*(9* f4 +19* f3 -5* f2 + f1 ) ; disp ( y4= ) ; disp ( y4 ) ; f4 = f (0.4 , y4 ) ; disp ( f 4= ) ; disp ( f4 ) ;
157
Scilab code Exa 27.18 solving simultanious ODE using picards method
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
// q u e s 1 8 clc disp ( P i c a r d s method ) ; x0 =0; y0 =2; z0 =1; syms x function yy = f (x ,y , z ) yy = x + z ; endfunction function yy = g (x ,y , z ) yy =x - y ^2; endfunction disp ( f i r s t a p p r o x i m a t i o n ) ; y1 = y0 + integ ( f (x , y0 , z0 ) ,x , x0 , x ) ; disp ( y1= ) ; disp ( y1 ) ; z1 = z0 + integ ( g (x , y0 , z0 ) ,x , x0 , x ) ; disp ( z 1= ) ; disp ( z1 ) ; disp ( s e c o n d a p p r o x i m a t i o n ) ; y2 = y0 + integ ( f (x , y1 , z1 ) ,x , x0 , x ) ; disp ( y2= ) ; disp ( y2 ) ; z2 = z0 + integ ( g (x , y1 , z1 ) ,x , x0 , x ) ; disp ( z 2= ) ; disp ( z2 ) ; disp ( t h i r d a p p r o x i m a t i o n ) ; y3 = y0 + integ ( f (x , y2 , z2 ) ,x , x0 , x ) ; disp ( y3= ) ; disp ( y3 ) ; z3 = z0 + integ ( g (x , y2 , z2 ) ,x , x0 , x ) ; disp ( z 3= ) ; 158
37 disp ( z3 ) ; 38 x =0.1; 39 disp ( y ( 0 . 1 ) = ) ; 40 disp ( float ( eval ( y3 ) ) ) ; 41 disp ( z ( 0 . 1 ) = ) ; 42 disp ( float ( eval ( z3 ) ) ) ;
Scilab code Exa 27.19 solving ssecond ODE using runge kutta method
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
// q u e s 1 9 clc syms x function yy = f (x ,y , z ) yy = z ; endfunction function yy = g (x ,y , z ) yy = x * y ^2 - y ^2; endfunction x0 =0; y0 =1; z0 =0; h =0.2; disp ( u s i n g k1 k2 . . f o r f and l 1 l 2 . . . f o r g r u n g a k u t t a f o r m u l a e becomes ) ; h =0.2; k1 = h * f ( x0 , y0 , z0 ) ; l1 = h * g ( x0 , y0 , z0 ) ; k2 = h * f ( x0 +1/2* h , y0 +1/2* k1 , z0 +1/2* l1 ) ; l2 = h * g ( x0 +1/2* h , y0 +1/2* k1 , z0 +1/2* l1 ) ; k3 = h * f ( x0 +1/2* h , y0 +1/2* k2 , z0 +1/2* l2 ) ; l3 = h * g ( x0 +1/2* h , y0 +1/2* k2 , z0 +1/2* l2 ) ; k4 = h * f ( x0 +h , y0 + k3 , z0 + l3 ) ; l4 = h * g ( x0 +h , y0 + k3 , z0 + l3 ) ; k =1/6*( k1 +2* k2 +2* k3 + k4 ) ; l =1/6*( l1 +2* l2 +2* l3 +2* l4 ) ; 159
26 27 28 29 30 31 32 33 34 35
// a t x =0.2 x =0.2; y = y0 + k ; y1 = z0 + l ; disp ( y= ) ; disp ( float ( y ) ) ; disp ( y1= ) ; disp ( float ( y1 ) ) ; y
Scilab code Exa 27.20 solving ODE using milnes method

1 // q u e s 2 0 2 clc
160
Chapter 28 Numerical Solution of Partial Dierential Equations
Scilab code Exa 28.1 classication of partial dierential equation

1 // q u e s 2 8 . 1 2 clear 3 clc 4 disp ( D=B2 4AC ) ; 5 disp ( i f D<0 t h e n e l l i p t i c 6 7 8 9 10 11 12
i f D=0 t h e n p a r a b o l i c s o The
i f D>0 t h e n h y p e r b o i c ) ; disp ( ( i ) A=x 2 , B1y 2 D=42 4 1 4=0 e q u a t i o n i s PARABOLIC ) ; disp ( ( i i ) D=4x 2 ( y 2 1) ) ; disp ( f o r i n f <x< i n f and 1<y <1 D<0 ) ; disp ( So t h e e q u a t i o n i s ELLIPTIC ) ; disp ( ( i i i ) A=1+x 2 ,B=5+2x 2 ,C=4+x 2 ) ; disp ( D=9>0 ) ; disp ( So t h e e q u a t i o n i s HYPERBOLIC ) ;
Scilab code Exa 28.2 solving elliptical equation 161
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
// q u e s 2 8 . 2 disp ( S e e f i g u r e i n q u e s t i o n ) ; disp ( From symmetry u7=u1 , u8=u2 , u9=u3 , u3=u1 , u6=u4 , u9=u7 ) ; disp ( u5 =1/4 (2000+2000+1000+1000) =1500 ) ; u5 =1500; disp ( u1 =1/4(0=1500+1000+2000) =1125 ) ; u1 =1125; disp ( u2 =1/4 (1125+1125+1000+1500) =1188 ) ; u2 =1188; disp ( u4 =1/4(2000+1500+1125+1125) =1438 ) ; u4 =1438; disp ( u1 , u2 , u4 , u5 ) disp ( I t e r a t i o n s : ) ; // n=i n p u t ( I n p u t t h e number o f i t e r a t i o n s r e q u i r e d : ) ; for i =1:6 u11 =1/4*(1000+ u2 +500+ u4 ) ; u22 =1/4*( u11 + u1 +1000+ u5 ) ; u44 =1/4*(2000+ u5 + u11 + u1 ) ; u55 =1/4*( u44 + u4 + u22 + u2 ) ; disp ( ) ; disp ( u55 , u44 , u22 , u11 ) ; u1 = u11 ; u2 = u22 ; u4 = u44 ; u5 = u55 ; end
Scilab code Exa 28.3 evaluating function satisfying laplace equation

1 // q u e s 3 2 clear 3 clc 4 disp ( S e e f i g u r e
i n q u e s t i o n ); 162
5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
disp ( To f i n d t h e i n i t i a l v a l u e s o f u1 u2 u3 u4 we assume u4=0 ) ; disp ( u1 =1/4 (1000+0+1000+2000) =1000 ) ; u1 =1000; disp ( u2 =1/4(1000+500+1000+500) =625 ) ; u2 =625; disp ( u3 =1/4 (2000+0+1000+500) =875 ) ; u3 =875; disp ( u4 =1/4(875+0+625+0) =375 ) ; u4 =375; disp ( u1 , u2 , u3 , u4 ) disp ( I t e r a t i o n s : ) ; // n=i n p u t ( I n p u t t h e number o f i t e r a t i o n s r e q u i r e d : ) ; for i =1:6 u11 =1/4*(2000+ u2 +1000+ u3 ) ; u22 =1/4*( u11 +500+1000+ u4 ) ; u33 =1/4*(2000+ u4 + u11 +500) ; u44 =1/4*( u33 +0+ u22 +0) ; disp ( ) ; disp ( u44 , u33 , u22 , u11 ) ; u1 = u11 ; u2 = u22 ; u4 = u44 ; u3 = u33 ; end
Scilab code Exa 28.4 solution of poissons equation

1 // q u e s 4 2 clear 3 clc 4 disp ( S e e f i g u r e i n q u e s t i o n ) ; 5 disp ( u s i n g n u m e r i c a l p o i s s o n s e q u a t i o n u ( i 1) ( j )+u (
i +1) ( j )+u ( i ) ( j 1)+u ( i ) ( j +1)=h 2 f ( i h , j h ) ) ; 163
6 7 8 9 10 11 12 13 14 15 16 17
18 19 20 21 22 23 // c r e a t e u p p e r t r i a n g u l a r m a t r i x 24 s =0; 25 for j =1: n -1 26 if A (j , j ) ==0 27 k=j; 28 for k = k +1: n 29 if A (k , j ) ==0 30 continue 31 end 32 break 33 end 34 B = A (j ,:) ; C = r ( j ) ; 35 A (j ,:) = A (k ,:) ; r ( j ) = r ( k ) ; 36 A (k ,:) = B ; r ( k ) = C ; 37 end 38 for i =1+ s :n -1 39 L = A ( i +1 , j ) / A (j , j ) ; 40 A ( i +1 ,:) = A ( i +1 ,:) -L * A (j ,:) ; 41 r ( i +1) = r ( i +1) -L * r ( j ) ;
disp ( Here f ( x , y ) = 10(x2+y 2+10 ) ; disp ( Here f o r u1 i =1 , j =2 p u t t i n g i n e q u a t i o n t h i s g i v e s : ); disp ( u1 =1/4( u2+u3 +150 ) ; disp ( s i m i l a r l y ) ; disp ( u2 =1/4( u1+u4 +180 ) ; disp ( u3 =1/4( u1+u4 +120 ) ; disp ( u4 =1/4( u2+u3 +150 ) ; disp ( r e d u c i n g t h e r s e e q u a t i o n s s i n c e u4=u1 ) ; disp ( 4 u1u2u3 150=0 ) ; disp ( u1 2u2+90=0 ) ; disp ( u1 2u3+60=0 ) ; disp ( S o l v n g t h e s e e q u a t i o n s by Gauss j o r d o n method ); A =[4 -1 -1;1 -2 0;1 0 -2]; r =[150; -90; -60]; D=A;d=r; n =3;
164
42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58
end s = s +1; end // S o l u t i o n o f e q u a t i o n s x ( n ) = r ( n ) / A (n , n ) ; for i =n -1: -1:1 sum =0; for j = i +1: n sum = sum + A (i , j ) * x ( j ) ; end x ( i ) =(1/ A (i , i ) ) *( r ( i ) - sum ) ; end // h e c k i n g w i t h s c i l a b f u n c t i o n s p = inv ( D ) * d ; // Output disp ( @ @) disp ( Output [ B ] [ x ] = [ b ] ) disp ( Upper r i a n g u l a r M a t r i x [ B ] = ) ; disp ( A ) disp ( M a t r i x [ b ] = ) ; disp ( r ) disp ( s o l u t i o n o f l i n e a r e q u a t i o n s : ) ; disp (x )
59 60 61 62
Scilab code Exa 28.5 solving parabolic equation

1 // q u e s 5 2 clear 3 clc 4 disp ( Here c 2=4 , h=1 , k=1/8 , 5
t h e r e f o r e a l p h a =( c
2) k /( h 2) ); disp ( U s i n g b e n d r e s c h m i d i t s r e c u r r e n c e r e l a t i o n i e u ( i ) ( j +1)=t u ( i 1) ( j )+t u ( i +1) ( j ) +(1 2 t ) u ( i , j ) ) ; 6 disp ( Now s i n c e u ( 0 , t )=0=u ( 8 , t ) t h e r e f o r e u ( 0 , i ) =0 165
7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
and u ( 8 , j ) =0 and u ( x , 0 ) =4x 1/2x 2 ) ; c =2; h =1; k =1/8; t =( c ^2) * k /( h ^2) ; A = ones (9 ,9) ; for i =1:9 for j =1:9 A (1 , i ) =0; A (9 , i ) =0; A (i ,1) =4*( i -1) -1/2*( i -1) ^2; end end // i =2; // j =2; for i =2:8 for j =2:7 // A( i , j ) =1/2 (A( i 1 , j 1)+A( i +1 , j 1) ) ; A (i , j ) = t * A (i -1 ,j -1) + t * A ( i +1 ,j -1) +(1 -2* t ) * A (i -1 ,j -1) ; end end for i =2:8 j =2; disp ( A (i , j ) ) ;
27 28 29 30 31 32 33 end
Scilab code Exa 28.6 solving heat equation

1 // q u e s 5 2 clear 3 clc
166
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
disp ( Here c 2=1 , h=1/3 , k =1/36 , t h e r e f o r e t =( c 2 ) k / ( h 2 ) =1/4 ) ; disp ( So b e n d r e s c h m i d i t s r e c u r r e n c e r e l a t i o n i e u ( i ) ( j +1) =1/4( u ( i 1) ( j )+u ( i +1) ( j ) +2u ( i , j ) ) ; disp ( Now s i n c e u ( 0 , t )=0=u ( 1 , t ) t h e r e f o r e u ( 0 , i ) =0 and u ( 1 , j ) =0 and u ( x , 0 ) =s i n ( %pi ) x ) ; c =1; h =1/3; k =1/36; t =( c ^2) * k /( h ^2) ; A = ones (9 ,9) ; for i =1:9 for j =1:9 A (1 , i ) =0; A (2 , i ) =0; A (i ,1) = sin ( %pi /3*( i -1) ) ; end end //A( 2 , 1 ) = 0 . 8 6 6 ; //A( 3 , 1 ) = 0 . 8 6 6 ; for i =2:8 for j =2:8 // A( i , j ) =1/4 (A( i 1 , j 1)+A( i +1 , j 1)+2A( i 1 , j 1) ) ; A (i , j ) = t * A (i -1 ,j -1) + t * A ( i +1 ,j -1) +(1 -2* t ) * A (i -1 , j -1) ; end end for i =2:8 j =2; disp ( A (i , j ) ) ;
27 28 29 30 31 32 33 end
167
Scilab code Exa 28.7 solving wave equation

1 // q u e s 7 2 clear 3 clc 4 disp ( Here c 2=16 , t a k i n g h=1 , f i n d i n g k s u c h t h a t 5
6 7 8 9 10 11 12 13
14 15 16 17 18 19 20 21 end 22 end 23 for i =1:4 24 A ( i +1 ,2) =1/2*( A (i ,1) + A ( i +2 ,1) ) ; 25 26 end
c 2 t 2=1 ) ; disp ( So b e n d r e s c h m i d i t s r e c u r r e n c e r e l a t i o n i e u ( i ) ( j +1) =(16 t 2 ( u ( i 1) ( j )+u ( i +1) ( j ) ) +2(1 16 t 2 u ( i , j )u ( i ) ( j 1) ) ; disp ( Now s i n c e u ( 0 , t )=0=u ( 5 , t ) t h e r e f o r e u ( 0 , i ) =0 and u ( 5 , j ) =0 and u ( x , 0 ) =x 2(5 x ) ) ; c =4; h =1; k =( h / c ) ; t=k/h; A = zeros (6 ,6) ; disp ( A l s o from 1 s t d e r i v a t i v e ( u ( i ) ( j +1)u ( i , j 1) ) /2 k=g ( x ) and g ( x ) =0 i n t h i s c a s e ) ; disp ( So i f j =0 t h i s g i v e s u ( i ) ( 1 ) =1/2 ( u ( i 1) ( 0 )+u ( i +1) ( 0 ) ) ) for i =0:5 for j =2:9 A (1 , i +1) =0; A (6 , i +1) =0; A ( i +1 ,1) =( i ) ^2*(5 - i ) ;
168
27 28 29 30 31 32 33 34 35 36 37 38
for i =3:5 for j =3:5 A (i -1 , j ) =( c * t ) ^2*( A (i -2 ,j -1) + A (i ,j -1) ) +2*(1 -( c * t ) ^2) * A (i -1 ,j -1) -A (i -1 ,j -2) ; end end for i =1:5 for j =1:5 disp ( A (i , j ) ) ; end end
Scilab code Exa 28.8 solving wave equation

1 // q u e s 8 2 clear 3 clc 4 disp ( Here c 2=4 , t a k i n g h=1 , f i n d i n g k s u c h t h a t 5
6 7 8 9 10 11 12 13
c 2 t 2=1 ) ; disp ( So b e n d r e s c h m i d i t s r e c u r r e n c e r e l a t i o n i e u ( i ) ( j +1) =(16 t 2 ( u ( i 1) ( j )+u ( i +1) ( j ) ) +2(1 16 t 2 u ( i , j )u ( i ) ( j 1) ) ; disp ( Now s i n c e u ( 0 , t )=0=u ( 4 , t ) t h e r e f o r e u ( 0 , i ) =0 and u ( 4 , j ) =0 and u ( x , 0 ) =x(4 x ) ) ; c =2; h =1; k =( h / c ) ; t=k/h; A = zeros (6 ,6) ; disp ( A l s o from 1 s t d e r i v a t i v e ( u ( i ) ( j +1)u ( i , j 1) ) /2 k=g ( x ) and g ( x ) =0 i n t h i s c a s e ) ; disp ( So i f j =0 t h i s g i v e s u ( i ) ( 1 ) =1/2 ( u ( i 1) ( 0 )+u ( i +1) ( 0 ) ) ) 169
14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
for i =0:5 for j =2:9 A (1 , i +1) =0; A (5 , i +1) =0; A ( i +1 ,1) =( i ) *(4 - i ) ;
end end for i =1:4 A ( i +1 ,2) =1/2*( A (i ,1) + A ( i +2 ,1) ) ; end for i =3:5 for j =3:5 A (i -1 , j ) =( c * t ) ^2*( A (i -2 ,j -1) + A (i ,j -1) ) +2*(1 -( c * t ) ^2) * A (i -1 ,j -1) -A (i -1 ,j -2) ; end end for i =1:5 for j =1:5 disp ( A (i , j ) ) ; end end
170
Chapter 34 Probability and Distributions
Scilab code Exa 34.1 Calculating probability

1 clear 2 clc 3 disp ( from t h e
p r i n c i p l e o f c o u n t i n g , t h e r e q u i r e d no . o f ways a r e 12 11 10 9= ) 4 12*11*10*9
Scilab code Exa 34.2.1 Calculating the number of permutations

1 clear 2 clc 3 disp ( no . o f p e r m u t a t i o n s = 9 ! / ( 2 ! 2 ! 2 ! ) ) 4 factorial (9) /( factorial (2) * factorial (2) * factorial (2)
Scilab code Exa 34.2.2 Number of permutations
171
1 clear 2 clc 3 disp ( no . o f p e r m u t a t i o n s = 9 ! / ( 2 ! 2 ! 3 ! 3 ! ) ) 4 factorial (9) /( factorial (2) * factorial (2) * factorial (3)
* factorial (3) )
Scilab code Exa 34.3.1 Calculating the number of committees

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( no . o f c o m m i t t e e s=C ( 6 , 3 ) C ( 5 , 2 )= ) 7 C (6 ,3) * C (5 ,2)
Scilab code Exa 34.3.2 Finding the number of committees

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( no . o f c o m m i t t e e s=C ( 4 , 1 ) C ( 5 , 2 )= ) 7 C (4 ,1) * C (5 ,2)
Scilab code Exa 34.3.3 Finding the number of committees

1 clear 2 clc
172
3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( no . o f c o m m i t t e e s=C ( 6 , 3 ) C ( 4 , 2 )= ) 7 C (6 ,3) * C (4 ,2)
Scilab code Exa 34.4.1 Finding the probability of getting a four in a single throw of a die
1 clear 2 clc 3 disp ( t h e 4 1/6
p r o b a b i l i t y o f g e t t i n g a f o u r i s 1/6= )
Scilab code Exa 34.4.2 Finding the probability of getting an even number in a single throw of a die
p r o b a b i l i t y o f g e t t i n g an e v e n no . 1/2= )
Scilab code Exa 34.5 Finding the probability of 53 sundays in a leap year
p r o b a b i l i t y o f 53 s u n d a y s i s 2/7= )
173
Scilab code Exa 34.6 probability of getting a number divisible by 4 under given conditions
1 clear 2 clc 3 disp ( t h e f i v e d i g i t s can be a r r a n g e d i n 5 ! ways = ) 4 factorial (5) 5 disp ( o f which 4 ! w i l l b e g i n w i t h 0= ) 6 factorial (4) 7 disp ( so , t o t a l no . o f f i v e d i g i t numbers =5! 4!= ) 8 factorial (5) - factorial (4) 9 disp ( t h e numbers e n d i n g i n 0 4 , 1 2 , 2 0 , 2 4 , 3 2 , 4 0 w i l l 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
be d i v i s i b l e by 4 ) disp ( numbers e n d i n g i n 04=3! ) factorial (3) disp ( numbers e n d i n g i n 12=3! 2! factorial (3) - factorial (2) disp ( numbers e n d i n g i n 20=3! ) factorial (3) disp ( numbers e n d i n g i n 24=3! 2! factorial (3) - factorial (2) disp ( numbers e n d i n g i n 32=3! 2! factorial (3) - factorial (2) disp ( numbers e n d i n g i n 40=3! ) factorial (3) disp ( so , t o t a l no . o f f a v o u r a b l e 6+4+6+4+4+6 disp ( p r o b a b i l i t y =30/96= ) 30/96
) )
ways=6+4+6+4+4+6= )
Scilab code Exa 34.7 Finding the probability
174
1 2 3 4 5 6 7 8 9 10 11
clear clc function [ x ]= C (a , b ) x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) endfunction disp ( t o t a l no . o f p o s s i b l e c a s e s=C ( 4 0 , 4 ) ) C (40 ,4) disp ( f a v o u r a b l e o u t c o m e s=C ( 2 4 , 2 ) C ( 1 5 , 1 )= ) C (24 ,2) * C (15 ,1) disp ( p r o b a b i l i t y = ) ( C (24 ,2) * C (15 ,1) ) / C (40 ,4)

1 2 3 4 5 6 7 8 9 10 11
clear clc function [ x ]= C (a , b ) x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) endfunction disp ( t o t a l no . o f p o s s i b l e c a s e s=C ( 4 0 , 4 ) ) C (15 ,8) disp ( f a v o u r a b l e o u t c o m e s=C ( 2 4 , 2 ) C ( 1 5 , 1 )= ) C (5 ,2) * C (10 ,6) disp ( p r o b a b i l i t y = ) ( C (5 ,2) * C (10 ,6) ) / C (15 ,8)
Scilab code Exa 34.9.1 Finding the probability

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction
175
6 disp ( t o t a l no . o f p o s s i b l e c a s e s=C ( 9 , 3 ) ) 7 C (9 ,3) 8 disp ( f a v o u r a b l e o u t c o m e s=C ( 2 , 1 ) C ( 3 , 1 ) C ( 4 , 1 )= ) 9 C (2 ,1) * C (3 ,1) * C (4 ,1) 10 disp ( p r o b a b i l i t y = ) 11 ( C (2 ,1) * C (3 ,1) * C (4 ,1) ) / C (9 ,3)

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( t o t a l no . o f p o s s i b l e c a s e s=C ( 9 , 3 ) ) 7 C (9 ,3) 8 disp ( f a v o u r a b l e o u t c o m e s=C ( 2 , 2 ) C ( 7 , 1 )+C ( 3 , 2 ) C
( 6 , 1 )+C ( 4 , 2 ) C ( 5 , 1 )= ) 9 C (2 ,2) * C (7 ,1) + C (3 ,2) * C (6 ,1) + C (4 ,2) * C (5 ,1) 10 disp ( p r o b a b i l i t y = ) 11 ( C (2 ,2) * C (7 ,1) + C (3 ,2) * C (6 ,1) + C (4 ,2) * C (5 ,1) ) / C (9 ,3)

1 2 3 4 5 6 7 8 9
clear clc function [ x ]= C (a , b ) x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) endfunction disp ( t o t a l no . o f p o s s i b l e c a s e s=C ( 9 , 3 ) ) C (9 ,3) disp ( f a v o u r a b l e o u t c o m e s=C ( 3 , 3 )+C ( 4 , 3 )= ) C (3 ,3) + C (4 ,3) 176
10 11
disp ( p r o b a b i l i t y = ) 5/84
Scilab code Exa 34.13 probability of drawing an ace or spade from pack of 52 cards
1 clear 2 clc 3 disp ( p r o b a b i l i t y
o f d r a w i n g an a c e o r s p a d e o r b o t h from pack o f 52 c a r d s =4/52+13/52 1/52= ) 4 4/52+13/52 -1/52

1 clear 2 clc 3 disp ( p r o b a b i l i t y 4 4/52 5 disp ( p r o b a b i l i t y 6 7
o f f i r s t c a r d b e i n g a k i n g =4/52 ) o f s e c o n d c a r d b e i n g a q u e e n =4/52
) 4/52 disp ( p r o b a b i l i t y o f d r a w i n g b o t h c a r d s i n s u c c e s s i o n =4/52 4/52= ) 8 4/52*4/52

1 clear 2 clc
177
3 4 5 6 7 8
disp ( p r o b a b i l i t y o f g e t t i n g 7 i n f i r s t t o s s and n o t g e t t i n g i t i n s e c o n d t o s s =1/6 5/6 ) 1/6*5/6 disp ( p r o b a b i l i t y o f n o t g e t t i n g 7 i n f i r s t t o s s and g e t t i n g i t i n s e c o n d t o s s =5/6 1/6 ) 5/6*1/6 disp ( r e q u i r e d p r o b a b i l i t y =1/6 5/6+5/6 1/6 ) 1/6*5/6+5/6*1/6

o f not g e t t i n g 7 in e i t h e r t o s s
=5/6 5/6 ) 4 5/6*5/6 5 disp ( p r o b a b i l i t y o f g e t t i n g 7 a t l e a s t o n c e =1 5/6 5/6 ) 6 1 -5/6*5/6

1 clear 2 clc 3 disp ( p r o b a b i l i t y 4 1/6*1/6
o f g e t t i n g 7 t w i c e =1/6 1/6 )

1
clear 178
2 clc 3 disp ( p r o b a b i l i t y
of engineering subject being c h o o o s e n = ( 1 / 3 3 / 8 ) + ( 2 / 3 5 / 8 )= ) 4 (1/3*3/8) +(2/3*5/8)

of white b a l l being choosen =2/6 6/13+4/6 5/13= ) 4 2/6*6/13+4/6*5/13

1 clear 2 clc 3 disp ( c h a n c e s o f w i n n i n g o f A=1/2+(1/2) 2 ( 1 / 2 )
+(1/2) 4 ( 1 / 2 ) +(1/2) 6 ( 1 / 2 ) +..= ) 4 ( 1 / 2 ) /(1 (1/2) 2) 5 d i s p ( c h a n c e s o f w i n n i n g o f B=1 c h a n c e s o f w i n n i n g of A ) 6 1 2/3

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) )
179
5 endfunction 6 disp ( t o t a l no . o f p o s s i b l e o u t c o m e s=C ( 1 0 , 2 )= ) 7 C (10 ,2) 8 disp ( no . o f f a v o u r a b l e o u t c o m e s =55= ) 9 5*5 10 disp ( p= ) 11 25/49

1 clear 2 clc 3 disp ( t o t a l no . o f p o s s i b l e o u t c o m e s =10 9= ) 4 10*9 5 disp ( no . o f f a v o u r a b l e o u t c o m e s =5 5+5 5= ) 6 5*5+5*5 7 disp ( p= ) 8 50/90

1 clear 2 clc 3 disp ( t o t a l no . o f p o s s i b l e o u t c o m e s =10 9= ) 4 10*10 5 disp ( no . o f f a v o u r a b l e o u t c o m e s =5 5+5 5= ) 6 5*5+5*5 7 disp ( p= ) 8 50/100
180

1 2 3 4 5 6 7 8 9 10 11 12 13 14
clear clc A =1/4 B =1/3 AorB =1/2 AandB = A +B - AorB disp ( p r o b a b i l i t y AandB / B disp ( p r o b a b i l i t y AandB / A disp ( p r o b a b i l i t y A - AandB disp ( p r o b a b i l i t y (1/6) /(1 -1/3)
o f A/B=AandB/B= ) o f B/A=AandB/A= ) o f AandBnot=AAandB= ) o f A/ Bnot=AandBnot / Bnot= )

1 clear 2 clc 3 disp ( p r o b a b i l i t y 4 disp ( p r o b a b i l i t y 5 disp ( p r o b a b i l i t y 6 disp ( p r o b a b i l i t y
o f A h i t t i n g t a r g e t =3/5 ) o f B h i t t i n g t a r g e t =2/5 ) o f C h i t t i n g t a r g e t =3/4 ) t h a t two s h o t s h i t =3/5 2/5 (1 3/4) +2/5 3/4 (1 3/5) +3/4 3/5 (1 2/5) ) 7 3/5*2/5*(1 -3/4) +2/5*3/4*(1 -3/5) +3/4*3/5*(1 -2/5)

1 clear 2 clc
181
disp ( p r o b a b i l i t y o f p r o b l e m n o t g e t t i n g s o l v e d =1/2 2/3 3/4= ) 4 1/2*2/3*3/4 5 disp ( p r o b a b i l i t y o f p r o b l e m g e t t i n g s o l v e d =1 (1/2 2/3 3/4)= ) 6 1 -(1/2*2/3*3/4)
3
Scilab code Exa 34.25 nding the probability

1 clc 2 disp ( t o t a l f r e q u e n c y= i n t e g r a t e ( f , x , 0 , 2 )= ) 3 n = integrate ( x 3 , x ,0 ,1) + integrate ( (2 x ) 3 , x 4 5 6 7 8 9 10
11
,1 ,2) disp ( u1 a b o u t o r i g i n = ) u1 =(1/ n ) *( integrate ( ( x ) ( x 3 ) , x ,0 ,1) + integrate ( ( x ) ((2 x ) 3 ) , x ,1 ,2) ) disp ( u2 a b o u t o r i g i n = ) u2 =(1/ n ) *( integrate ( ( x 2 ) ( x 3 ) , x ,0 ,1) + integrate ( ( x 2 ) ((2 x ) 3 ) , x ,1 ,2) ) disp ( s t a n d a r d d e v i a t i o n =(u2u1 2 ) 0.5= ) ( u2 - u1 ^2) ^0.5 disp ( mean d e v i a t i o n a b o u t t h e mean =(1/ n ) ( i n t e g r a t e ( | x 1 | ( x 3 ) , x , 0 , 1 )+ i n t e g r a t e ( | x 1 | ((2 x ) 3 ) , x ,1 ,2 ) ) (1/ n ) *( integrate ( (1 x ) ( x 3 ) , x ,0 ,1) + integrate ( ( x 1) ((2 x ) 3 ) , x ,1 ,2) )

1 clear 2 clc 3 disp ( p r o b a b i l i t y = ( 0 . 4 5 0 . 0 3 )
/(0.450.03+0.250.05+0.30.04= ) 182
(0.45*0.03) /(0.45*0.03+0.25*0.05+0.3*0.04)

1 clear 2 clc 3 disp ( p r o b a b i l i t y = ( 1 / 3 2 / 6 3 / 5 ) 4
/(1/32/63/5+1/31/62/5+1/33/61/5 ) (1/3*2/6*3/5) / (1/ 3*2/ 6*3/ 5+1/ 3*1 /6*2 /5+1 /3*3 /6* 1/5)

1 clc 2 disp ( p r o b a b i l i t y o f no s u c c e s s =8/27 ) 3 disp ( p r o b a b i l i t y o f a s u c c e s s =1/3 ) 4 disp ( p r o b a b i l i t y o f one s u c c e s s =4/9 ) 5 disp ( p r o b a b i l i t y o f two s u c c e s s e s =2/9 ) 6 disp ( p r o b a b i l i t y o f t h r e e s u c c e s s e s =2/9 ) 7 A =[0 1 2 3;8/27 4/9 2/9 1/27] 8 disp ( mean=sum o f i p i= ) 9 A (1 ,1) * A (2 ,1) + A (1 ,2) * A (2 ,2) + A (1 ,4) * A (2 ,4) + A (1 ,3) * A 10 11 12 13
(2 ,3) disp ( sum o f i p i 2= ) A (1 ,1) ^2* A (2 ,1) + A (1 ,2) ^2* A (2 ,2) + A (1 ,4) ^2* A (2 ,4) + A (1 ,3) ^2* A (2 ,3) disp ( v a r i a n c e =(sum o f i p i 2 ) 1= ) A (1 ,1) ^2* A (2 ,1) + A (1 ,2) ^2* A (2 ,2) + A (1 ,4) ^2* A (2 ,4) + A (1 ,3) ^2* A (2 ,3) -1
Scilab code Exa 34.29 nding the probability 183
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
clc syms k A =[0 1 2 3 4 5 6; k 3* k 5* k 7* k 9* k 11* k 13* k ] disp ( sumof a l l p i =1 ) //A( 2 , 1 )+A( 2 , 2 )+A( 2 , 3 ) +(A( 2 , 4 )+A( 2 , 5 )+A( 2 , 6 )+A( 2 , 7 ) disp ( hence , ) k =1/49 disp ( p ( x < 4)= ) a = A (2 ,1) + A (2 ,2) + A (2 ,4) + A (2 ,3) eval ( a ) disp ( eval ( a ) ) disp ( p ( x >=5)= ) b = A (2 ,6) + A (2 ,7) eval ( b ) disp ( eval ( b ) ) disp ( p(3 < x <=6)= ) c = A (2 ,5) + A (2 ,6) + A (2 ,7) eval ( c ) disp ( eval ( c ) ) disp ( p ( x <=2)= ) c = A (2 ,1) + A (2 ,2) + A (2 ,3)

1 clc 2 syms k 3 A =[0 1 2 3 4 5 6 7;0 k 2* k 2* k 3* k k ^2 2* k ^2 7* k ^2+ k 4 5 6 7 8 9 10
] disp ( sumof a l l p i =1 ) //A( 2 , 1 )+A( 2 , 2 )+A( 2 , 3 ) +(A( 2 , 4 )+A( 2 , 5 )+A( 2 , 6 )+A( 2 , 7 ) disp ( hence , ) k =1/10 disp ( p ( x < 6)= ) a = A (2 ,1) + A (2 ,2) + A (2 ,4) + A (2 ,3) + A (2 ,4) + A (2 ,5) + A (2 ,6) eval ( a ) 184
11 disp ( eval ( a ) ) 12 disp ( p ( x >=6)= ) 13 b = A (2 ,7) + A (2 ,8) 14 eval ( b ) 15 disp ( eval ( b ) ) 16 disp ( p(3 < x < 5)= ) 17 c = A (2 ,2) + A (2 ,3) + A (2 ,4) + A (2 ,5) 18 eval ( c ) 19 disp ( eval ( c ) )

1 clc 2 syms x ; 3 f = %e ^( - x ) 4 disp ( c l e a r l y , f >0 f o r 5 6
e v e r y x i n ( 1 , 2 ) and i n t e g r a t e
( f , x , 0 , % i n f )= ) integrate ( %e( y ) , y ,0 , %inf ) disp ( r e q u i r e d p r o b a b i l i t y =p(1<=x <=2)= i n t e g r a t e ( f , x , 1 , 2 )= ) 7 integrate ( %e( y ) , y ,1 ,2) 8 disp ( c u m u l a t i v e p r o b a b i l i t y f u n c t i o n f ( 2 )= i n t e g r a t e ( f , x , %inf , 2 ) = ) 9 integrate ( %e( y ) , y ,0 ,2)

1 2 3 4 5 6
clc syms k ; disp ( t o t a l p = integrate q = integrate r = integrate
p r o b a b i l i t y = i n t e g r a t e ( f , x , 0 , 6 )= ) ( k x , x ,0 ,2) ( 2 k , x ,2 ,4) ( k x+6 k , x ,4 ,6) 185

1 clc 2 A =[ -3 6 9;1/6 1/2 1/3] 3 disp ( f i r s t row o f A d i s p l a y s t h e v a l u e o f x ) 4 disp ( t h e s e c o n d row o f x d i s p l a y s t h e p r o b a b i l i t y
of corresponding to x )
5 disp ( E( x )= ) 6 c = A (1 ,1) * A (2 ,1) + A (1 ,2) *(2 ,2) + A (1 ,3) * A (2 ,3) 7 disp ( E( x ) 2= ) 8 b = A (1 ,1) ^2* A (2 ,1) + A (1 ,2) ^2*(2 ,2) + A (1 ,3) ^2* A (2 ,3) 9 disp ( E( 2 x +1)2=E( 4 x 2+4 x +1) 10 4* b +4* c +1

1 clc 2 disp ( t o t a l f r e q u e n c y= i n t e g r a t e ( f , x , 0 , 2 )= ) 3 n = integrate ( x 3 , x ,0 ,1) + integrate ( (2 x ) 3 , x 4 5 6 7 8 9 10
,1 ,2) disp ( u1 a b o u t o r i g i n = ) u1 =(1/ n ) *( integrate ( ( x ) ( x 3 ) , x ,0 ,1) + integrate ( ( x ) ((2 x ) 3 ) , x ,1 ,2) ) disp ( u2 a b o u t o r i g i n = ) u2 =(1/ n ) *( integrate ( ( x 2 ) ( x 3 ) , x ,0 ,1) + integrate ( ( x 2 ) ((2 x ) 3 ) , x ,1 ,2) ) disp ( s t a n d a r d d e v i a t i o n =(u2u1 2 ) 0.5= ) ( u2 - u1 ^2) ^0.5 disp ( mean d e v i a t i o n a b o u t t h e mean =(1/ n ) ( i n t e g r a t e ( | x 1 | ( x 3 ) , x , 0 , 1 )+ i n t e g r a t e ( | x 1 | ((2 x ) 3 ) , x ,1 ,2 ) ) 186
11 (1/ n ) *( integrate ( (1 x ) ( x 3 ) , x ,0 ,1) + integrate (
( x 1) ((2 x ) 3 ) , x ,1 ,2) )

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( p r o b a b i l i t y t h a t e x a c t l y two w i l l be d e f e c t i v e
=C ( 1 2 , 2 ) ( 0 . 1 ) 2 ( 0 . 9 ) 10= ) l e a s t two w i l l be d e f e c t i v e =1(C ( 1 2 , 0 ) ( 0 . 9 ) 12+C ( 1 2 , 1 ) ( 0 . 1 ) ( 0 . 9 ) 1 1 )= ) 9 1 -( C (12 ,0) *(0.9) ^12+ C (12 ,1) *(0.1) *(0.9) ^11) 10 disp ( t h e p r o b a b i l i t y t h a t none w i l l be d e f e c t i v e =C ( 1 2 , 1 2 ) ( 0 . 9 ) 12= ) 11 C (12 ,12) *(0.9) ^12
7 C (12 ,2) *(0.1) ^2*(0.9) ^10 8 disp ( p r o b a b i l i t y t h a t a t

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( p r o b a b i l i t y o f 8 h e a d s and 4 t a i l s i n 12
t r i a l s =p ( 8 )=C ( 1 2 , 8 ) ( 1 / 2 ) 8 ( 1 / 2 ) 4= )
7 C (12 ,8) *(1/2) ^8*(1/2) ^4 8 disp ( t h e e x p e c t e d no . o f s u c h c a s e s
i n 256 s e t s
=256 p ( 8 ) = ) 187
256*(495/4096)

1 clear 2 clc 3 function [ x ]= C (a , b ) 4 x = factorial ( a ) /( factorial ( b ) * factorial (a - b ) ) 5 endfunction 6 disp ( p r o b a b i l i t y o f a d e f e c t i v e p a r t =2/20=0.1 ) 7 disp ( p r o b a b i l i t y o f a non d e f e c t i v e p a r t =0.9 ) 8 disp ( p r o b a b a i l i t y o f a t l e a s t t h r e e d e f e c t i v e s i n a
sample = )
9 1 -( C (20 ,0) *(0.9) ^20+ C (20 ,1) *(0.1) *(0.9) ^19+ C (20 ,2)
*(0.1) ^2*(0.9) ^18 ) 10 disp ( no . o f s a m p l e s h a v i n g t h r e e d e f e c t i v e p a r t s =1000 0.323= ) 11 1000*0.323
188
Chapter 35 Sampling and Inference
Scilab code Exa 35.1 calculating the SD of given sample

1 clc 2 disp ( s u p p o s e t h e c o i n 3 disp ( t h e n p r o b a b i l i t y 4 5 6 7 8 9 10 11 12 13 14
i s unbiased ) o f g e t t i n g t h e head i n a t o s s
=1/2 ) disp ( then , e x p e c t e d no . o f s u c c e s s e s =a =1/2 400 ) a =1/2*400 disp ( o b s e r v e d no . o f s u c c e s s e s =216 ) b =216 disp ( t h e e x c e s s o f o b s e r v e d v a l u e o v e r e x p e c t e d v a l u e= ) b-a disp ( S . D . o f s i m p l e s a m p l i n g = ( n p q ) 0.5= c ) c =(400*0.5*0.5) ^0.5 disp ( hence , z =(ba ) / c= ) (b - a ) / c disp ( a s z < 1 . 9 6 , t h e h y p o t h e s i s i s a c c e p t e d a t 5% level of s i g n i f i c a n c e )
Scilab code Exa 35.2 Calculating SD of sample 189
1 clc 2 disp ( s u p p o s e t h e d i e i s u n b i a s e d ) 3 disp ( t h e n p r o b a b i l i t y o f g e t t i n g 5 o r 6 w i t h one
d i e =1/3 )
4 disp ( then , e x p e c t e d no . o f s u c c e s s e s =a =1/3 9000 ) 5 a =1/3*9000 6 disp ( o b s e r v e d no . o f s u c c e s s e s =3240 ) 7 b =3240 8 disp ( t h e e x c e s s o f o b s e r v e d v a l u e o v e r e x p e c t e d
v a l u e= )
9 b-a 10 disp ( S . D . o f s i m p l e s a m p l i n g = ( n p q ) 0.5= c ) 11 c =(9000*(1/3) *(2/3) ) ^0.5 12 disp ( hence , z =(ba ) / c= ) 13 (b - a ) / c 14 disp ( a s z > 2 . 5 8 , t h e h y p o t h e s i s h a s t o be r e j e c t e d
at
1% l e v e l o f s i g n i f i c a n c e )
Scilab code Exa 35.3 Analysis of sample

1 2 3 4 5 6
clc p =206/840 disp ( q=1p ) q =1 - p n =840 disp ( s t a n d a r d e r r o r o f t h e p o p u l a t i o n o f f a m i l i e s h a v i n g a monthly i n c o m e o f r s . 250 o r l e s s =(p q / n ) 0.5= ) 7 ( p * q / n ) ^0.5 8 disp ( h e n c e t a k i n g 1 0 3 / 4 2 0 t o be t h e e s t i m a t e o f f a m i l i e s h a v i n g a monthly i n c o m e o f r s . 250 o r l e s s , t h e l i m i t s a r e 20% and 29% a p p r o x i m a t e l y )
190
Scilab code Exa 35.4 Analysis of sample

1 2 3 4 5 6 7 8 9 10 11 12 13 14
clear clc n1 =900 n2 =1600 p1 =20/100 p2 =18.5/100 disp ( p=(n1 p1+n2 p2 ) / ( n1+n2 ) ) p =( n1 * p1 + n2 * p2 ) /( n1 + n2 ) disp ( q=1p ) q =1 - p disp ( e =(p q ( 1 / n1+1/n2 ) ) 0 . 5 ) e =( p * q *((1/ n1 ) +(1/ n2 ) ) ) ^0.5 z =( p1 - p2 ) / e disp ( a s z < 1 , t h e d i f f e r e n c e b e t w e e n t h e p r o p o r t i o n s i s not s i g n i f i c a n t . )
Scilab code Exa 35.5 Checking whether real dierence will be hidden
1 2 3 4 5 6 7 8 9 10 11 12 13 14
clear clc p1 =0.3 p2 =0.25 disp ( q1=1 p1 ) q1 =1 - p1 disp ( q2=1 p2 ) q2 =1 - p2 n1 =1200 n2 =900 disp ( e =(( p1 q1 / n1 ) +(p2 q2 / n2 ) ) 0 . 5 ) e =(( p1 * q1 / n1 ) +( p2 * q2 / n2 ) ) ^0.5 z =( p1 - p2 ) / e disp ( hence , i t i s l i k e l y t h a t r e a l d i f f e r e n c e w i l l be h i d d e n . ) 191
Scilab code Exa 35.6 Checking whether given sample can be regarded as a random sample
1 clear 2 clc 3 disp ( m and n r e p r e s e n t s mean and number o f 4 5 6 7 8 9 10
objects
in sample r e s p e c t i v e l y ) m =3.4 n =900 M =3.25 d =1.61 disp ( z =(m M) / ( d / ( n 0 . 5 ) ) z =( m - M ) /( d /( n ^0.5) ) disp ( a s z > 1 . 9 6 , i t c a n n o t be r e g a r d e d a s a random sample )
Scilab code Exa 35.9 Checking whethet samples can be regarded as taken from the same population
1 clc 2 disp ( m1 and n1 r e p r e s e n t s mean and no . 3 4 5 6 7 8 9
of objects
in sample 1 ) disp ( m2 and n2 r e p r e s e n t s mean and no . o f o b j e c t s in sample 2 ) m1 =67.5 m2 =68 n1 =1000 n2 =2000 d =2.5 disp ( on t h e h y p o t h e s i s t h a t t h e s a m p l e s a r e drawn from t h e same p o p u l a t i o n o f d = 2 . 5 , we g e t ) 192
10 z =( m1 - m2 ) /( d *((1/ n1 ) +(1/ n2 ) ) ^0.5) 11 disp ( s i n c e | z | > 1 . 9 6 , t h u s s a m p l e s c a n n o t be
r e g a r d e d a s drawn from t h e same p o p u l a t i o n )
Scilab code Exa 35.10 calculating SE of dierence of mean hieghts

1 clc 2 disp ( m1 , d1 and n1 d e n o t e s mean , d e v i a t i o n and no .
of
o b j e c t s in f i r s t sample )
3 m1 =67.85 4 d1 =2.56 5 n1 =6400 6 disp ( m2 , d2 and n2 d e n o t e s mean , d e v i a t i o n and no . 7 8 9 10 11 12 13
of
o b j e c t s in second sample ) m2 =68.55 d2 =2.52 n2 =1600 disp ( S . E . o f t h e d i f f e r e n c e o f t h e mean h e i g h t s i s ) e =(( d1 ^2/ n1 ) +( d2 ^2/ n2 ) ) ^0.5 m1 - m2 disp ( | m1m2 | > 10 e , t h i s i s h i g h l y s i g n i f i c a n t . hence , t h e d a t a i n d i c a t e s t h a t t h e s a i l o r s a r e on t h e a v e r a g e t a l l e r than the s o l d i e r s . )
Scilab code Exa 35.12 Mean and standard deviation of a given sample
1 clear 2 clc 3 n =9 4 disp ( f i r s t
o f row d e n o t e s t h e d i f f e r e n t v a l u e s o f
sample ) 5 A (1 ,:) =[45 47 50 52 48 47 49 53 51]; 193
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
disp ( t h e s e c o n d row d e n o t e s t h e c o r r e s p o n d i n g deviation ) for i =1:9 A (2 , i ) = A (1 , i ) -48; end disp ( t h e t h i r d row d e n o t e s t h e c o r r e s p o n d i n g s q u a r e of deviation ) for i =1:9 A (3 , i ) = A (2 , i ) ^2; end disp ( t h e sum o f s e c o n d row e l e m e n t s = ) a =0; for i =1:9 a = a + A (2 , i ) ; end a disp ( t h e sum o f t h i r d row e l e m e n t s ) b =0; f o r i =1:9 b=b+A( 3 , i ) ; end b d i s p ( let m be the mean ) m=48+a / n d i s p ( let d be the standard deviation ) d =(( b /n ) (a / n ) 2 ) 0 . 5 t =(m 4 7 . 5 ) ( n 1) 0 . 5 / d
Scilab code Exa 35.13 Mean and standard deviation of a given sample
1 clc 2 disp ( d and n r e p r e s e n t s t h e d e v i a t i o n and no .
of
o b j e c t s in given sample ) 3 n =10 4 d =0.04 194
5 m =0.742 6 M =0.700 7 disp ( t a k i n g t h e h y p o t h e s i s t h a t t h e p r o d u c t
i n f e r i o r i . e . t h e r e i s no s i g n i f i c a n t b e t w e e n m and M ) 8 t =( m - M ) *( n -1) ^0.5/ d 9 disp ( d e g r e e s o f f r e e d o m= ) 10 f =n -1
i s not differene
Scilab code Exa 34.15 Standard deviation of a sample

1 clear 2 clc 3 n =11 4 disp ( t h e f i r s t row d e n o t e s t h e boy no . ) 5 A (1 ,:) =[1 2 3 4 5 6 7 8 9 10 11]; 6 disp ( t h e s e c o n d row d e n o t e s t h e marks i n t e s t 7 8 9 10 11 12 13 14 15 16 17 18 19 20
I ( x1
) ) A (2 ,:) =[23 20 19 21 18 20 18 17 23 16 19]; disp ( t h e t h i r d row d e n o t e s t h e marks i n t e s t I ( x2 ) ) A (3 ,:) =[24 19 22 18 20 22 20 20 23 20 17]; disp ( t h e f o u r t h row d e n o t e s t h e d i f f e r e n c e o f marks i n two t e s t s ( d ) ) for i =1:11 A (4 , i ) = A (3 , i ) -A (2 , i ) ; end disp ( t h e f i f t h row d e n o t e s t h e ( d 1) ) for i =1:11 A (5 , i ) = A (4 , i ) -1; end disp ( t h e s i x t h row d e n o t e s t h e s q u a r e o f e l e m e n t s o f f o u r t h row ) for i =1:11 A (6 , i ) = A (4 , i ) ^2; 195
21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
end A a =0; disp ( t h e sum o f e l e m e n t s o f f o u r t h row= ) for i =1:11 a = a + A (4 , i ) ; end a b =0; disp ( t h e sum o f e l e m e n t s o f s i x t h row= ) for i =1:11 b = b + A (6 , i ) ; end b disp ( s t a n d a r d d e v i a t i o n ) d =( b /( n -1) ) ^0.5 t =(1 -0) *( n ) ^0.5/2.24
196

Higher Engineering Mathematics - B. S. Grewal Companion Text

Uploaded by

Copyright:

Available Formats

Higher Engineering Mathematics - B. S. Grewal Companion Text

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Higher Engineering Mathematics - B. S. Grewal Companion Text

Uploaded by

Copyright:

Available Formats

Scilab Textbook Companion for Higher Engineering Mathematics by B. S.

161 171 189

List of Scilab Codes

82 82 83 85 85 86 86 87 87 88 88 89 89 90 90 91 92 94 94 94 95 95 95 96 96 96 97 97 97 98 98 98 99 99 99 100 100 100

Chapter 1 Solution of equation and curve tting

Scilab code Exa 1.1 nding the roots of quadratic equations

Scilab code Exa 1.11 forming an equation with known roots

Scilab code Exa 1.12 forming an equation under restricted conditions

14 x1 = -0.7784571 15 x2 =2.2891685 16 x3 =4.4892886 17 disp ( now , s i n c e we want e q u a t i o n whose sum o f

Figure 1.1: Finding the solution of equation by drawing graphs

disp ( from t h e graph , i t i s c l e a r t h a t t h e p o i n t o f i n t e r s e c t i o n i s n e a r l y x =1.43 )

Figure 1.2: Finding the solution of equation by drawing graphs

Figure 1.3: Finding the solution of equation by drawing graphs

Chapter 2 Determinants and Matrices

Scilab code Exa 2.1 Calculating Determinant

Scilab code Exa 2.2 Calculating Determinant

Scilab code Exa 2.3 Calculating Determinant

Scilab code Exa 2.4 Calculating Determinant

Scilab code Exa 5.8 Partial derivative of given function

Scilab code Exa 2.16 product of two matrices

Scilab code Exa 2.17 Product of two matrices

Scilab code Exa 2.18 Product and inverse of matrices

5 disp ( AB=C >B=i n v (A) C ) 6 B = inv ( A ) * C

Scilab code Exa 2.19 Solving equation of matrices

Scilab code Exa 2.20 Nth power of a given matrix

clc A =[11 -25;4 -9] n = input ( E n t e r t h e v a l u e o f n ) ; d i s p ( calculating A ^ n ) ; A n

Scilab code Exa 2.23 Inverse of matrix

Scilab code Exa 2.24.1 Rank of a matrix

Scilab code Exa 2.24.2 Rank of a matrix

Scilab code Exa 2.25 Inverse of matrix

Scilab code Exa 2.26 eigen values vectors rank of matrix

Scilab code Exa 2.28 Inverse of a matrix

Scilab code Exa 2.31 Solving equation using matrices

Scilab code Exa 2.32 Solving equation using matrices

disp ( e q u a t i o n s have i n f i n i t e no . o f s o l u t i o n s . ) end

Scilab code Exa 2.38 Inverse of a matrix

Scilab code Exa 2.39 Transpose and product of matrices

1 clear 2 clc 3 A =[5 4;1 2] 4 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f

2 clc 3 A =[3 1 4;0 2 6;0 0 5] 4 disp ( l e t R r e p r e s e n t s t h e m a t r i x o f

Scilab code Exa 2.45 eigen values and characteristic equation

Scilab code Exa 2.46 eigen values and characteristic equation

Scilab code Exa 2.47 eigen values and characteristic equation

and P r e p r e s e n t s a d i a g o n a l m a t r i x whose v a l u e s are the e i g e n v a l u e s o f A. then ) 7 [ R P ]= spec ( A ) 8 disp ( so , c a n o n i c a l form i s 2 x2+y2+ z 2 )

Scilab code Exa 2.52 Hermitian matrix

clc A =[2+ %i 3 -1+3* %i ; -5 %i 4 -2* %i ] disp ( A= ) A disp ( AA= ) A *A disp ( c l e a r l y ,AA i s h e r m i t i a n m a t r i x )

Scilab code Exa 2.53 tranpose and inverse of complex matrix

%i ) ] disp ( A= ) A disp ( AA= ) A *A disp ( AA= ) A * A disp ( i n v e r s e o f A i s ) inv ( A )

Scilab code Exa 2.54 Unitary matrix

Chapter 4 Dierentiation and Applications

Scilab code Exa 4.4.1 nding nth derivative

Scilab code Exa 4.5 nding nth derivative

Scilab code Exa 4.6 nding nth derivative

Scilab code Exa 4.7 nding nth derivative

Scilab code Exa 4.8 proving the given dierential equation

^2+ a ^2) * diff (y ,x , n ) ; disp ( n ) ; disp ( t h e e x p r e s s i o n f o r yn i s disp ( F ) ; disp ( Which i s e q u a l t o 0 ) ;

Scilab code Exa 4.9 proving the given dierential equation

n2+a 2 ) yn ) ; 12 // n=i n p u t ( E n t e r t h e o r d e r o f d i f f e r e n t i a t i o n ) ; 13 disp ( c a l c u l a t i n g yn f o r v a r i o u s v a l u e s o f n ) ; 14 for n =1:4

18 disp ( n ) ; 19 disp ( t h e e x p r e s s i o n f o r yn i s 20 disp ( F ) ; 21 disp ( Which i s e q u a l t o 0 ) ; 22 23 end 24 disp ( Hence p r o v e d ) ;

%i ) ] disp ( A= ) A disp ( AA= ) A A disp ( AA= ) A A disp ( i n v e r s e o f A i s ) inv ( A )

// q u e s 3 4 disp ( E q u a t i o n o f t a n g e n t ) ; syms x a t y xo = a ( cos ( t ) + t sin ( t ) ) ; yo = a ( sin ( t ) -t cos ( t ) ) ; s = diff ( xo , t ) / diff ( yo , t ) ; y = yo + s *( x - xo ) ; 50