Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

11 files changed, 3729 insertions, 0 deletions

diff --git a/Concise_Physics_by_H_Matyaka/1-Basics.ipynb b/Concise_Physics_by_H_Matyaka/1-Basics.ipynb
new file mode 100644
index 0000000..6c666dc
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/1-Basics.ipynb
@@ -0,0 +1,73 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 1: Basics"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.1: Neutral_temperature.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"x=(0:50:550)//temperature difference in x axis\n",
+"y=[0,0.43,0.79,1.10,1.36,1.54,1.69,1.77,1.80,1.78,1.70,1.54]//emf in y axis\n",
+"//calculation\n",
+"title('a graph of E vs teta')//setting title for graph\n",
+"xlabel('temperature difference teta')//setting x label\n",
+"ylabel('emf E')//setting y label\n",
+"plot(x,y)//plotting the graph\n",
+"printf('from the grapgh it can be determined that neutral temperature is 400deg C')\n",
+"x=(50:50:550)//temperature difference in x axis\n",
+"y=[8.6,7.9,7.3,6.8,6.2,5.6,5.1,4.5,4.0,3.4,2.8] //E/theta in y axis\n",
+"plot(x,y,'+-')//plotting the graph\n",
+"title('a graph of E/teta vs teta')//set title\n",
+"xlabel('temperature difference teta')//set x label\n",
+"ylabel(' E/teta')//set y label\n",
+"legend('E Vs Theta','E/theta Vs theta')\n",
+"b=-(4.5*10^-6)/400//gradient of graph is b \n",
+"a=4.5*10^-6-(b*400)//finding the intercept on y axis by substututing the points(400,4.5) in line equation\n",
+"printf('\n the value of b is %3.3e VdegC^-2',b)\n",
+"printf('\n the value of a is %3.3e VdegC^-1',a)\n",
+"\n",
+" "
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/10-Physical_Optics.ipynb b/Concise_Physics_by_H_Matyaka/10-Physical_Optics.ipynb
new file mode 100644
index 0000000..7557142
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/10-Physical_Optics.ipynb
@@ -0,0 +1,277 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 10: Physical Optics"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.1: plancks_theory.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"h=6.6*10^-34 //plancks constant\n",
+"c=3*10^8 //velocity of light\n",
+"e1=12.34//excited state\n",
+"e2=14.19//ground state\n",
+"//calculation\n",
+"l=(h*c)/((e2-e1)*1.6*10^-19)//conservation of energy and plancks theory\n",
+"//output\n",
+"printf('the wavelength is %3.3e m',l)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.2: wavelength_and_prism_angle.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"la=0.535*10^-6//wavelength\n",
+"nb=1.51//refractive index\n",
+"dmin=34 //minimum deviation\n",
+"//calculation\n",
+"l=la/nb//wavelength of light\n",
+"x=(nb-cosd(dmin/2))/sind(dmin/2)//refractive index of prism\n",
+"y=acotd(x)\n",
+"z=y*2\n",
+"//output\n",
+"printf('the wavelength of light is %3.3e m',l)\n",
+"printf('\nthe angle of prism is %3.0d deg',z)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.3: thin_film_interference.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"n=7//order of fringe\n",
+"l=0.63*10^-6 //wavelength\n",
+"x=24.8*10^-3 //seperation of bands\n",
+"d=1.5\n",
+"//calculation\n",
+"a=n*d*l/x//slit seperation\n",
+"//output\n",
+"printf('the slit seperation is %3.3e m',a)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.4: fringe_width_determination.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"n=6//order of fringe\n",
+"l=0.63*10^-6 //wavelength\n",
+"x=24.8*10^-3 //seperation of bands\n",
+"d=1.5\n",
+"a=2.7*10^-4\n",
+"//calculation\n",
+"x=d*(6+1/2)*l/a//distance between centre and sixth fringe\n",
+"w=l*1.6/a//fringe width\n",
+"//output\n",
+"printf('the distance between centre and sixth fringe is %3.3e m',x)\n",
+"printf('\nthe fringe width is %3.3e m',w)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.5: increasing_thickness_effect.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"a=4//widge dimension\n",
+"b=64//edge of tissue\n",
+"c=33//bright fringes\n",
+"l=0.53*10^-6 //wavelength\n",
+"//calculation\n",
+"m=b*c/a//number of bright fringes\n",
+"t=m*l/2//thickness\n",
+"//output\n",
+"printf('the thickness is %3.3e m and hence number of fringes also increases',t)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.6: wavelength_and_angular_displacement.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"n1=6//6th order image\n",
+"n2=5//5th order image\n",
+"n=3000//lines per cm\n",
+"//calculation\n",
+"l=n2*0.11*10^-6/(6-5)//applying dsinx=nl\n",
+"l1=l+(0.11*10^-6)//applying dsinx=nl\n",
+"d=1/(n*100)//applying dsinx=nl ,grating space calculation\n",
+"x=n1*l/d \n",
+"y=asind(x)\n",
+"//output\n",
+"printf('the wavenlength of first wave is %3.3e m',l)\n",
+"printf('\nthe wavenlength of second wave is %3.3e m',l1)\n",
+"printf('\n the angular displacement is %3.3f deg',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.7: wavelength_and_diffraction_angle.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"n2=1.36//refractive index\n",
+"N=5000*100 //number of lines per m\n",
+"t=23 //angle of diffraction\n",
+"//calculation\n",
+"l=sind(t)/(n2*N)//applying dsinx=nl,calculating wavelength \n",
+"x=N*l//angle of diffraction\n",
+"y=asind(x)\n",
+"//output\n",
+"printf('the wavelength of light in methanol is %3.3e m',l)\n",
+"printf('\n the angle of diffraction is %3.3f degrees',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.8: telescope_angular_magnification.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"fo=1.5//objective's focal length\n",
+"fc=0.04//eyepiece focal length\n",
+"//calculation\n",
+"m=fo/fc//angular magnification\n",
+"v=fc*(fc+fo)/fo//distance of eye ring from eyepiece\n",
+"//output\n",
+"printf('the angular magnification is %3.2f',m)\n",
+"printf('\n the distance of eye ring from eyepiece is %3.3f m',v)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/11-Semiconductors.ipynb b/Concise_Physics_by_H_Matyaka/11-Semiconductors.ipynb
new file mode 100644
index 0000000..bafe760
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/11-Semiconductors.ipynb
@@ -0,0 +1,151 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 11: Semiconductors"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.1: rms_current_and_peak_pd.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"vp=50//ac source supply\n",
+"r1=35\n",
+"r2=1450 //resistors\n",
+"//calculation\n",
+"vs=4*vp//transformer equation\n",
+"i=100/(r1+r2)//peak current\n",
+"irms=i/sqrt(2)//rms current\n",
+"v0=100*r1/(r1+r2)\n",
+"pp=100-v0//peak pd\n",
+"//output\n",
+"printf('the rms value of current is %3.3f A',irms)\n",
+"printf('\n the peak pd is %3.3f V',pp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.2: common_emittor_transistor.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"vbe=1.2//pd across emitter\n",
+"ib=120*10^-6//base current\n",
+"v1=1.5//final voltafe\n",
+"i2=175*10^-6//increased current\n",
+"//calculation\n",
+"r=vbe/ib//static input resistence\n",
+"h=(v1-vbe)/(i2-ib)//input hybrid parameter\n",
+"//output\n",
+"printf('the static input resistence is %3.0e ohm',r)\n",
+"printf('\nthe input hybrid parameter is %3.3e ohm',h)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.3: common_base_transistor.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v1=7.5//initial voltag\n",
+"v2=11.5//final voltage\n",
+"ic=18*10^-6//collector current\n",
+"//calculation\n",
+"r=(v2-v1)/ic//output resistance\n",
+"//output\n",
+"printf('the output resistance is %2.2e ohm ',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.4: common_emittor_amplifier.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"vbe=2.5//voltage across base-emitter\n",
+"hfe=75//current gain\n",
+"rb=75*10^3 //base current\n",
+"//calculation\n",
+"rc=5*rb/(vbe*hfe)//collector load resistance\n",
+"//output\n",
+"printf('the collector load resistance is %2.2e ohm',rc)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/2-Mechanics.ipynb b/Concise_Physics_by_H_Matyaka/2-Mechanics.ipynb
new file mode 100644
index 0000000..28799cc
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/2-Mechanics.ipynb
@@ -0,0 +1,886 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 2: Mechanics"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.10: resultant_force.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc \n",
+"clear\n",
+"//input magnitude of forces\n",
+"f1=40\n",
+"f2=50\n",
+"//calculation\n",
+"d=50^2+40^2+2*50*40*cosd(50)//finding the diagonal\n",
+"r=50^2+40^2+2*50*(40)*cosd(130)//reversing the side and finding diagonlprintf('the resultant is %3.3f',d1)\n",
+"r1=sqrt(r)//resultant sum\n",
+"d1=sqrt(d)// resultant when smaller force is subtracted from larger\n",
+"//output\n",
+"printf('1. the resultant sum is %3.3f N',d1)\n",
+"printf('\n 2. the resultant when smaller force is subtracted from larger is %3.3f N',r1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11: components_of_velocity.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=380//velocity\n",
+"//calculation\n",
+"vh=v*cosd(60)//horizontal component\n",
+"vv=v*sind(60)//vertical component\n",
+"//output\n",
+"printf('the horizontal component is %3.3f ms^-1',vh)\n",
+"printf('\nthe vertical component is %3.3f ms^-1',vv)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.12: components_of_force.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"fc=50//force applied by magnet\n",
+"x=90-20 //angle of force\n",
+"//calculation\n",
+"fb=fc*sind(70)//force due to b\n",
+"fa=fc*cosd(70)//force due to a\n",
+"//output\n",
+"printf('the force due to b is %3.3f N',fb)\n",
+"printf('\nthe force due to b is %3.3f N',fa)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.13: inelastic_collission.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m1=1\n",
+"v1=25\n",
+"m2=2\n",
+"v2=0\n",
+"//calculation\n",
+"v=(m1*v1)+(m2*v2)//applying princilpe of conservation of linear momentum\n",
+"v4=v/(m1+m2)\n",
+"//output\n",
+"printf('the velocity with which both will move is %3.3f ms^-1',v4)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.14: Inelastic_collission.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m1=1//mass of object 1\n",
+"v1=25//velocity of object 1\n",
+"m2=2//mass of object 2\n",
+"v2=0//body at rest,velocity =0\n",
+"v3=10\n",
+"//caclulation\n",
+"u=((m1*v1)+(m2*v2)-(m2*v3))/2//applying princilpe of conservation of linear momentum\n",
+"//output\n",
+"printf('\n the value of u is %3.3f ms^-1',-u)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.15: angularvelocity_and_centripetal_force.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=2//mass\n",
+"r=4//radius\n",
+"v=6//uniform speed\n",
+"//calculation\n",
+"w=v/r//angular velocity\n",
+"f=m*r*w*w//centripetal force\n",
+"//output\n",
+"printf('the angular velocity is %3.3f rads^-1',w)\n",
+"printf('\n the centripetal force is %3.3f N',f)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.16: tension_in_arm.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=140//mass\n",
+"v=8//speed\n",
+"r=5//radius\n",
+"g=9.8//acceleration due to gravity\n",
+"//calculation\n",
+"t=((m*v^2/5)^2)+(140*9.8)^2 //applying parallelogram of vectors\n",
+"t1=sqrt(t)\n",
+"//output\n",
+"printf('the tension in arm is %3.3f N',t1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.17: inclination_and_reaction.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=15//velocity\n",
+"m=70//mass\n",
+"r=50//radius\n",
+"//calculation\n",
+"x=v*v/(r*10)//applying parallelogram of vectors,then for equilibrium\n",
+"y=atand(x)\n",
+"r1=(m*10)/cosd(y)\n",
+"//output\n",
+"printf('the inclination is %3.3f deg',y)\n",
+"printf('\n the reaction is %3.3f N',r1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.18: planet_mean_density.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"r=5500//radius\n",
+"g1=6.7*10^-11\n",
+"g=7//acceleration due to gravity\n",
+"//calculation of mean density\n",
+"p=3*g/(4*%pi*r*10^3*g1)//mean density\n",
+"//output\n",
+"printf('the mean density of planet is %3.3f kgm^-3',p)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.19: orbit_radius_and_linearvelocity.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=5*10^24//mass of earth\n",
+"g1=6.7*10^-11\n",
+"//calculation\n",
+"r=((6.7*10^-11*5*10^24*(3600*24)^2)/(4*%pi^2))^(1/3)//orbit radius\n",
+"v=(2*%pi*r)/(3600*24)//linear velocity\n",
+"//output\n",
+"printf('the orbit radius is %3.3f',r)\n",
+"printf('\n the linear velocity is %3.3f ms^-1',v)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.1: acceleration_and_distance.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input from given graph\n",
+"//calculation of initial accleration\n",
+"ia=18/4\n",
+"// calculation of final accleration\n",
+"fa=-18/10\n",
+"decel=-(fa)//calculation of deceleration\n",
+"//calculation of total distance covered\n",
+"d=0.5*(4*18)+(8*18)+0.5*(10*18)//area under velocity time graph\n",
+"//output\n",
+"printf('\n the initial acceleration is %3.3f m/s^2',ia)\n",
+"printf('\n the final acceleration is %3.3f m/s^2',decel)\n",
+"printf('\n the distance covered is is %3.3f m',d)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.20: mass_of_galaxy.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=3*10^5//orbit speed\n",
+"r=4.6*10^20//distance\n",
+"g1=6.7*10^-11\n",
+"//calculation of mass\n",
+"m=v*v*r/g1 //Newtons law\n",
+"//output\n",
+"printf('the mass is %2.3e kg',m)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.21: total_kinetic_energy.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=0.6//speed\n",
+"m=0.3//mass\n",
+"//calculation\n",
+"e=0.75*m*v*v//total kinetic energy is kinetic energy+moment of inertia\n",
+"//output\n",
+"printf('the total kinetic energy is %3.3f J',e)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.22: time_taken_to_move.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"t1=34\n",
+"u=0//starts from rest\n",
+"x=3//distance to move\n",
+"//calculation\n",
+"t=(3*3/(10*sind(t1)))^0.5//from law of conservation of energy\n",
+"//output\n",
+"printf('the time taken is %3.3f s',t)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.23: angular_velocity_ratio.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"i1=53//inertia when it spins with panels carrying solar cells\n",
+"i2=37//inertia about axis of rotation\n",
+"//calculation\n",
+"r=i1/i2//law of conservation of angular momentum\n",
+"//output\n",
+"printf('the ratio of angular velocities is %3.3f',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.25: attributes_of_shm.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f=9//frequency\n",
+"x=0//at midpoint of stroke x=0\n",
+"//calculation\n",
+"t=1/f\n",
+"a=-4*%pi^2*f^2*x//acceleration for shm\n",
+"v=2*%pi*f*0.05//velocity for shm\n",
+"a1=-4*%pi^2*9^2*0.05//acceleration at amplitude\n",
+"v1=0//velocity at amplitude is 0\n",
+"//output\n",
+"printf('the period of oscillation is %3.3f s^-1',t)\n",
+"printf('\n the velocity  at midpoint of stroke is %3.3f ms^-1',v)\n",
+"printf('\n the acceleration at midpoint of stroke is %3.3f ms^-2',a)\n",
+"\n",
+"printf('\n the velocity at amplitude is %3.3f  ms^-1',v1)\n",
+"printf('\n the acceleration at amplitude is %3.3f ms^-2',a1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.26: simple_harmonic_motion.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"g=10\n",
+"t=0.3//period of shm\n",
+"//calculation\n",
+"x=g*t^2/(4*%pi^2)//for shm maximum amplitude\n",
+"//output\n",
+"printf('the maximum amplitude for bead to be in contact is %3.3f m',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.27: attrbutes_simple_pendulum.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p1=2.3//period of pendulum\n",
+"p2=3.1//period when pendulum is lengthened\n",
+"//calculation\n",
+"g=4*%pi^2/(p2^2-p1^2)//acceleration of free fall\n",
+"l=p1^2*g/(4*%pi^2)//length of pendulum\n",
+"//output\n",
+"printf('the acceleration of free fall is %3.3f m/s^2',g)\n",
+"printf('\n the length of pendulum is %3.3f m',l)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.28: maximum_displacement_shm.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"f=55 //frequency\n",
+"a=7*10^-3 //amplitude\n",
+"\n",
+"\n",
+"//calculation\n",
+"a=(-2*%pi*f)^2*a\n",
+"\n",
+"//output\n",
+"printf('the acceleration of the body when it is at its maximum displacement from its zero position is -%3.1f ms^-2',a)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.29: maximum_potential_energy_shm.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f=55//frequency\n",
+"amp=7*10^-3//amplitude\n",
+"m=1.2//mass\n",
+"//calculation\n",
+"e=0.5*m*4*%pi^2*f^2*amp^2//maximum pe occurs at zero position\n",
+"//output\n",
+"printf('the maximum pe is %3.3f J',e)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.2: acceleration_and_distance.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=0 //car stops => final velocity=0\n",
+"u=29 //initial velocity\n",
+"t=11 //time \n",
+"//calculation of acceleration\n",
+"a=(v-u)/t//eqn of uniformly accelerated body\n",
+"//calculating distance travelled during this period\n",
+"d=(v+u)*t*0.5//eqn of uniformly accelerated body\n",
+"//output\n",
+"printf('the accleration is %3.3f ms^-2 ',a)\n",
+"printf('\nthe distance travelled is %3.3f m',d)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.30: extension_of_steel_wire.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"l=6.5//length\n",
+"m=0.06//mass of wire\n",
+"m1=10//mass attached\n",
+"g=9.8//acceleration due to gravity\n",
+"e=2.1*10^11//youngs modulus\n",
+"ro=8*10^3//density of steel\n",
+"//calculation\n",
+"e1=m1*g*ro*l*l/(e*m)//extension caused \n",
+"pe=0.5*g*m1*e1//potential energy \n",
+"//output\n",
+"printf('the extension caused is %3.3e m',e1)\n",
+"printf('\n the potential energy is %3.3f J',pe)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.31: Youngs_modulus.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"w=250*10^3\n",
+"s=0.00003//strain\n",
+"a=0.04//area\n",
+"w1=320*10^3\n",
+"//calculation\n",
+"e=w/(a*s)//youngs module\n",
+"st=w1/a//stress\n",
+"//output\n",
+"printf('the youngs modulus is %3.3e N/m^2',e)\n",
+"printf('\n the stress is %3.0e N/m^2',st)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.32: wire_length_change.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=40//mass\n",
+"g=9.8//acceleration due to gravity\n",
+"E=2*10^11//youngs modulus\n",
+"//calculation\n",
+"t1=m*g/5//principle of momentum\n",
+"t2=4*m*g/5 //principle of momentum\n",
+"d=4*(t2-t1)/(4*%pi*10^-6*E)//difference in length\n",
+"//output\n",
+"printf('the difference is %3.0e m',d)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.3: time_to_reach_aircraft.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=120 //velocity\n",
+"a=75 //accleration\n",
+"//calculation of time\n",
+"t=2*v/(a*cosd(45))//eqn of uniformly accelerated body\n",
+"//output\n",
+"printf('the time taken is %3.3f s',t)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.4: resultant_force.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f1=50\n",
+"f2=50\n",
+"//calculation of net force\n",
+"f=f1+f2 // the two forces act in same direction\n",
+"//output\n",
+"printf('the resultant force is %3.3f N',f)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.5: car_and_wind.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"vc=25 //velocity of car\n",
+"va=10 //velocity of wind\n",
+"va1=15 //velocity of wind westward\n",
+"//calculation\n",
+"v1=vc+va//resultant velocity for a tail of wind\n",
+"v2=vc-va//when wind blows westward at 10 m/s^resultant velocity \n",
+"v3=vc-va1//resultant velocity when wind blows westward at 15m/s^2\n",
+"//output\n",
+"printf('1. the resultant velocity of wind is %3.3f ms^-1 eastwards ',v1)\n",
+"printf('\n2. the resultant velocity of wind is %3.3f ms^-1 westwards ',v2)\n",
+"printf('\n3. the resultant velocity of wind is %3.3f  ms^-1westwards ',v3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.7: velocity_of_speedboat.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc \n",
+"clear\n",
+"//input\n",
+"v=30 //velocity of speedboat\n",
+"vw=40 //velocity of wind\n",
+"//calculation\n",
+"x=(30/40)//angle between original velocity of boat and resultant velocity\n",
+"y=atand(x)//applying trigonometry\n",
+"b=90+y//bearing of boat\n",
+"//output\n",
+"printf('the bearing of speedboat is %3.3f deg',b)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.8: tension_on_string.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc \n",
+"clear\n",
+"//input\n",
+"f1=6 //tension on string AB\n",
+"f2=6 //tension on string BC\n",
+"//calculation of tension\n",
+"t=2*f1*sind(55)// the resultant tension is the diagonal of rhombus formed\n",
+"//output\n",
+"printf('/n the resultant tension is %3.3f N',t)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/3-Waves.ipynb b/Concise_Physics_by_H_Matyaka/3-Waves.ipynb
new file mode 100644
index 0000000..d4a2733
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/3-Waves.ipynb
@@ -0,0 +1,211 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 3: Waves"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.1: refraction_and_incidence_angle.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//calculation of angle of refraction\n",
+"rj=(sind(6)/0.76)//from snells law\n",
+"x=asind(rj)\n",
+"printf('the refractive index of jelly is %3.3f deg',x)\n",
+"// calculating angle of incidence\n",
+"printf('\nsince angle of refraction and angle of incidence are alternate angles , angle of incidence is %3.3f deg',x)\n",
+"//calculating angle of refraction\n",
+"np=0.59/0.76 // according to relationship of media\n",
+"jnp=sind(7.9)/0.78\n",
+"rp=asind(jnp)\n",
+"printf('\nthe angle of refraction  is %3.3f deg',rp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.2: critical_angle.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input data\n",
+"a=1.28 //refractive index of X\n",
+"b=1.41 //refractive index of Y\n",
+"//calculation of condition for total internal reflection\n",
+"x=(a/b) \n",
+"c=asind(x) // calculating critical angle\n",
+"//output\n",
+"printf('light incident with an angle greater than %3.3f degrees will be totally internally reflected',c)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.3: wavespeed_in_medium.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input data\n",
+"nb=0.67 //refractive index\n",
+"va=3.45*10^3\n",
+"//calculation\n",
+"vb=va/nb //snells law\n",
+"//output\n",
+"printf('the speed of the wave in medium b is %3.3f m/s',vb)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.4: frequency_for_antinode.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input data\n",
+"f=120 //lowest frequency \n",
+"//calculation\n",
+"x=3*f // the next higher frequency is thrice the lowest frequency\n",
+"//output\n",
+"printf('the next higher frequency where the antinode is formed is at %3.3f Hz',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.5: wave_frequency_speed.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input data\n",
+"amp=3.4*10^-5 //amplitude of the wave\n",
+"af=5.7*10^2 //angular frequency\n",
+"k=20 //wavenumber\n",
+"//calculation\n",
+"//wave frequency\n",
+"f=af/(2*%pi)\n",
+"l=(2*%pi)/k\n",
+"v=f*l\n",
+"printf('the wave frequency is %3.3f and the speed is %3.3f  m/s',f,v)\n",
+"//calculating greatest speed for the wave to pass through\n",
+"vmax=af*amp //greatest speed\n",
+"//output\n",
+"printf('\nthe greatest value of speed for the wave to pass through is %3.3f m/s',vmax)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.6: wave_attributes.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"k=16\n",
+"w=23\n",
+"//calculation\n",
+"//1.wavelength\n",
+"l=2*%pi/k\n",
+"//output\n",
+"printf('the wavelength is %3.3f m',l)\n",
+"//2.wavespeed\n",
+"v=(l*w)/(2*%pi)\n",
+"printf('\nthe wavespeed is %3.3f m/s',v)\n",
+"//3.pase difference\n",
+"pha=(0.5*2*%pi)/0.39 // phase difference of molecules 0.5m apart \n",
+"printf('\n the phase difference is %3.3f radians',pha)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/4-Waves.ipynb b/Concise_Physics_by_H_Matyaka/4-Waves.ipynb
new file mode 100644
index 0000000..07bb5ed
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/4-Waves.ipynb
@@ -0,0 +1,212 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 4: Waves"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.1: amplitude_and_pressure_change.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"w=1.8 //wavelength\n",
+"//calulation\n",
+"y=sind(15*360/180)//displacement at 15cm from reflector\n",
+"//output\n",
+"printf('1. at 45cm, antinode occurs and hence pressure is minimum')\n",
+"printf('\n 2. at 90cm node arises and hence pressure is maximum')\n",
+"printf('\n 3. at 15cm frm reflector the displacement is %3.3f',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.3: length_of_tube.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"f=520 //frequency\n",
+"t2=293 //air temperature to produce fundamental +273\n",
+"t1=273// 0deg C\n",
+"v1=330//speed of sound waves\n",
+"//calculation\n",
+"v2=330*(293/273)^0.5 //speed at 20 deg C\n",
+"l=v2/f//wavelength\n",
+"len=l/4 - 0.01 //length\n",
+"//output\n",
+"printf('the length of tube is %3.3f m',len)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.4: frequency_of_beats.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"v1=330 //speed of sound\n",
+"t3=303 //fundamental temperature for the air\n",
+"t1=273// 0deg C\n",
+"//calculation\n",
+"v3=v1*(t3/t1)^0.5 //new speed of sound\n",
+"f=v3/0.66 //frequency\n",
+"fb=f-520 //frequency of beats\n",
+"//output\n",
+"printf('the frequency of beats is %3.3f Hz',fb)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.5: fundamental_frequency.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"T=100 //tension\n",
+"l=1.5 //length\n",
+"m=0.3*10^-6 //mass\n",
+"//calculation\n",
+"f=(T/(m/l))^0.5/(2*l)//fundamental frequency produced \n",
+"//output\n",
+"printf('the fundamental frequency produced is %3.3f Hz',f)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.6: doppler_effect.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"f=150 //frequency \n",
+"v=320 //speed of sound\n",
+"ul=11 //speed with which listener approaches\n",
+"us=7 //speed of source\n",
+"//calculation\n",
+"fa=f*v/(v-us)//doppler effect\n",
+"fa1=(v+ul)*f/(v)//doppler effect\n",
+"fa2=(v+ul)*f/(v-us)//doppler effect\n",
+"//output\n",
+"printf('frequency when source moves at 7ms^-1 %3.3f Hz',fa)\n",
+"\n",
+"printf('\n frequency when listener moves at 11ms^-1 %3.3f Hz',fa1)\n",
+"printf('\n frequency when source moves at 7ms^-1 and listener at 11ms^-1 %3.3f Hz',fa2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.7: apparent_frequency_change.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"us=264 //speed of jet fighter\n",
+"x=71.7\n",
+"v=340 //velocity of sound\n",
+"f=1*10^3 //frequency\n",
+"//calculation\n",
+"usd=us*cosd(x)//horizontal component of velocity\n",
+"fr= (v*f)/(v-usd) -((v*f)/(v+usd))//frequency range ,doppler effect\n",
+"//output\n",
+"printf('the frequency range is %3.3f Hz ',fr)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/5-Light.ipynb b/Concise_Physics_by_H_Matyaka/5-Light.ipynb
new file mode 100644
index 0000000..d8aaa79
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/5-Light.ipynb
@@ -0,0 +1,217 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 5: Light"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.1: minimum_deviation.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"np=1.39 //refractive index of prism\n",
+"nl=1.29 //refractive index of liquid\n",
+"a=62 //refracting angle of prism\n",
+"//calculation\n",
+"x=np*sind(62/2)/nl//snells law\n",
+"y=asind(x)\n",
+"d=(y*2)-a//minimum deviation\n",
+"//output\n",
+"printf('the minimum deviation is %3.3f degree',d)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.2: incidence_and_prism_angle.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"np=1.39 //refractive index in air\n",
+"a=62 //refracting angle of prism\n",
+"//calculation\n",
+"x=1/np\n",
+"c=asind(x)//critical angle\n",
+"r=a-c\n",
+"i= np* sind(r)//snells law\n",
+"i1=asind(i)\n",
+"A=2*c//greatest prism angle allowing refraction\n",
+"//output\n",
+"printf('angle of incidence producing maximum deviation is %3.3f deg',r)\n",
+"printf(' \n greatest prism angle allowing refraction is %3.3f deg',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.3: position_and_nature_of_image.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f=0.15 //focal length\n",
+"u=0.2 //distance of object\n",
+"//calculation\n",
+"x=(1/-f)-(1/u)//lens formula\n",
+"y=1/x\n",
+"m=y/u//linear magnification\n",
+"//output\n",
+"printf('the position of image is %3.3f m',y)\n",
+"printf('\n linear magnification is %3.3f hence image is diminished',m)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.4: position_of_image.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f1=0.25 //focal length of diverging lens\n",
+"f2=0.2 //focal length of converging lens\n",
+"//calculation\n",
+"x=(1/-f1)+(1/f2)//lens formula\n",
+"y=1/x\n",
+"a=(1/y)-(1/0.15)//lens formula\n",
+"b=1/a\n",
+"//output\n",
+"printf('the position of image is %3.3f m hence the image is virtual',b)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.5: position_and_nature_of_image.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"f=0.5 //focal length \n",
+"u=0.8 //distance of object \n",
+"f1=0.2 //focal length of converging lens\n",
+"d=1 //distance behind the first lens\n",
+"//calculation\n",
+"x=(1/f)-(1/u)//lens formula\n",
+"y=1/x\n",
+"u1=-(y-d)//second lens \n",
+"a=1/f1 +(1/-u1)//lens formula\n",
+"b=1/a\n",
+"//output\n",
+"printf('the image lies %3.3f m behind second lens',b)\n",
+"printf('\n the image is %3.3f  m behind first lens',b+d)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.6: lens_values.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"F=5 //power of lenses\n",
+"f1=0.45 //focal length\n",
+"//calculation\n",
+"x=F-(1/f1)//lens formula\n",
+"f2=1/x\n",
+"//output\n",
+"printf('the focal length is %3.3f m',f2)\n",
+"printf('\n the power is %3.3f dioptre',x)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/6-Heat.ipynb b/Concise_Physics_by_H_Matyaka/6-Heat.ipynb
new file mode 100644
index 0000000..fd85183
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/6-Heat.ipynb
@@ -0,0 +1,446 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 6: Heat"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.10: gas_external_work.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=3*10^-2 //mass of water\n",
+"r1=1*10^3//density of water\n",
+"r2=0.5//density of steam\n",
+"p=1.01*10^5//atmospheric pressure\n",
+"//calculation\n",
+"v1=m/r1//volume of water\n",
+"v2=m/r2//volume of gas\n",
+"w=(v2-v1)*p//external work done by gas\n",
+"//output\n",
+"printf('the work done is %3.0f J',w)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.12: platinum_resistance_theromoeter.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"r100=6.9//resistence of steam\n",
+"r0=5.8 //resistece of ice\n",
+"t=550 //temperature\n",
+"//calculation\n",
+"r=(t*(r100-r0))/100 +5.8//platinum resistance thermometre\n",
+"//output\n",
+"printf('the resistence is %3.3f ohm',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.14: length_at_temperature.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"l=11.7//length of thermometer at steam\n",
+"l0=3.4*10^-2//length of thermometer at ice\n",
+"//calculation\n",
+"x=0.034+0.034*(0.244*10^-3*45^2)//length of temperature on standard scale\n",
+"//output\n",
+"printf('thread length is %3.3f m',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.15: heat_transfer_rate.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"a=5 //area\n",
+"k=0.07 //thermal conductivity\n",
+"dt=21 //temperature difference\n",
+"x= 4*10^-3 //thickness of wood\n",
+"//calculation\n",
+"y=-(k*a*dt/x)//steady state equation\n",
+"//output\n",
+"printf('the rate of transfer is %3.3f Js^-1',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.16: temperature_gradient.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"d=3*10^-3//thickness of sheet\n",
+"l=12*10^-3//seperated distance\n",
+"//calculation\n",
+"x=1/40//law of conservation of energy\n",
+"y=x*d/l//from x\n",
+"//output\n",
+"printf('the ratio of temperature gradient in rubber to polystyrene is %3.3f0',x)\n",
+"printf('\nthe ratio of temperature difference across rubber and polystyrene is %3.3e',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.1: heat_given_out.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"m=0.5 //mass\n",
+"c=460 //specific heat capacity of iron\n",
+"t1=70//initial temperature\n",
+"t2=10//final temperature\n",
+"//calculation\n",
+"q=m*c*(t1-t2)//heat required \n",
+"//output\n",
+"printf('the heat required is %3.0f J',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.2: potential_difference_heater.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"T=100 //rise in temperature\n",
+"i=2.7 //current\n",
+"t=950 //time taken\n",
+"mc=0.15//mass of calorimeter \n",
+"cy=3*10^3//specific heat capacity of y\n",
+"cc=2*10^3//specific heat capacity of calorimeter\n",
+"my=160*10^-3//mass of liquid\n",
+"//calculation\n",
+"v=((my*cy)+(mc*cc))*T/(i*t)//law of conservation of heat\n",
+"//output\n",
+"printf('the potential difference is %3.0f V',v)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.3: heat_loss_and_specific_heat.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"iw=4.5 //current\n",
+"vw=5.2 //pd of water\n",
+"mw=6*10^-2 //flow of water\n",
+"cw=4.18*10^3 //heat capacity of water\n",
+"ix=5.5//current of x\n",
+"iv=7.7//pd of x\n",
+"im=18*10^-2//flow of x\n",
+"//calculation\n",
+"x=(iw*vw)-((mw*cw*5)/60)//rate of heat loss\n",
+"cx=(6*4180)/18 +1263//specific heat capacity of x\n",
+"//output\n",
+"printf('the rate of heat loss is %3.3f W',x)\n",
+"printf('\n the specific heat of x is %3.3e Jkg^-1K^-1',cx)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.4: Boyles_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+" clc\n",
+"clear\n",
+"//input\n",
+"v1=0.52 //volume of ideal gas\n",
+"p1=2.3*10^5 //pressure of ideal gas\n",
+"p2=6.7*10^5 //pressure changed\n",
+"//calculation\n",
+"v2=p1*v1/p2//boyle's law\n",
+"//output\n",
+"printf('the volume is %3.3f m^3',v2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.5: Charles_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v2=11.3 //final volume\n",
+"v1=7.8//initial volume\n",
+"t1=67+273 //initial temperature\n",
+"//calculation\n",
+"t2=v2*t1/v1//charles law\n",
+"//output\n",
+"printf('the final temperature is %3.0d K',t2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.6: pressure_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p1=1.01*10^5//initial pressure\n",
+"t2=135+273//final temperature\n",
+"t1=273//initial temperature\n",
+"d=2.8 //density\n",
+"//calculation\n",
+"p2=p1*t2/t1//pressure law\n",
+"p=(3*p2/2.8)^0.5//kinetic theory\n",
+"//output\n",
+"printf('rms speed of gas molecule is %3.0f m/s',p)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.7: KE_and_rms_velocity.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"t1=273//initial tenperature\n",
+"t2=408//final temperature\n",
+"//calculation\n",
+"e=t1/t2//ratio of mean molecuar KE\n",
+"c1=402*sqrt(0.67)//rms speed\n",
+"//output\n",
+"printf('the ratio of kinetic energy is %3.3f',e)\n",
+"printf('\n  the rms speed of gas molecule is %3.0f ms^-1',c1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.8: ideal_gas_equatio.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p=1.01*10^7 //pressure of gas\n",
+"v=0.1 //volume of gas\n",
+"R=8.3\n",
+"T=280//temperature\n",
+"g=0.017//mass of 1 mole\n",
+"d=1100//density\n",
+"//calculation\n",
+"n=p*v/(R*T)//ideal gas equation\n",
+"m=n*g//mass of gas\n",
+"v=m/d//volume occupied\n",
+"//output\n",
+"printf('the volume is %3.3e m^3',v)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.9: Boyles_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p1=9*10^4//total pressure\n",
+"x=1*10^4//water pressure\n",
+"//calculation\n",
+"p2=(p1-x)/2//boyles law\n",
+"p=p2+x//adding vapour pressure\n",
+"//output\n",
+"printf('the final pressure is %3.0e Pa',p)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/7-Electricity.ipynb b/Concise_Physics_by_H_Matyaka/7-Electricity.ipynb
new file mode 100644
index 0000000..7542c66
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/7-Electricity.ipynb
@@ -0,0 +1,272 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: Electricity"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.1: Electric_potential_strength.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"e=1.6*10^-19 //charge of electron\n",
+"r=0.075*10^-3 // radius of electron\n",
+"ep=8.85*10^-12 //permittivity of free space\n",
+"//calculation\n",
+"v=-e/(4*%pi*ep*r)//electric potential\n",
+"e=-e/(4*%pi*ep*r*r)//electric field strength\n",
+"//output\n",
+"printf('resultant potential is %3.3e V',v)\n",
+"printf('\n resultant electric field strength %3.3f V/m',e)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.2: ratio_of_force.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"q=2.4*10^-19 //charge1\n",
+"Q=3.8*10^-19//charge2\n",
+"ep=8.85*10^-12//permittivity of free space\n",
+"G=6.7*10^-11\n",
+"m=8.9*10^-31//mass 1\n",
+"M=1.5*10^-30//mass 2\n",
+"//calculation\n",
+"x=q*Q/(4*%pi*ep*m*M*G)//coulumbs law\n",
+"//output\n",
+"printf('the ratio of electrostatic force between charges %3.3e',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.3: emf_and_internal_resistance.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"i=0.5 //current in circuit\n",
+"R=6 //resistance of circuit\n",
+"i1=0.3//dropped current\n",
+"//calculation\n",
+"r=1.2/0.2\n",
+"e=i*(r+R)//ohms law\n",
+"//output\n",
+"printf('the battery emf is %3.3f V',e)\n",
+"printf('\n the internal resistence is %3.3f ohm',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.4: power_output.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"d=8.2*10^-7 //resistivity of coil\n",
+"l=15 //length of wire\n",
+"r=0.3*10^-3 //radius of wires\n",
+"v=160 //power output\n",
+"//calculations\n",
+"R=d*l/(%pi*r*r)\n",
+"p=v*v/R //for one coil\n",
+"p1=v*v/(R+R) //for two coils in series\n",
+"rp=(R*R)/(R+R)//total resistence \n",
+"pp=(v*v)/rp//total power\n",
+"//output\n",
+"printf('the power when one coil is %3.3f W',p)\n",
+"printf('\nthe power when two coils in series is %3.3f W',p1)\n",
+"printf('\n the power when coils in parallel is %3.3f W',pp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: percent_of_pd.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"r1=40//resistance 1\n",
+"r2=20//resistamce 2\n",
+"r3=10//resistance 3\n",
+"v=1.6//voltage\n",
+"//calculation\n",
+"R=r1+r2+r3//total resistance in series\n",
+"x=((v*r1)*70)/((2*50)*(1.6*40))//fraction of pd\n",
+"x=x*100//percentage pd\n",
+"//output\n",
+"printf('the percentage of pd is %3.0f percent',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: final_resistance_calculation.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"a=4.3*10^-3//temperature co-efficient of resistance\n",
+"//calculation\n",
+"r2=((60*a+1)/(20*a+1))*10//resistance \n",
+"//output\n",
+"printf('the final resistence is %3.3f ohm',r2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.7: internal_resistance_calculation.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"l1=82.3//balance length  with switch open\n",
+"l2=75.8//balance length with switch closed\n",
+"R=9//resistance\n",
+"//calculation\n",
+"r=(R*l1/l2)-R//internal resistance\n",
+"//output\n",
+"printf('the internal resistence is %3.3f ohm',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.8: calculation_of_resistance.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p=2*10^-6//pd across wire\n",
+"v=1.5//voltage \n",
+"l=1.5*10^3//length of potentiometer\n",
+"R=7//resistance\n",
+"//calculation\n",
+"vw=p*l//pd across the wire\n",
+"x=(7*v/vw)-R//resistace of x\n",
+"//output\n",
+"printf('the resistance of x is %3.0f ohm',x)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/8-Magnetisn_and_ac_theroy.ipynb b/Concise_Physics_by_H_Matyaka/8-Magnetisn_and_ac_theroy.ipynb
new file mode 100644
index 0000000..cb654bf
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/8-Magnetisn_and_ac_theroy.ipynb
@@ -0,0 +1,573 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 8: Magnetisn and ac theroy"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.10: power_loss_ratio.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=15*10^3 //voltage\n",
+"p=80*10^3 //power\n",
+"r=430 //resistence\n",
+"v1=150*10^3//stepped value\n",
+"//calculation\n",
+"i=p/v//cable current\n",
+"i1=p/v1//stepped up cable current\n",
+"k=i*i/(i1*i1)//ratio of power loss\n",
+"//output\n",
+"printf('the ratio of power loss is %d',k)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.11: secondary_power_output.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"ep=150*10^3 //electric energy to primary\n",
+"e=0.69 //efficieny\n",
+"t=70 //time\n",
+"//calculation\n",
+"es=e*ep//transformer equation\n",
+"ps=es/t//power\n",
+"//output\n",
+"printf('the power output is %3.3e W',ps)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.12: charge_produced.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=250 //dc voltage\n",
+"s=0.22 //length\n",
+"d=4*10^-3 //diameter\n",
+"//calculation\n",
+"q=8.9*10^-12*1*0.22*0.22*250/(4*10^-3)//for air\n",
+"q1=8.9*10^-12*6.8*0.22*0.22*250/(4*10^-3)//for material\n",
+"//output\n",
+"printf('the permittivity  for air is %3.3e C',q)\n",
+"printf('\n the relative permittivity for material is %3.3e C',q1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.13: relative_permittivity.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"d=6*10^-5\n",
+"w=0.1\n",
+"er=9.4 //relative permittivity of medium\n",
+"c=1*10^-6 //capacitance\n",
+"//calculation\n",
+"l=c*d/(8.9*10^-12*er*w)//parallel plate capacitor formula\n",
+"//output\n",
+"printf('the length of wax paper is %3.3f m',l)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.14: charge_in_capacitors.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=3 //voltage\n",
+"c1=2.5*10^-6 //capacitance\n",
+"c2=2.5*10^-6\n",
+"c3=2.5*10^-6\n",
+"//calculation\n",
+"q=v/((1/c1)+(1/c2)+(1/c3))//capacitors in series\n",
+"q1=c1*v//capacitors in parallel\n",
+"//output\n",
+"printf('the pd when capacitors are in series is %3.3e C',q)\n",
+"printf('\n the pd when capacitors are in parallel is %3.3e C',q1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.15: rms_and_peak_voltage.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=14 //voltage\n",
+"//calculation\n",
+"v0=v*sqrt(2)//rms value\n",
+"//output\n",
+"printf('rms value of ac is 14 V')\n",
+"printf('\n the peak value of ac is %3.3f V',v0)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.16: Qmax_and_rms_current.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"c=65*10^-6 //capcacitor\n",
+"v=12 //voltage\n",
+"f=90 //frequency\n",
+"//calculation\n",
+"vmax=v*sqrt(2)//peak pd\n",
+"qmax=c*vmax//from eqn Q=CV\n",
+"irms=v*2*%pi*f*c//maximum charge from capacitor reactance\n",
+"//output\n",
+"printf('the maximum charge is %3.3f A',irms)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.17: capacitance_of_C.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"r=200 //resistence\n",
+"v=14 //voltage\n",
+"vr=9//pd across each component\n",
+"f=90 //frequency\n",
+"//calculation\n",
+"c=vr/(2*%pi*f*vr*r)//capacitor connected\n",
+"//output\n",
+"printf('the capacitor connected is %3.3e F',c)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.18: rate_of_change_of_pd.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=4 //voltage\n",
+"r=200 //resistence\n",
+"c=8.8*10^-6 //capacitance\n",
+"//calculation\n",
+"x=v/(r*c)//calculating V/t\n",
+"//output\n",
+"printf('the initial rate is %3.3e Vs^-1',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.19: determine_resistance_and_capacitance.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=14 //voltage\n",
+"f=90 //frequency\n",
+"i=0.4 //current\n",
+"t=55 //phase\n",
+"//calculation\n",
+"r=v/(i*sqrt(1+tand(t)^2))// value of resistance\n",
+"l=r*tand(t)/(2*f*%pi)//value of inductance\n",
+"c=1/(4*%pi*%pi*f*f*l)//value of capacitance for resonance to occur\n",
+"//output\n",
+"printf('the value of resistance is %3.3f ohm',r)\n",
+"printf('\nthe value of inductance is %3.3f H',l)\n",
+"printf('\nthe value of capacitor is %3.3e F',c)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.1: force_on_field.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"B=4.3*10^-4//magnetic flux density\n",
+"I=6.4//current \n",
+"L=4.8//length of wire\n",
+"t=24//inclination through the field\n",
+"//calculation\n",
+"f=B*I*L//force on wire when it is perpendicular\n",
+"f1=B*I*L*sind(t)//force on wire when it is inclined at t degrees\n",
+"//output\n",
+"printf('the force on wire is %3.3f N',f)\n",
+"printf('\nthe force at an anglr 24 deg is %3.3e N',f1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.2: flux_density.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"i=3.4 //current passing\n",
+"a=0.04 //distance from centre of cconductor\n",
+"//calcution\n",
+"b=(4*%pi*10^-7*5)/(2*%pi*a)//magnetic flux density\n",
+"//output\n",
+"printf('the flux density is %3.3e T',b)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.4: permeability_of_free_space.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//INPUT DATA\n",
+"Ix=1 //current in first wire\n",
+"Iy=1 //current in second wire\n",
+"FbyL=2*10^-7 //according to the definition of ampere\n",
+"a=1 //distance between the wires\n",
+"\n",
+"\n",
+"//calculation\n",
+"\n",
+"m=(2*%pi*a*FbyL)/(Ix*Iy)\n",
+"\n",
+"\n",
+"\n",
+"//output\n",
+"printf('the permeability of free space is %3.3e H/m  ',m)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.5: faraday_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"n=10 //number of rounds\n",
+"B=2*10^-2 //flux density\n",
+"a=5*10^-4 //areaof cross section\n",
+"t=10//time\n",
+"//calculation\n",
+"c=n*B*a //change in flux\n",
+"emf=c/t //induced emf\n",
+"//output\n",
+"printf('the flux changed is  %3.3e Wb ',c)\n",
+"printf('\n the induced emf is %3.3e V',emf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.6: moment_of_couple.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"N=250 //number of turns\n",
+"B=8.6*10^-4 //flux density\n",
+"I=5 //current\n",
+"A=16*10^-4//area\n",
+"t=35\n",
+"//calculation\n",
+"c=B*I*A*N*sind(t)//moment of couple\n",
+"x=c/(B*I*2*A*N)//doubling the area\n",
+"y=asind(x)\n",
+"//output\n",
+"printf('the moment of couple is %3.3e Nm',c)\n",
+"printf('\n the new angle produced  is %3.3f deg',y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.7: maximum_emf_power.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"a=20*10^-4 //area\n",
+"n=900 //number of turns\n",
+"b=5*10^-2 //flux density\n",
+"i=4.5 //current\n",
+"//calculation\n",
+"e=b*a*n*2*%pi*30//emf induced\n",
+"p=e*i//power output\n",
+"//output\n",
+"printf('the emf induced is %3.3f V',e)\n",
+"printf('\n the power output is %3.3f W',p)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.8: pd_across_motor.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"R=68 //resistence\n",
+"i=4.5 //current\n",
+"e=17 //emf\n",
+"//calculation\n",
+"v=(i*R)+e//supply pd\n",
+"//output\n",
+"printf('the supply of pd across motor is %3.0f V',v)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.9: transformer_equation.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"ns=330 //number of turns of secondary\n",
+"np=450 //number of turns in primary\n",
+"e=0.65 //efficiency\n",
+"vp=240 //ac supply of primary\n",
+"//calculation\n",
+"vs=e*(vp*ns)/np//transformer equation\n",
+"//output\n",
+"printf('the pd across secondary is %3.0f V',vs)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Concise_Physics_by_H_Matyaka/9-The_Atom.ipynb b/Concise_Physics_by_H_Matyaka/9-The_Atom.ipynb
new file mode 100644
index 0000000..a0e207a
--- /dev/null
+++ b/Concise_Physics_by_H_Matyaka/9-The_Atom.ipynb
@@ -0,0 +1,411 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 9: The Atom"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.10: count_rate_determination.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"w=0.004//weight of manganese\n",
+"a=6*10^23\n",
+"t=303*24*3600//half time\n",
+"//calculation\n",
+"N=w*a/0.054//number of moles\n",
+"x=0.693*N/(303*24*3600)//count rate from decay law\n",
+"//output\n",
+"printf('the count rate is %3.3e counts per second',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.11: determination_of_attributes.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=400//pd\n",
+"d=4*10^-3 //distance of seperation\n",
+"B=0.52//flux density\n",
+"na=6*10^23//avagadro number\n",
+"//calcuation\n",
+"E=v/d//electric field strength\n",
+"v1=E/B// speed of ions \n",
+"m=24*10^-3/na//mass of each ion\n",
+"ke=m*v1*v1/2//kinetic energy \n",
+"W=1.6*10^-19*1\n",
+"KE=ke/W//kinetic energy in electron volts\n",
+"//output\n",
+"printf('the electric field strength is %3.3e Vm^-1',E)\n",
+"printf('\n the speed of ions is %3.3e m/s',v1)\n",
+"printf('\n the kinetic energy is %3.3e J',ke)\n",
+"printf('\n the kinetic energy in electron volts is %3.3f ev',KE)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.12: velocity_selectio.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=400//pd\n",
+"d=4*10^-3 //distance of seperation\n",
+"B=0.52//flux density\n",
+"na=6*10^23//avagadro number\n",
+"//calculation\n",
+"x=2*1.6*10^-19/(4*10^-26)//specific charge of ions\n",
+"r=1*10^5/(8*10^6*B*B)// path radius\n",
+"//output\n",
+"printf('the specific charge of ions is %3.0e C/kg',x)\n",
+"printf('\n the path radius is %3.3e m',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.1: electric_field_effect.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=400 //voltage\n",
+"d=0.18 //distance of screen from centre\n",
+"e=1.6*10^-19 //electronic charge\n",
+"m=9.1*10^-31 //mass\n",
+"l=0.03 //length of parallel plates\n",
+"s=0.01 //air gap\n",
+"//calculation\n",
+"w=e*v//work done\n",
+"v1=sqrt(2*e*v/m)//speed of electron \n",
+"e1=v/s//electric field strength\n",
+"d1=d*6*10^3*l/(2*v)//vertical displacement\n",
+"//output\n",
+"printf('the work done is %3.3e J',w)\n",
+"printf('\n the speed of electron is %3.3e ms^-1',v1)\n",
+"printf('\n the displacement is %3.3f m',d1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.2: Millikan_experiment.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"v=5.7*10^-4 //velocity\n",
+"ro=830 //density\n",
+"d=4*10^-3\n",
+"V=3.2*10^3 //pd\n",
+"g=9.8 //acceleration due to gravity\n",
+"k=4.2*10^-4 //resistive force of air\n",
+"//calculation\n",
+"r=sqrt(3*k*v/(4*%pi*ro*g))//equating the forces on drop\n",
+"q=4*%pi*r^3*ro*g/(3*V/d)//electric firld between plates\n",
+"//output\n",
+"printf('the radius of oil drop is %3.3e m',r)\n",
+"printf('\n the value of electric firld between plates is %3.3e C',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.3: Stephan_Boltzmann_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"sig=6//stephans constant\n",
+"//calculation\n",
+"x=3^4*6*2^2/6//ratio of rate of emission \n",
+"//output\n",
+"printf('the ratio of rate of emission is %d and hence larger cube emits faster than smaller',x)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.4: working_temperature.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"p=900 //power\n",
+"d=4*10^-3 //diameter\n",
+"l=0.87//length\n",
+"sig=5.7*10^-8 //stephans constant\n",
+"//calculation\n",
+"t=(p/(%pi*d*l*sig))^0.25//temperature\n",
+"//output\n",
+"printf('the working temperature is %d K',t)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.5: stephan_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"e1=350//heat per second\n",
+"t=7+273 //teperature\n",
+"sig=5.7*10^-8//stephans constant\n",
+"//calculation\n",
+"e2=e1*4//stephans law\n",
+"E=sig*(t^4-t^4)//stephans law\n",
+"//output\n",
+"printf('the rate of emission is %3.3f W',e2)\n",
+"printf('\nthe rate of emission when outer temperature is increased is %d W',E)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.6: incereased_temperature_effect.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"t1=280\n",
+"t2=290//temperature of surroundings\n",
+"sig=5.7*10^-8 //stephans constant\n",
+"//calculation\n",
+"e3=sig*(t1^4-t2^4)//stephans law\n",
+"e1=6.2*10^9*sig \n",
+"e3=0.15*e1\n",
+"//output\n",
+"printf('the absorbing rate is %d W',e3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.7: plancks_theory.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"c=3*10^8 //velocity of speed\n",
+"w=5.1*10^-7 //wavelength of green light\n",
+"w1=0.7 //wavelength of radio waves\n",
+"w2=1.3*10^-13 //wavelength of gamma\n",
+"h=6.6*10^-34\n",
+"//calculation\n",
+"e1=h*c/w//plancks theory for greeen light\n",
+"e2=h*c/w1//plancks theory for radio waves\n",
+"e3=h*c/w2//plancks theory for gamma waves\n",
+"//output\n",
+"printf('energy carried by green light is %3.3e J',e1)\n",
+"printf('\nenergy carried by radio waves is %3.3e J',e2)\n",
+"printf('\nenergy carried by gamma waves is %3.3e J',e3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.8: quantities_of_metal.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"c=3*10^8//speed of light\n",
+"m=9.1*10^-31//mass of electron\n",
+"tw=5.12*10^-7//threshhold wavelength\n",
+"w1=4.52*10^-8 //radiation wavelength\n",
+"h=6.6*10^-34//stephans constant\n",
+"//calculation\n",
+"f0=c/tw//threshhold frequency\n",
+"w=h*f0//work function\n",
+"a=h*c/w1//einsteins photo electric equation\n",
+"v=sqrt((2*(a-w))/m)//photoelectric energy \n",
+"emax=0.5*m*v*v\n",
+"//output\n",
+"printf('threshhold frequency is %3.3e Hz',f0)\n",
+"printf('\n the work function is %3.3e J',w)\n",
+"printf('\n the maximum photoelectric speed is %3.3e ms^-1',v)\n",
+"printf('\n the maximum photoelectric energy is %3.3e J',emax)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.9: decay_law.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//input\n",
+"t=2.14*10^6*365*24*60*60//half time\n",
+"//calculation\n",
+"l=0.693/t//decay constant\n",
+"t1=1.1097/l//decay law\n",
+"t2=t1/(365*60*60*24)//time in yrs\n",
+"//output\n",
+"printf('time taken is %3.3e yrs',t2)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}