Reflection and transmission of light at a curved interface: coherent state approach. N.I. Petrov

Size: px
Start display at page:

Download "Reflection and transmission of light at a curved interface: coherent state approach. N.I. Petrov"

Transcription

1 Reflecto ad trasmsso of lght at a curved terface: coheret state approach N.I. Petrov Lea Street, 9-39, Istra, Moscow Rego 435, Russa Phase-space procedure based o coheret state represetato s proposed for vestgato of reflecto ad trasmsso of lght beams at a curved delectrc boudary. Numercal smulatos of reflecto ad trasmsso of lght at varous boudares separatg two dfferet delectrcs are carred out. Sgfcat fluece of wave-frot curvature ad polarzato of cdet beam o the reflectace ad trasmttace s show. Itroducto: Mcro-optcal elemets are wdely used moder optcal systems, such as lght homogezers, mcro-les arrays, etc. For derato of large aperture mcro-optcal systems, covetoal methods, such as physcal optcs, become overly tme-umg. Beam-mode represetatos do ot provde rght alteratve, sce the beam modes caot be tracked smply through the arbtrarly curved surfaces. Dscrete phase-space methods [, ] have bee proposed as a effcet alteratve. These methods represet felds as dscrete ad fte superpostos of elemetary Gaussa beams that ca be traced easly a complcated evromet. However, the method usually proposed to dscrete felds, the so-called Gabor represetato, has bee show umercally ustable [3, 4].

2 I ths paper, the approach based o coheret state represetato s proposed for aalyss of reflecto ad trasmsso of lght beam at a curved delectrc surface profle. Coheret states are used may dfferet areas of physcs [5-7]. I [8] the coheret state method was used for derato of oparaxal propagato ad focusg of wave beams a graded-dex medum. The term coheret states was troduced by Glauber [5, 6] for a oedmesoal steady-state quatum oscllator coecto wth problems quatum optcs. Such states were tructed ad vestgated already by Schrodger 96 [9] order to establsh a relatoshp betwee the classcal ad quatum approaches. Formulato of the problem: Cosder the curved boudary betwee two dfferet delectrc meda Fg.. For smplcty, the two-dmesoal perodcally corrugated terface y depedet wth the perod d >> λ s dered, but the exteto to the arbtrary profle of 3D case s straghtforward. The calculato procedure of the reflected ad trasmtted powers sts of the followg steps. At frst, the cdet beam feld s expaded to coheret states, represetg elemetary Gaussa beams wth axs dsplacemet ad tlt. Note, that a sese such Gaussa beams are smlar to the complex rays used for smulato of reflecto ad trasmsso of beams at a curved terface []. Coheret states CS form a full set of fuctos π d, so the arbtrary cdet feld Ex, ca be expaded to a set of CS: x, π E d x f, where x s gve by the expresso 3, d d Re d Im s the elemetary phasespace volume, f dx x E x, are the ampltudes of the expaso.

3 Coheret state ca be determed as the egefucto of the ahlato operator a : a, ˆ x kw where a + pˆ x, px. w k x A explct form of CS s gve by a Gaussa beam fucto πw 4 x exp w x + w, 3 x kw where the complex egevalues + px exp ϑ determe the tal coordates w x of the ceter of the elemetary beam ad the agle of ts clato p x s to the z- axs, s the refractve dex of medum, w s the elemetary beam wdth, k π/λ s the waveumber. Note, cotrast to Fourer-expaso, there s o requremet for orthogoalty of fuctos, owg to CS form overfull fucto system. The square of the modules f E determe the cdet beam power dstrbuto betwee the elemetary beams CS. For the 4 cdet Gaussa beam E x, / a x a π exp the ampltudes of expaso have the form w / a w / a f exp + ϑ. 4 + w / + / a w a The Gaussa elemetary beams CS pass through the terface as determed by the correspodg Fresel coeffcets that vary accordg to the agle of cdece. Usually the Fresel formulae are kow from the plae wave lmt. However, these formulae ca be used also for localzed wave beams wth the beam wasts w > λ [, ]. 3

4 4 Fally, total reflected ad trasmtted powers are determed by a summato of powers of all elemetary beams. The reflectace ad trasmttace are defed as the ratos of reflected ad trasmtted powers to the cdet power, accordgly: f d R f d P P r r, f d T f d P P t t, 5 where R ad T are the reflecto ad trasmsso coeffcets, s the cdet agle. The reflecto ad trasmsso coeffcets for TE ad TM learly polarzed cdet beams, accordgly, are gve by the expressos s s E R +, s s 4 E T + ; s s + H R, s s 4 + H T, 6 where ad are the refractve dexes of meda ad, accordgly. Results: For the surface-relef profle zsx, the ray CS wth tal coordates x, strkes the terface at the pot x, z, where the correspodg cdece agle s uquely determed. For example, the cdet agle ca be expressed as ϕ, where [ ] arcta x s ϕ, dx x ds x s s the dervatve of the surface profle fucto wth respect to the x coordate.

5 Results of smulato are preseted for dfferet surface profles, wavefrot curvature raduses ad polarzatos of cdet beam. The parameters for cdet beam, surface-relef, ad elemetary beam are the rato a > d >> w > λ, where d s the dameter of the sgle elemet corrugated surface Fg.. Fg. ad Fg. 3 show the reflectace ad trasmttace as fucto of sag h of the surface relef s x h πx / d wth the perod d5µm for TE ad TM polarzed beams wth dfferet wavefrot curvature raduses. Aalogcal depedeces are obtaed for the parabolc surface profle s x x x m s, where Rsc RL s R, R L d/ sc ad R sc are the radus ad curvature radus of the sgle elemet, x m s the ceter coordate of the sgle elemet ot show. Lesser sestvty of the reflecto ad trasmsso to the sag h ad wavefrot curvature radus R f chages s observed for parabolc surface profle. Reflectace creases ad trasmttace decreases wth the crease of sag h. The reflectace ad trasmttace are sestve to the wavefrot curvature radus ad polarzato of the cdet beam. For lower-hgher dex terface the reflectace s lower ad trasmttace s hgher for TM polarzed beam f h < d. For hgher-lower dex terface there s o evdet dfferece betwee TE ad TM polarzatos. It follows from the smulatos that the decrease of the reflectace ad crease of the trasmttace take place wth the crease of the sag h for TM polarzato owg to Brewster agle effect. The proposed method s a effcet alteratve to model aperture fuctos ad flat-topped laser beams. Ulke the superposto of off-axs Gaussa fucto compoets used [3], CS decomposto represets the superposto of learly shfted ad spatally rotated beams, formg the full set of fuctos. Trasfer-matrx method [4] ca be used to calculate the 5

6 propagato of off-axs Gaussa beams optcal systems wth tlted, dsplaced ad curved optcal elemets. Cocluso: The CS approach ca be used for smulato of testy dstrbuto of lght dffracted by the mcro-les array [5, 6]. As llustrated by smulatos, the optcal effcecy trasmttace of such systems strogly depeds o the aspect rato h/d. For h/d < the trasmttace of mcro-les arrays exceeds 9%, whch s good agreemet wth the exstg expermetal data. The proposed method s a effcet alteratve to model aperture fuctos ad flat-topped laser beams. Thus, the CS based cotuous decomposto has bee preseted for smulato of lght beam reflecto ad trasmsso at arbtrarly curved surfaces. Reflectace ad trasmttace of lght beam at dfferet surface-relefs wth depths h >> λ are vestgated. Ths method provdes effectve algorthm to model complcated structures delectrc ateas, leses, multple terfaces, etc. wth the help of combato of wave propagato ad ray-tracg procedures. Refereces. P.D. Ezger, S. Raz, ad M. Shapra, J. Opt. Soc. Am. A 3, B.Z. Steberg, E. Heyma, ad L.B. Felse, J. Opt. Soc. Am. A 8, I. Daubeches, IEEE Tras. IT-36, D. Lugara, ad C. Letrau, Electrocs Letters 34, R.J. Glauber, Phys. Rev. 3, R.J. Glauber, Phys. Rev. Lett., J.R. Klauder ad B.S. Skagerstam, Coheret states, Sgapore: World Scetfc,

7 8. N.I. Petrov, Optcs Express 9, E.S. Schrodger, Naturwss. Bd. 4, S Y.Z. Rua, ad L.B. Felse, J. Opt. Soc. Am. A 3, N.I. Petrov, Optcs Letters 9, N.I. Petrov, J. Mod. Opt. 5, A.A. Tovar, J. Opt. Soc. Am. A 8, A.A. Tovar, ad L.W. Casperso, J. Opt. Soc. Am. A, N. Ldle, J. Opt. A: Pure Appl. Opt. 4, S. 6. H. Urey, ad K.D. Powell, Appl. Opt. 44,

8 Fgure captos Fg. Geometrcal cofgurato ad coordate system for a boudary. Fg. Reflectace r sold curves ad trasmttace t dashed curves versus depth h at lower/hgher dex terface for dfferet values of cdet wavefrot curvature raduses: left - TE polarzed beam, rght TM polarzed beam; curves, R f 5 µm; curves, plae wavefrot. Fg. 3 Reflectace r sold curves ad trasmttace t dashed curves versus depth h at hgher/lower dex terface for dfferet values of cdet wavefrot curvature raduses: left - TE polarzed beam, rght TM polarzed beam; curves, R f 5 µm; curves, plae wavefrot. 8

9 Fgure z x, zsx h d x 9

10 Fgure r, t,,8,6,4,..5 TE, h, µm r, t..5 TM,,8,6,4,, h, µm

11 Fgure 3,.5. TE,.5. TM,8,8,6,6 r, t r, t,4,4,, h, µm,, h, µm

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg

More information

Ellipsometry Overview

Ellipsometry Overview llpsometry Overvew ~ R Δ p ρ = ta( Ψ) e = ~ Rs ñ(λ) = (λ) + k(λ) ε = ñ 2 p-plae s-plae p-plae plae of cdece s-plae llpsometry buldg-blocks Lght ad Polarzato Materals / Optcal Costats Iteracto of Lght wth

More information

Evolution Operators and Boundary Conditions for Propagation and Reflection Methods

Evolution Operators and Boundary Conditions for Propagation and Reflection Methods voluto Operators ad for Propagato ad Reflecto Methods Davd Yevck Departmet of Physcs Uversty of Waterloo Physcs 5/3/9 Collaborators Frak Schmdt ZIB Tlma Frese ZIB Uversty of Waterloo] atem l-refae Nortel

More information

CHAPTER VI Statistical Analysis of Experimental Data

CHAPTER VI Statistical Analysis of Experimental Data Chapter VI Statstcal Aalyss of Expermetal Data CHAPTER VI Statstcal Aalyss of Expermetal Data Measuremets do ot lead to a uque value. Ths s a result of the multtude of errors (maly radom errors) that ca

More information

Lecture 07: Poles and Zeros

Lecture 07: Poles and Zeros Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto

More information

Introduction to local (nonparametric) density estimation. methods

Introduction to local (nonparametric) density estimation. methods Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest

More information

MOLECULAR VIBRATIONS

MOLECULAR VIBRATIONS MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal

More information

On modifications of the CSC GM code (R364n)

On modifications of the CSC GM code (R364n) O modfcatos of the CSC GM code (R364) Y. Satou Jue 3, 202 Abstract The code CSC GM was modfed to hadle proto ad eutro destes separately for both target ad core ucle so as to allow drect use of proto-proto

More information

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory ROAD MAP... AE301 Aerodyamcs I UNIT C: 2-D Arfols C-1: Aerodyamcs of Arfols 1 C-2: Aerodyamcs of Arfols 2 C-3: Pael Methods C-4: Th Arfol Theory AE301 Aerodyamcs I Ut C-3: Lst of Subects Problem Solutos?

More information

AN EULER-MC LAURIN FORMULA FOR INFINITE DIMENSIONAL SPACES

AN EULER-MC LAURIN FORMULA FOR INFINITE DIMENSIONAL SPACES AN EULER-MC LAURIN FORMULA FOR INFINITE DIMENSIONAL SPACES Jose Javer Garca Moreta Graduate Studet of Physcs ( Sold State ) at UPV/EHU Address: P.O 6 890 Portugalete, Vzcaya (Spa) Phoe: (00) 3 685 77 16

More information

L5 Polynomial / Spline Curves

L5 Polynomial / Spline Curves L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a

More information

ENGI 4421 Propagation of Error Page 8-01

ENGI 4421 Propagation of Error Page 8-01 ENGI 441 Propagato of Error Page 8-01 Propagato of Error [Navd Chapter 3; ot Devore] Ay realstc measuremet procedure cotas error. Ay calculatos based o that measuremet wll therefore also cota a error.

More information

Stationary states of atoms and molecules

Stationary states of atoms and molecules Statoary states of atos ad olecules I followg wees the geeral aspects of the eergy level structure of atos ad olecules that are essetal for the terpretato ad the aalyss of spectral postos the rotatoal

More information

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices. 4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal

More information

Carbonyl Groups. University of Chemical Technology, Beijing , PR China;

Carbonyl Groups. University of Chemical Technology, Beijing , PR China; Electroc Supplemetary Materal (ESI) for Physcal Chemstry Chemcal Physcs Ths joural s The Ower Socetes 0 Supportg Iformato A Theoretcal Study of Structure-Solublty Correlatos of Carbo Doxde Polymers Cotag

More information

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract Numercal Smulatos of the Complex Moded Korteweg-de Vres Equato Thab R. Taha Computer Scece Departmet The Uversty of Georga Athes, GA 002 USA Tel 0-542-2911 e-mal thab@cs.uga.edu Abstract I ths paper mplemetatos

More information

Functions of Random Variables

Functions of Random Variables Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,

More information

Summary of the lecture in Biostatistics

Summary of the lecture in Biostatistics Summary of the lecture Bostatstcs Probablty Desty Fucto For a cotuos radom varable, a probablty desty fucto s a fucto such that: 0 dx a b) b a dx A probablty desty fucto provdes a smple descrpto of the

More information

Lecture Notes Types of economic variables

Lecture Notes Types of economic variables Lecture Notes 3 1. Types of ecoomc varables () Cotuous varable takes o a cotuum the sample space, such as all pots o a le or all real umbers Example: GDP, Polluto cocetrato, etc. () Dscrete varables fte

More information

X-ray vortices from nonlinear inverse Thomson scattering

X-ray vortices from nonlinear inverse Thomson scattering JLab semar 8/7/6 X-ray vortces from olear verse Thomso scatterg Yoshtaka Tara Natoal Isttute of Advaced Idustral Scece ad Techology (AIST) Vstg scetst: Msssspp State Uversty ad Jefferso Lab. Optcal vortex

More information

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace

More information

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions Iteratoal Joural of Computatoal Egeerg Research Vol, 0 Issue, Estmato of Stress- Stregth Relablty model usg fte mxture of expoetal dstrbutos K.Sadhya, T.S.Umamaheswar Departmet of Mathematcs, Lal Bhadur

More information

Statistical characteristics of the normalized Stokes parameters

Statistical characteristics of the normalized Stokes parameters Scece Cha Seres F: Iformato Sceces 008 SCIENCE IN CHINA PRESS Sprger www.sccha.com fo.sccha.com www.sprgerlk.com Statstcal characterstcs of the ormalzed Stokes parameters LIU Tao 1 WANG XueSog 1 & XIAO

More information

Risk management of hazardous material transportation

Risk management of hazardous material transportation Maagemet of atural Resources, Sustaable Developmet ad Ecologcal azards 393 Rs maagemet of hazardous materal trasportato J. Auguts, E. Uspuras & V. Matuzas Lthuaa Eergy Isttute, Lthuaa Abstract I recet

More information

Recovery of graded index profile by cubic spline function

Recovery of graded index profile by cubic spline function Recovery of graded dex profle by cubc sple fucto Weju Lao, Xafeg Che*, Yupg Che, Yuxg Xa ad Ygl Che Departmet of Physcs, The State Key Laboratory o Fber-Optc Local Area Commucato Networs ad Advaced Optcal

More information

Reflection of Plane Waves from a Free Surface of an Initially Stressed Transversely Isotropic Dissipative Medium

Reflection of Plane Waves from a Free Surface of an Initially Stressed Transversely Isotropic Dissipative Medium ppled Mathematcs 0 9- do:0.6/am.0.956 Publshed Ole September 0 (http://www.scrp.org/joural/am) Reflecto of Plae Waves from a Free Surface of a Itally Stressed Trasversely Isotropc Dsspatve Medum bstract

More information

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law Module : The equato of cotuty Lecture 5: Coservato of Mass for each speces & Fck s Law NPTEL, IIT Kharagpur, Prof. Sakat Chakraborty, Departmet of Chemcal Egeerg 2 Basc Deftos I Mass Trasfer, we usually

More information

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs

More information

4 Inner Product Spaces

4 Inner Product Spaces 11.MH1 LINEAR ALGEBRA Summary Notes 4 Ier Product Spaces Ier product s the abstracto to geeral vector spaces of the famlar dea of the scalar product of two vectors or 3. I what follows, keep these key

More information

TESTS BASED ON MAXIMUM LIKELIHOOD

TESTS BASED ON MAXIMUM LIKELIHOOD ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal

More information

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM Jose Javer Garca Moreta Ph. D research studet at the UPV/EHU (Uversty of Basque coutry) Departmet of Theoretcal

More information

A Method for Damping Estimation Based On Least Square Fit

A Method for Damping Estimation Based On Least Square Fit Amerca Joural of Egeerg Research (AJER) 5 Amerca Joural of Egeerg Research (AJER) e-issn: 3-847 p-issn : 3-936 Volume-4, Issue-7, pp-5-9 www.ajer.org Research Paper Ope Access A Method for Dampg Estmato

More information

AN ANALYTIC EXPRESSION FOR NEAR-FIELD ANGU- LAR GLINT PREDICTION OF RADAR SENSOR USING FAR-FIELD SCATTERING CENTERS MODELS

AN ANALYTIC EXPRESSION FOR NEAR-FIELD ANGU- LAR GLINT PREDICTION OF RADAR SENSOR USING FAR-FIELD SCATTERING CENTERS MODELS Progress I Electromagetcs Research M, Vol. 3, 5 38, 13 AN ANALYTIC EXPRESSION FOR NEAR-FIELD ANGU- LAR GLINT PREDICTION OF RADAR SENSOR USING FAR-FIELD SCATTERING CENTERS MODELS Japeg Fa, Shje Fa, Hogq

More information

KLT Tracker. Alignment. 1. Detect Harris corners in the first frame. 2. For each Harris corner compute motion between consecutive frames

KLT Tracker. Alignment. 1. Detect Harris corners in the first frame. 2. For each Harris corner compute motion between consecutive frames KLT Tracker Tracker. Detect Harrs corers the frst frame 2. For each Harrs corer compute moto betwee cosecutve frames (Algmet). 3. Lk moto vectors successve frames to get a track 4. Itroduce ew Harrs pots

More information

Beam Warming Second-Order Upwind Method

Beam Warming Second-Order Upwind Method Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet

More information

MEASURES OF DISPERSION

MEASURES OF DISPERSION MEASURES OF DISPERSION Measure of Cetral Tedecy: Measures of Cetral Tedecy ad Dsperso ) Mathematcal Average: a) Arthmetc mea (A.M.) b) Geometrc mea (G.M.) c) Harmoc mea (H.M.) ) Averages of Posto: a) Meda

More information

ESS Line Fitting

ESS Line Fitting ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here

More information

The E vs k diagrams are in general a function of the k -space direction in a crystal

The E vs k diagrams are in general a function of the k -space direction in a crystal vs dagram p m m he parameter s called the crystal mometum ad s a parameter that results from applyg Schrödger wave equato to a sgle-crystal lattce. lectros travelg dfferet drectos ecouter dfferet potetal

More information

Fractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications

Fractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications IOS Joural of Mathematcs (IOS-JM e-iss: 78-578. Volume 5, Issue 4 (Ja. - Feb. 3, PP -8 www.osrourals.org Fractoal Order Fte Dfferece Scheme For Sol Mosture Dffuso quato Ad Its Applcatos S.M.Jogdad, K.C.Takale,

More information

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific CIS 54 - Iterpolato Roger Crawfs Basc Scearo We are able to prod some fucto, but do ot kow what t really s. Ths gves us a lst of data pots: [x,f ] f(x) f f + x x + August 2, 25 OSU/CIS 54 3 Taylor s Seres

More information

Optical Design with Zemax for PhD - Basics

Optical Design with Zemax for PhD - Basics Optcal Desg wth Zema for PhD - Bascs Semar 5 : Phscal Modellg II 4-- Herbert Gross Wter term 4 www.ap.u-jea.de Prelmar Schedule No Date Subject Detaled cotet.. Repetto Correcto, hadlg, mult-cofgurato 9..

More information

x y exp λ'. x exp λ 2. x exp 1.

x y exp λ'. x exp λ 2. x exp 1. egecosmcd Egevalue-egevector of the secod dervatve operator d /d hs leads to Fourer seres (se, cose, Legedre, Bessel, Chebyshev, etc hs s a eample of a systematc way of geeratg a set of mutually orthogoal

More information

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever. 9.4 Sequeces ad Seres Pre Calculus 9.4 SEQUENCES AND SERIES Learg Targets:. Wrte the terms of a explctly defed sequece.. Wrte the terms of a recursvely defed sequece. 3. Determe whether a sequece s arthmetc,

More information

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION Malasa Joural of Mathematcal Sceces (): 95-05 (00) Fourth Order Four-Stage Dagoall Implct Ruge-Kutta Method for Lear Ordar Dfferetal Equatos Nur Izzat Che Jawas, Fudzah Ismal, Mohamed Sulema, 3 Azm Jaafar

More information

CHAPTER 3 POSTERIOR DISTRIBUTIONS

CHAPTER 3 POSTERIOR DISTRIBUTIONS CHAPTER 3 POSTERIOR DISTRIBUTIONS If scece caot measure the degree of probablt volved, so much the worse for scece. The practcal ma wll stck to hs apprecatve methods utl t does, or wll accept the results

More information

DKA method for single variable holomorphic functions

DKA method for single variable holomorphic functions DKA method for sgle varable holomorphc fuctos TOSHIAKI ITOH Itegrated Arts ad Natural Sceces The Uversty of Toushma -, Mamhosama, Toushma, 770-8502 JAPAN Abstract: - Durad-Kerer-Aberth (DKA method for

More information

Finite Difference Approximations for Fractional Reaction-Diffusion Equations and the Application In PM2.5

Finite Difference Approximations for Fractional Reaction-Diffusion Equations and the Application In PM2.5 Iteratoal Symposum o Eergy Scece ad Chemcal Egeerg (ISESCE 5) Fte Dfferece Appromatos for Fractoal Reacto-Dffuso Equatos ad the Applcato I PM5 Chagpg Xe, a, Lag L,b, Zhogzha Huag,c, Jya L,d, PegLag L,e

More information

Stochastic dynamics and the dynamic phase transition in thin ferromagnetic films

Stochastic dynamics and the dynamic phase transition in thin ferromagnetic films tochastc dyamcs ad the dyamc phase trasto th ferromagetc flms Hyubum Jag 1, Malcolm J. Grmso, ad Thomas B. Woolf 1 1 Departmet of Physology, Johs Hopks Uversty, Baltmore, MD 5 Departmet of Physcs, Uversty

More information

Descriptive Statistics

Descriptive Statistics Page Techcal Math II Descrptve Statstcs Descrptve Statstcs Descrptve statstcs s the body of methods used to represet ad summarze sets of data. A descrpto of how a set of measuremets (for eample, people

More information

Online Publication Date: 12 December, 2011 Publisher: Asian Economic and Social Society

Online Publication Date: 12 December, 2011 Publisher: Asian Economic and Social Society Ole Publcato Date: December, Publsher: Asa Ecoomc ad Socal Socety Soluto Of A System Of Two Partal Dfferetal Equatos Of The Secod Order Usg Two Seres Hayder Jabbar Abood (Departmet of Mathematcs, College

More information

Quantum Description of Optical Devices Used in Interferometry

Quantum Description of Optical Devices Used in Interferometry ADIOENGINEEING, VO. 6, NO., SEPTEME 7 Quatum Descrpto of Optcal Devces Used Iterferometr Petr KUČEA, Dept. of ado Electrocs, ro Uverst of Techolog, Purkňova 8, 6 ro, Czech epublc Optoelectrocs-aser Metrolog

More information

The solution of Euler-Bernoulli beams using variational derivative method

The solution of Euler-Bernoulli beams using variational derivative method Scetfc Research ad Essays Vol. 5(9), pp. 9-4, 4 May Avalable ole at http://www.academcjourals.org/sre ISSN 99-48 Academc Jourals Full egth Research Paper The soluto of Euler-Beroull beams usg varatoal

More information

MODELING AND OPTIMIZATION OF TWO-LAYER POROUS ASPHALT ROADS ABSTRACT 1- INTRODUCTION

MODELING AND OPTIMIZATION OF TWO-LAYER POROUS ASPHALT ROADS ABSTRACT 1- INTRODUCTION 143.DOC/1 MODELING AND OPTIMIZATION OF TWO-LAYER POROUS ASPHALT ROADS A. KUIJPERS, G. VAN BLOKLAND M+P Nose & Vbrato Cosultats, Brustesgel 3, NL-531 AD's, Hertogebosch, Netherlads Tél.: +31 73 648851 /

More information

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh

More information

1 Convergence of the Arnoldi method for eigenvalue problems

1 Convergence of the Arnoldi method for eigenvalue problems Lecture otes umercal lear algebra Arold method covergece Covergece of the Arold method for egevalue problems Recall that, uless t breaks dow, k steps of the Arold method geerates a orthogoal bass of a

More information

Multivariate Transformation of Variables and Maximum Likelihood Estimation

Multivariate Transformation of Variables and Maximum Likelihood Estimation Marquette Uversty Multvarate Trasformato of Varables ad Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Assocate Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Copyrght 03 by Marquette Uversty

More information

Analysis of Variance with Weibull Data

Analysis of Variance with Weibull Data Aalyss of Varace wth Webull Data Lahaa Watthaacheewaul Abstract I statstcal data aalyss by aalyss of varace, the usual basc assumptos are that the model s addtve ad the errors are radomly, depedetly, ad

More information

Newton s Power Flow algorithm

Newton s Power Flow algorithm Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other

More information

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity ECONOMETRIC THEORY MODULE VIII Lecture - 6 Heteroskedastcty Dr. Shalabh Departmet of Mathematcs ad Statstcs Ida Isttute of Techology Kapur . Breusch Paga test Ths test ca be appled whe the replcated data

More information

COMPROMISE HYPERSPHERE FOR STOCHASTIC DOMINANCE MODEL

COMPROMISE HYPERSPHERE FOR STOCHASTIC DOMINANCE MODEL Sebasta Starz COMPROMISE HYPERSPHERE FOR STOCHASTIC DOMINANCE MODEL Abstract The am of the work s to preset a method of rakg a fte set of dscrete radom varables. The proposed method s based o two approaches:

More information

Lecture 3 Probability review (cont d)

Lecture 3 Probability review (cont d) STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto

More information

MA/CSSE 473 Day 27. Dynamic programming

MA/CSSE 473 Day 27. Dynamic programming MA/CSSE 473 Day 7 Dyamc Programmg Bomal Coeffcets Warshall's algorthm (Optmal BSTs) Studet questos? Dyamc programmg Used for problems wth recursve solutos ad overlappg subproblems Typcally, we save (memoze)

More information

Investigating Cellular Automata

Investigating Cellular Automata Researcher: Taylor Dupuy Advsor: Aaro Wootto Semester: Fall 4 Ivestgatg Cellular Automata A Overvew of Cellular Automata: Cellular Automata are smple computer programs that geerate rows of black ad whte

More information

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler Mathematcal ad Computatoal Applcatos, Vol. 8, No. 3, pp. 293-300, 203 BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS Aysegul Ayuz Dascoglu ad Nese Isler Departmet of Mathematcs,

More information

ABOUT ONE APPROACH TO APPROXIMATION OF CONTINUOUS FUNCTION BY THREE-LAYERED NEURAL NETWORK

ABOUT ONE APPROACH TO APPROXIMATION OF CONTINUOUS FUNCTION BY THREE-LAYERED NEURAL NETWORK ABOUT ONE APPROACH TO APPROXIMATION OF CONTINUOUS FUNCTION BY THREE-LAYERED NEURAL NETWORK Ram Rzayev Cyberetc Isttute of the Natoal Scece Academy of Azerbaa Republc ramrza@yahoo.com Aygu Alasgarova Khazar

More information

Binary plasmonic waveguide arrays. Modulational instability and gap solitons.

Binary plasmonic waveguide arrays. Modulational instability and gap solitons. Bary plasmoc wavegude arrays. Modulatoal stablty ad gap soltos. Aldo Audtore Matteo Cofort Costato De Agels Dpartmeto d Igegera dell Iformaoe Uverstà degl Stud d Bresca Alejadro B. Aceves Souther Methodst

More information

Stress Wave propagation in Electro-Magneto-Elastic plate of arbitrary cross-sections

Stress Wave propagation in Electro-Magneto-Elastic plate of arbitrary cross-sections Iteratoal Joural of Scetfc & Egeerg Research Volume 3, Issue 8, August-1 1 ISSN 9-5518 Stress Wave propagato Electro-Mageto-Elastc plate of arbtrary cross-sectos P. Pousamy Departmet of Mathematcs Govermet

More information

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Joural of Appled Matematcs ad Computatoal Mecacs 04, 3(4), 7-34 FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Ata Cekot, Stasław Kukla Isttute of Matematcs, Czestocowa Uversty of Tecology Częstocowa,

More information

A Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line

A Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line HUR Techcal Report 000--9 verso.05 / Frak Borg (borgbros@ett.f) A Study of the Reproducblty of Measuremets wth HUR Leg Eteso/Curl Research Le A mportat property of measuremets s that the results should

More information

h-analogue of Fibonacci Numbers

h-analogue of Fibonacci Numbers h-aalogue of Fboacc Numbers arxv:090.0038v [math-ph 30 Sep 009 H.B. Beaoum Prce Mohammad Uversty, Al-Khobar 395, Saud Araba Abstract I ths paper, we troduce the h-aalogue of Fboacc umbers for o-commutatve

More information

Layered structures: transfer matrix formalism

Layered structures: transfer matrix formalism Layered tructure: trafer matrx formalm Iterface betwee LI meda Trafer matrx formalm Petr Kužel Practcally oly oe formula to be kow order to calculate ay tructure Applcato: Atreflectve coatg Delectrc mrror,

More information

Analysis of Lagrange Interpolation Formula

Analysis of Lagrange Interpolation Formula P IJISET - Iteratoal Joural of Iovatve Scece, Egeerg & Techology, Vol. Issue, December 4. www.jset.com ISS 348 7968 Aalyss of Lagrage Iterpolato Formula Vjay Dahya PDepartmet of MathematcsMaharaja Surajmal

More information

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat

More information

Non-degenerate Perturbation Theory

Non-degenerate Perturbation Theory No-degeerate Perturbato Theory Proble : H E ca't solve exactly. But wth H H H' H" L H E Uperturbed egevalue proble. Ca solve exactly. E Therefore, kow ad. H ' H" called perturbatos Copyrght Mchael D. Fayer,

More information

Chapter 4 Multiple Random Variables

Chapter 4 Multiple Random Variables Revew for the prevous lecture: Theorems ad Examples: How to obta the pmf (pdf) of U = g (, Y) ad V = g (, Y) Chapter 4 Multple Radom Varables Chapter 44 Herarchcal Models ad Mxture Dstrbutos Examples:

More information

6.867 Machine Learning

6.867 Machine Learning 6.867 Mache Learg Problem set Due Frday, September 9, rectato Please address all questos ad commets about ths problem set to 6.867-staff@a.mt.edu. You do ot eed to use MATLAB for ths problem set though

More information

Entropy ISSN by MDPI

Entropy ISSN by MDPI Etropy 2003, 5, 233-238 Etropy ISSN 1099-4300 2003 by MDPI www.mdp.org/etropy O the Measure Etropy of Addtve Cellular Automata Hasa Aı Arts ad Sceces Faculty, Departmet of Mathematcs, Harra Uversty; 63100,

More information

Third handout: On the Gini Index

Third handout: On the Gini Index Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The

More information

ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections

ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections ENGI 441 Jot Probablty Dstrbutos Page 7-01 Jot Probablty Dstrbutos [Navd sectos.5 ad.6; Devore sectos 5.1-5.] The jot probablty mass fucto of two dscrete radom quattes, s, P ad p x y x y The margal probablty

More information

ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS

ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS Şl uv dr g Ja-Crsta GRIGORE, Uverstatea d Pteşt, strtîrgu dvale Nr Prof uv dr g Ncolae PANDREA, Uverstatea d Pteşt, strtîrgu dvale Nr Cof

More information

Maximum Likelihood Estimation

Maximum Likelihood Estimation Marquette Uverst Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Coprght 08 b Marquette Uverst Maxmum Lkelhood Estmato We have bee sag that ~

More information

Engineering Vibration 1. Introduction

Engineering Vibration 1. Introduction Egeerg Vbrato. Itroducto he study of the moto of physcal systems resultg from the appled forces s referred to as dyamcs. Oe type of dyamcs of physcal systems s vbrato, whch the system oscllates about certa

More information

2018 TIGP: Advanced Nanotechnology (A) Part 1: Photonic Crystals and Devices

2018 TIGP: Advanced Nanotechnology (A) Part 1: Photonic Crystals and Devices 08 TIGP: Advaced Naotecholog A Part : Photoc Crstals ad Devces Lecture # M-Hsug Shh Research Ceter for Appled Sceces RCAS Academa Sca Tawa Mar. 4 08 Comparso of Quatum Mechacs ad lectrdamcs Table teded

More information

EECE 301 Signals & Systems

EECE 301 Signals & Systems EECE 01 Sgals & Systems Prof. Mark Fowler Note Set #9 Computg D-T Covoluto Readg Assgmet: Secto. of Kame ad Heck 1/ Course Flow Dagram The arrows here show coceptual flow betwee deas. Note the parallel

More information

Some Notes on the Probability Space of Statistical Surveys

Some Notes on the Probability Space of Statistical Surveys Metodološk zvezk, Vol. 7, No., 200, 7-2 ome Notes o the Probablty pace of tatstcal urveys George Petrakos Abstract Ths paper troduces a formal presetato of samplg process usg prcples ad cocepts from Probablty

More information

Module 7: Probability and Statistics

Module 7: Probability and Statistics Lecture 4: Goodess of ft tests. Itroducto Module 7: Probablty ad Statstcs I the prevous two lectures, the cocepts, steps ad applcatos of Hypotheses testg were dscussed. Hypotheses testg may be used to

More information

1 Lyapunov Stability Theory

1 Lyapunov Stability Theory Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may

More information

Lecture 2 - What are component and system reliability and how it can be improved?

Lecture 2 - What are component and system reliability and how it can be improved? Lecture 2 - What are compoet ad system relablty ad how t ca be mproved? Relablty s a measure of the qualty of the product over the log ru. The cocept of relablty s a exteded tme perod over whch the expected

More information

ONE GENERALIZED INEQUALITY FOR CONVEX FUNCTIONS ON THE TRIANGLE

ONE GENERALIZED INEQUALITY FOR CONVEX FUNCTIONS ON THE TRIANGLE Joural of Pure ad Appled Mathematcs: Advaces ad Applcatos Volume 4 Number 205 Pages 77-87 Avalable at http://scetfcadvaces.co. DOI: http://.do.org/0.8642/jpamaa_7002534 ONE GENERALIZED INEQUALITY FOR CONVEX

More information

A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures

A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures A Helmholtz eergy equato of state for calculatg the thermodyamc propertes of flud mxtures Erc W. Lemmo, Reer Tller-Roth Abstract New Approach based o hghly accurate EOS for the pure compoets combed at

More information

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model Lecture 7. Cofdece Itervals ad Hypothess Tests the Smple CLR Model I lecture 6 we troduced the Classcal Lear Regresso (CLR) model that s the radom expermet of whch the data Y,,, K, are the outcomes. The

More information

to the estimation of total sensitivity indices

to the estimation of total sensitivity indices Applcato of the cotrol o varate ate techque to the estmato of total sestvty dces S KUCHERENKO B DELPUECH Imperal College Lodo (UK) skuchereko@mperalacuk B IOOSS Electrcté de Frace (Frace) S TARANTOLA Jot

More information

Physical Optics. Lecture 12: Polarization Herbert Gross.

Physical Optics. Lecture 12: Polarization Herbert Gross. Phscal Optcs Lecture : Polarato 7-5-3 Herbert Gross www.ap.u-jea.de Phscal Optcs: Cotet No Date Subject Ref Detaled Cotet 5.4. Wave optcs G Comple felds, wave equato, k-vectors, terferece, lght propagato,

More information

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen.

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen. .5 x 54.5 a. x 7. 786 7 b. The raked observatos are: 7.4, 7.5, 7.7, 7.8, 7.9, 8.0, 8.. Sce the sample sze 7 s odd, the meda s the (+)/ 4 th raked observato, or meda 7.8 c. The cosumer would more lkely

More information

Stability For a stable numerical scheme, the errors in the initial condition will not grow unboundedly with time.

Stability For a stable numerical scheme, the errors in the initial condition will not grow unboundedly with time. .3.5. Stablty Aalyss Readg: Taehll et al. Secto 3.6. Stablty For a stable umercal scheme, the errors the tal codto wll ot grow uboudedly wth tme. I ths secto, we dscuss the methods for determg the stablty

More information

A NEW FINITE ELEMENT CONSIDERING SHEAR LAG

A NEW FINITE ELEMENT CONSIDERING SHEAR LAG Bullet of the raslvaa Uversty of Braşov CIBv 204 Vol. 7 (56) Specal Issue No. - 204 A NEW FINIE EEMEN CONSIDERING SHEAR AG A. PROIC M. VOJNIC PURCAR D. UIC Abstract: A ew model of descrbg the shear lag

More information

15. Nanoparticle Optics in the Electrostatic Limit

15. Nanoparticle Optics in the Electrostatic Limit 5 aopartcle Optcs the lectrostatc Lmt I the prevous chapter, we were restrcted to aospheres because a full soluto of Maxwell s equato was eeded Such a aalytc calculato s oly possble for very smple geometres

More information

Physical Optics. Lecture 7: Polarization Herbert Gross.

Physical Optics. Lecture 7: Polarization Herbert Gross. Phscal Optcs Lecture 7: Polarato 8--6 Herbert Gross www.ap.u-jea.de Phscal Optcs: Cotet No Date Subject Ref Detaled Cotet 8.. Wave optcs G Comple felds, wave equato, k-vectors, terferece, lght propagato,

More information

[ L] υ = (3) [ L] n. Q: What are the units of K in Eq. (3)? (Why is units placed in quotations.) What is the relationship to K in Eq. (1)?

[ L] υ = (3) [ L] n. Q: What are the units of K in Eq. (3)? (Why is units placed in quotations.) What is the relationship to K in Eq. (1)? Chem 78 Spr. M. Wes Bdg Polyomals Bdg Polyomals We ve looked at three cases of lgad bdg so far: The sgle set of depedet stes (ss[]s [ ] [ ] Multple sets of depedet stes (ms[]s, or m[]ss All or oe, or two-state

More information

Analysis of a Repairable (n-1)-out-of-n: G System with Failure and Repair Times Arbitrarily Distributed

Analysis of a Repairable (n-1)-out-of-n: G System with Failure and Repair Times Arbitrarily Distributed Amerca Joural of Mathematcs ad Statstcs. ; (: -8 DOI:.593/j.ajms.. Aalyss of a Reparable (--out-of-: G System wth Falure ad Repar Tmes Arbtrarly Dstrbuted M. Gherda, M. Boushaba, Departmet of Mathematcs,

More information