Numerical Solution of a non-linear Volterra Integrodifferential Equation via Runge-Kutta-Verner Method

Similar documents
Approximation of Functions Belonging to. Lipschitz Class by Triangular Matrix Method. of Fourier Series

LINEAR 2 nd ORDER DIFFERENTIAL EQUATIONS WITH CONSTANT COEFFICIENTS

1. Introduction and notations.

EE Control Systems LECTURE 11

Right Angle Trigonometry

A L A BA M A L A W R E V IE W

Advanced Engineering Mathematics, K.A. Stroud, Dexter J. Booth Engineering Mathematics, H.K. Dass Higher Engineering Mathematics, Dr. B.S.

On the Existence and uniqueness for solution of system Fractional Differential Equations

T h e C S E T I P r o j e c t

Erlkönig. t t.! t t. t t t tj "tt. tj t tj ttt!t t. e t Jt e t t t e t Jt

EEE 303: Signals and Linear Systems

x, x, e are not periodic. Properties of periodic function: 1. For any integer n,

FL/VAL ~RA1::1. Professor INTERVI of. Professor It Fr recru. sor Social,, first of all, was. Sys SDC? Yes, as a. was a. assumee.

The Heat and Mass Transfer Modeling with Time Delay

Pupil / Class Record We can assume a word has been learned when it has been either tested or used correctly at least three times.

(A) 1 (B) 1 + (sin 1) (C) 1 (sin 1) (D) (sin 1) 1 (C) and g be the inverse of f. Then the value of g'(0) is. (C) a. dx (a > 0) is

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 4, Issue 1, July 2014

Approximate Integration. Left and Right Endpoint Rules. Midpoint Rule = 2. Riemann sum (approximation to the integral) Left endpoint approximation

Opening. Monster Guard. Grades 1-3. Teacher s Guide

NEW FLOODWAY (CLOMR) TE TE PIN: GREENS OF ROCK HILL, LLC DB: 12209, PG: ' S67 46'18"E APPROX. FLOODWAY NEW BASE FLOOD (CLOMR)

Face Detection and Recognition. Linear Algebra and Face Recognition. Face Recognition. Face Recognition. Dimension reduction

DIFFERENCE EQUATIONS

Jonathan Turner Exam 2-10/28/03

P a g e 3 6 of R e p o r t P B 4 / 0 9

AE57/AC51/AT57 SIGNALS AND SYSTEMS DECEMBER 2012

Trigonometric Formula

Dec. 3rd Fall 2012 Dec. 31st Dec. 16th UVC International Jan 6th 2013 Dec. 22nd-Jan 6th VDP Cancun News

Gavilan JCCD Trustee Areas Plan Adopted October 13, 2015

1- I. M. ALGHROUZ: A New Approach To Fractional Derivatives, J. AOU, V. 10, (2007), pp

ENJOY ALL OF YOUR SWEET MOMENTS NATURALLY

Quantum Properties of Idealized GW Detector

176 5 t h Fl oo r. 337 P o ly me r Ma te ri al s

Mathcad Lecture #4 In-class Worksheet Vectors and Matrices 1 (Basics)

P a g e 5 1 of R e p o r t P B 4 / 0 9

TABLE 3 - HOUSING MATERIAL TABLE 5 - MICRON CODE TABLE 1 - PRESSURE TABLE 2 - PORT SIZE, TYPE TABLE 4 - ELEMENT CODE TABLE 8 - ACCESSORIES

FOURIER ANALYSIS Signals and System Analysis

Chapter4 Time Domain Analysis of Control System

Global Solutions of the SKT Model in Population Dynamics

Single Correct Type. cos z + k, then the value of k equals. dx = 2 dz. (a) 1 (b) 0 (c)1 (d) 2 (code-v2t3paq10) l (c) ( l ) x.

Chapter 5 Transient Analysis

OH BOY! Story. N a r r a t iv e a n d o bj e c t s th ea t e r Fo r a l l a g e s, fr o m th e a ge of 9

Special Curves of 4D Galilean Space

CONSTACYCLIC CODES OF LENGTH OVER A FINITE FIELD

Using the Printable Sticker Function. Using the Edit Screen. Computer. Tablet. ScanNCutCanvas

Chapter 3 Fourier Series Representation of Periodic Signals

CSE 373: More on graphs; DFS and BFS. Michael Lee Wednesday, Feb 14, 2018

Root behavior in fall and spring planted roses...

EXERCISE - 01 CHECK YOUR GRASP

BMM3553 Mechanical Vibrations

1 Introduction to Modulo 7 Arithmetic

UNIT I FOURIER SERIES T

Why the Junction Tree Algorithm? The Junction Tree Algorithm. Clique Potential Representation. Overview. Chris Williams 1.

Maturity in Christ. The. Rock

82A Engineering Mathematics

Thermal Stresses of Semi-Infinite Annular Beam: Direct Problem

I M P O R T A N T S A F E T Y I N S T R U C T I O N S W h e n u s i n g t h i s e l e c t r o n i c d e v i c e, b a s i c p r e c a u t i o n s s h o

MEDWAY SPORTS DEVELOPMENT

, R we have. x x. ) 1 x. R and is a positive bounded. det. International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:10 No:06 11

Classical Theory of Fourier Series : Demystified and Generalised VIVEK V. RANE. The Institute of Science, 15, Madam Cama Road, Mumbai

BULLETIN THE BULLETIN OCTOBER VICTOR VALLEY NEWSLETTER OF MINERAL CLUB GEM AND THE PH. (760) TOR VALLEY GEM & MINERAL

Madad Khan, Saima Anis and Faisal Department of Mathematics, COMSATS Institute of Information Technology, Abbottabad, Pakistan

Overview. Splay trees. Balanced binary search trees. Inge Li Gørtz. Self-adjusting BST (Sleator-Tarjan 1983).

Q.28 Q.29 Q.30. Q.31 Evaluate: ( log x ) Q.32 Evaluate: ( ) Q.33. Q.34 Evaluate: Q.35 Q.36 Q.37 Q.38 Q.39 Q.40 Q.41 Q.42. Q.43 Evaluate : ( x 2) Q.

DERIVING THE DEMAND CURVE ASSUMING THAT THE MARGINAL UTILITY FUNCTIONS ARE LINEAR

A Study of the Solutions of the Lotka Volterra. Prey Predator System Using Perturbation. Technique

( A) ( B) ( C) ( D) ( E)

, each of which is a tree, and whose roots r 1. , respectively, are children of r. Data Structures & File Management

SHINGLETON FOREST AREA Stand Level Information Compartment: 44 Entry Year: 2009

Preliminary Concept 3

BENEFITS OF COMPLETING COLLEGE Lesson Plan #3

THIS PAGE DECLASSIFIED IAW EO 12958

The Licking County Health Department 675 Price Rd., Newark, OH (740)

Inverse Fourier Transform. Properties of Continuous time Fourier Transform. Review. Linearity. Reading Assignment Oppenheim Sec pp.289.

page 11 equation (1.2-10c), break the bar over the right side in the middle

Extension Formulas of Lauricella s Functions by Applications of Dixon s Summation Theorem

Signals & Systems - Chapter 3

3+<6,&6([DP. September 29, SID (last 5 digits): --

CSE 421 Algorithms. Warmup. Dijkstra s Algorithm. Single Source Shortest Path Problem. Construct Shortest Path Tree from s

Module graph.py. 1 Introduction. 2 Graph basics. 3 Module graph.py. 3.1 Objects. CS 231 Naomi Nishimura

ME 522 PRINCIPLES OF ROBOTICS. FIRST MIDTERM EXAMINATION April 19, M. Kemal Özgören

Boyce/DiPrima/Meade 11 th ed, Ch 4.1: Higher Order Linear ODEs: General Theory

Chapter 7 INTEGRAL EQUATIONS

NEWBERRY FOREST MGT UNIT Stand Level Information Compartment: 10 Entry Year: 2001

". :'=: "t',.4 :; :::-':7'- --,r. "c:"" --; : I :. \ 1 :;,'I ~,:-._._'.:.:1... ~~ \..,i ... ~.. ~--~ ( L ;...3L-. ' f.':... I. -.1;':'.

Air Filter 90-AF30 to 90-AF60

CS 688 Pattern Recognition. Linear Models for Classification

ECE COMBINATIONAL BUILDING BLOCKS - INVEST 13 DECODERS AND ENCODERS

Handout on. Crystal Symmetries and Energy Bands

Beechwood Music Department Staff

Lecture 21 : Graphene Bandstructure

Chapter 1 Fundamentals in Elasticity

Approximately Inner Two-parameter C0

c. What is the average rate of change of f on the interval [, ]? Answer: d. What is a local minimum value of f? Answer: 5 e. On what interval(s) is f

THE MIDWAY & GAMES GRADE 6 STEM STEP BY STEP POTENTIAL & KINETIC ENERGY MOVE THE CROWDS

Software Process Models there are many process model s in th e li t e ra t u re, s om e a r e prescriptions and some are descriptions you need to mode

THIS PAGE DECLASSIFIED IAW EO 12958

1 Finite Automata and Regular Expressions

HPLC/UHPLC Columns. Limit 5 columns per company. Excludes guard columns. Offer applies to columns with 4.6 mm ID. Specialty LC Columns

Introduction to Laplace Transforms October 25, 2017

Lecture 4: Laplace Transforms

Transcription:

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 ISSN -33 Nmril Solio of o-lir Volrr Igroiffril Eqio vi Rg-K-Vrr Mho Ali Filiz * * Dprm of Mhmis A Mrs Uivrsiy 9 AYDIN-TURKEY Asr- I his ppr highr-orr mril solio of o-lir Volrr igro-iffril qio is isss. Exmpl of his qsio hs solv mrilly sig h Rg-K-Vrr mho for Oriry Diffril Eqio (ODE) pr Nwo-Cos formls for igrl prs. Ix Trms- A highr-orr ry Lgrg irpolig qrr formls Rg-K mhos o-lir Volrr igro-iffril qio. I. INTRODUCTION A fiol qio i whih h ow fio pprs i h form of i is riviv s wll s r h igrl sig is ll igro-iffril qio (s [ 4 ]). I his ppr w will osir h o-lir Volrr igroiffril qio of h form (s [ 4 7 8]) ( ) ) ( ( )( ). () ' Eqio () solv mrilly sig vrios mhos (s [ 9 ]). I his ppr ) will o h x vl of h. W shll s ( ) or o o mril solio of. Howvr i his ppr w will osr highr-orr mril mho for qio (). Si h igrl o rmi xpliily i my pproxim sig fmilir mril igrio mhos. Th Nwo-Cos igrio forml whih il h -poi los Nwo-Cos forml is ll h rpzoil rl h 3-poi rl is ow s Simpso s /3 rl h 4-poi los rl is Simpso s 3/8 rl h -poi los rl is Bool s rl (Bo s rl Wl s rl highr rls il h -poi 7-poi 8-poi r wll si hr si hy s os whih wr giv i [ 37] [4 ]. II. THE NUMERICAL INTEGRATION OF A NON-LINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATION I grl forml for h mril solio of igro-iffril qios rly po forml for h rlyig Oriry Diffril Eqio (ODE omi wih xiliry qrr rls pproximio of z ( ) : h ( ) ) ( s. () Of ors whrs w hv fi pproximios z ( ) i rms of qrr rls h rfl h rlyig ODE mho i is i priipl possil o mix mh. Th omiios of forml hos o h sis of orr of ovrg. Th firs ivolvs pig Rg-K mhos. W will rqir o pproxim igrl rm z ( ) : h ( ) ( ) ( s sl vls. Eqio () solv svrl wys. I his ppr w shll fos o highr-orr mril mho for qio (). Th igrl my pproxim sig fmilir mril igrio mhos. Th Nwo-Cos igrio formls whih il lf righ rgl rls h rpzoil rl Simpso s /3 rl Simpso s 3/8 rl r wll si hr si hy s os whih wr prviosly ll [ ]: (

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 ISSN -33 ( s h ( ) ( ) whr r h ppropri offiis for h omposi igrio shms hos. A omiio of igrio mho my s. III. NUMERICAL ROUTINE FOR NON-LINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATION Now osir h o-imsiol prolm (). I orr o solv () mrilly w prpos h s of wo mhos fmilir o mos mhmiis. W osir mhos whih pproxim h solio h iiil vl prolm (IVP) im ' ( ) ) ( s ) h = 3. whr h is h os ol sp-siz i h Exmpl 3. ( r ( s. For xmpl h xplii Elr mho pproxims h solio o Exmpl 3. ( ) ( ). s h r s s Th xplii fii iffr mho giv i [] s ppli o qio () sily x o mor r prior-orror mho. Th prior sp ss ( ( h F z( )) ) o oi whih is follow y h orror sp whih ss highr orr rpzoil mho ( )) ( h F z F z ( ). This pror is ow s moifi Elr mho (so orr Rg-K-RK) is o orr mgi mor r h h xplii Elr mho. Th forh orr lssil Rg-K mho (RK4) lso p o h mril solio of qio (). Sppig from wih sp-siz h o oi h RK4 mho s ppli o his prolm i [ ]. Th sixh orr Rg-K-Vrr mhos [3] my s o rily si h irol vlio pois r iformly sp. Cosqly h igrls rig h irmi llios o sp from o my rqir h rpzoil rl or Lgrg polyomil irpolig igrio o o-iform priio ]. [ Rg-K-Vrr mho (RKV) lso p o h mril solio of (). Sppig from h RKV mho s ppli o his prolm my wri s: h z( ) (3) wih sp-siz h o oi / h / / / / z h h/ F/ / z / 4 4/ 7 7

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 3 ISSN -33 h 3 4/ 4/ 4/ 4/ z 4 h 3 h4/ F4/ 4/ z 4/ 3 8 3 / 3 3 h 4 / 3 / 3 / 3 / 3 z h 4 h/3 F/3 / 3 z / 3 43 84 / 4 4 9 h / / / / z h h/ F/ / z / 43 4 88 8 3 h z h h F z f 83 4 433 8 4 484 / 7 8 h f f 7 / / / / z f f f h 7 h / F/ / z / 3 g 3 3 9773 394 48 387 7 43 3 3 84 73 h g g 8 z g g g h 8 h F z 3 37 4 3 3 (4) 984 8 44 3 87 3 4 4 7 438 3. () 4 44 7 9 9 I his xmpl h rpzoil rl is s o pproxim z ( ) ( s o ] ] ] [ [ / [ 4/ [ / 3 / ] [ ] [ / ] [ ] i llig 3 4 7 8 rspivly. If sir h

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 4 ISSN -33 rpzoil rl my s o ] (givs so orr ry S Tl ); h rpzoil rl Simpso s /3 rl [ (givig hir orr ry s [ ]) my s o ]. [ I orr o g highr-orr ry h igrl rm ms vl mor rly o ] ] [ / [ 4/ [ / 3] [ / ] [ ] [ / ] [ ] i llig 3 4 7 8 s show i () (7 (8 (9 ( ( () low. Th -poi x los rl is Bool's mho my vis o ] s followig: z()= )= If = z(+)= z() + h( ) + +) ) / lsif == z(+)= z(-) + h( -) +4 ) + +) ) /3 lsif ==3 z(+)= z(-) +3h ( -) +3 -) +3 ) + +) ) / 8 lsif ==4 z(+)= z(-3) +h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 lsif == z(+)= z(-4) +h (9-4) +7-3) + -)+ -) + 7 ) +9 +) ) / 88 lsif == z(+)= z(-) + h (4-)+ -4) +7-3) +7 -)+7 -) + ) +4 +) ) / 4 lsif ==7 z(+)= z(-) + 7h (7-)+377-)+ 33-4) +989-3) +989 -)+33 -) + 377 ) +7 +) ) / 78 lsif ==8 z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 lsif mo(4)== z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 lsif mo(4)== z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 lsif mo(4)== z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 lsif mo(4)==3 z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 ls z(+)= z(-3) + h (7-3) +3 -) +-) + 3 ) +7 +) ) / 4 If w irpolig o / (spil forml rqir for h firs wo sps for xmpl w s (4) ()) Lgrg s forml for pois =- - / givs 3 h h h ) ( ) ( ) ( )( ) 3( )( )( ) ( )( ) /. 3 h h h h h h 3 7 9 h [

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 ISSN -33 If w igr h xprssio w h/ w g h / ( / 3 s h 8 84 Similrly w fi = =- - 4/ 34 s h 8 8 844 4 937 49 84 4h / ( 4/ fi = =- - /3 h / 3 ( /3 4 s h fi = =- - / 8 s h 4 8 84 8 4 37 8 3 8 38 84. 38 h / ( / fi = =- - h 3 s h 8 4 4 9 4 ( fi = =- - / h / ( 3 s h 9 filly fi = =- - h 3 s h 8 8 4 4 48 9 4 (. 79 8 Thrfor h Rg-K-Vrr forml om 3 (for srig vls w s qio (4) ()) h z( ) ) () (7) (8) (9) () () () / h / / / / z h 3 z h 8 84 844 49 84 / F/ / / 4 4/ 7 7 h 3 4/ 4/ 4/ 4/ h z 34 z h 8 8 4 937 38 3 4/ F4/ 4/ 4/

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 ISSN -33 8 3 / 3 3 h 4 / 3 / 3 / 3 / 3 z h 4 z h 8 8 4 37 8 4 / 3 F/ 3 / 3 / 3 43 84 / 4 4 9 h / / / / h z 8 z h 4 84 3 8 38 84 / F/ / / 43 4 88 8 3 h h z 3 z h 8 4 4 9 4 F 83 4 7 43 8 8 4 484 f 3 / h f f 7 / / / / z f f 3 79 f 7 / / / / 9 8 48 8 h F z h 3 7 3 43 977 3 394 48 3 84 38 73 g 3 7 h g g 8 z g g 3 9 8 h F z h 8 4 4 4 (3) h sixh-orr mho 3 87 3 4 4 7 438 3 is s o sim h rror i h fifh-orr mho 4 44 7 9 9 3 37 4 3 3. 984 8 44 I Exmpl 3. w hv s Rg-K-Vrr mhos mril qrr rpzoil rl h 3-poi rl is ow s Simpso s /3 rl h 4-poi los rl is Simpso s 3/8 rl h -poi los rl is Bool s rl (Bo s rl Wl s rl highr rls il h -poi 7-poi 8-poi hir omiios. Exmpl 3.: Cosir firs orr o-lir Volrr igro-iffril qio of h form ( ) ' ( ) ) r s s ; ). (4)

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 7 ISSN -33 For lyil solio of qio (4 i (4). I qio (4 if w hoos r. w will g x solio s ( ) ). Now wriig m( ) s his is h sm s h iffril qio ) ' ( ) )( r m( ) ' ( ) r ) m( ) If m( ) s h m' ( ) ) m' '( ) '( ). Aiiolly m( ) m' () ). m' '( ) r m'( )( r m( ) m' '( ) r m'( ) m'( ) m( Eqivl o h moifi logisi qio for m () m ( ) m'( ) r m( ) m'( m ( ) m'( ) r m( ) m ( ) m' ( ) r m( ) m ( ) ( m r m( ) ( ) ' m' ( ) ( m( ) )( m( ) () r m( ) whr r roos of m ( ). If w solv qio () wih iiil oiio ( ) m w g m( ) Afr rrrgig ov solio w oi ( m( ) ( ) ( ) ) ( m( ) ) whr ( )

Iriol Jorl of Siifi Rsrh Pliios Volm 3 Iss 9 Spmr 3 8 ISSN -33. W ow h or x solio ws (). Wh m' ( ) ) ) ( ) lyil solio of qio ( ) r m (4). Hr r roos of h qio m ( ) ( ) (whr ) so h is pproh y m lrg vls of. Th xpo is fi y ( ). Tl : Nmril solios (4) for RKV mho ( r. )... 3. 4... Nmril solio h=..9943.487.93937.498.87.3 Nmril solio h=.3.99739.443.94388.487.87.9 Al Solio Rsls.9897.444.944.447.348.3 mx Error wih h=..3-4.34-3.4-4.9-4.9-9.78- Rg-K-Vrr mho (RKV) mril qrr rls rsls. REFERENCES Error wih h=.3.87-4.7-3.3-4.3-4.78-4.- [] M. Armowiz I. A. Sg (Es Igrio..4 i Hoo of Mhmil Fios wih Formls Grphs Mhmil Tls 9 h priig Nw Yor: Dovr 97 pp. 88 887. [] A. Asov Uiqss of h solio of sysms of ovolio-yp Volrr igrl qios of h firs i I: Ivrs prolms for iffril qios of h mhmil physis (Rssi) Novsiirs: A. N SSSR Siirs. Ol. Vyhil. Tsr 978 Vol pp. 34. [3] R. L. Br J. D. Firs Nmril Alysis. Nw Yor: Broos/Col Plishig Compy USA 997 h.. [4] C. T. H. Br Th Nmril Trm of Igrl Eqios. Clro Prss; Oxfor Uivrsiy Prss 977. [] C. T. H. Br G. A. Bohorov A. Filiz N. J. For C. A. H. Pl F. A. Rih A. Tg R. M. Thoms H. Ti D. R. Will Nmril Mollig y Rr Fiol Diffril Eqios Nmril Alysis Rpor Mhsr Cr for Compiol Mhmis No:33 ISS 3-7998. [] C. T. H. Br G. A. Bohorov A. Filiz N. J. For C. A. H. Pl F. A. Rih A. Tg R. M. Thoms H. Ti D. R. Will Nmril Mollig y Dly Volrr Fiol Diffril Eqios Nmril Alysis Rpor I: Compr Mhmis is Apliios-Avs & Dvlopms (994- Elis A. Lipiis (Eior LEA Plishrs Ahs Gr pp. 33-. [7] R. Bllm A Srvy of h Thory of h Boss Siliy Asympoi Bhvior of Solios of Lir No-lir iffril iffr qios Wshigo D. C. 949. [8] K. L. Coo Fiol Diffril Eqios Clos o Diffril Eqio Amr. Mh. So. 9 Vol 7 pp. 8-88. [9] A. Filiz O h solio of Volrr Lo-Volrr Typ Eqios LMS sppor O Dy Mig i Dly Diffril qio (Livrpool UK h Mrh. [] A. Filiz Nmril Solio of Som Volrr Igrl Eqios PhD Thsis Th Uivrsiy of Mhsr. [] A. Filiz Forh-Orr Ros Nmril Mho for Igro-iffril Eqios Asi Jorl of Fzzy Appli Mhmis 3 Vol I pp. 8-33. [] P. Liz Alyil Nmril Mhos for Volrr Eqios SIAM Phillphi 98. [3] C. W. Urhr Nmril Compio : Mhos Sofwr lysis. Brli: Sprigr-Vrlg 997. [4] V. Volrr Lços Sr l Thori Mhmiq l L Por L Vi. Ghir-villrs Pris 93. [] V. Volrr Thory of Fiol of Igro-Diffril Eqios. Dovr Nw Yor 99. [] V. Volrr Sll Eqzioi Igro-iffrzili Dll Tori Dll lsi Ai Dll Rl Ami i Lii 8 (99 Rpri i Vio Volrr Opr Mhmih; Mmori No Vol 3 Ami i Lii Rom 97. [7] Wolfrm MhWorl Nwo-Cos Formls vill hp://mhworl.wolfrm.om/nwo-cosformls.hml AUTHORS Firs Ahor Ahor Nm: Dr. Ali Filiz BS (Eg Uivrsiy TR MS (Th Mhsr Uivrsiy UK PhD (Th Mhsr Uivrsiy UK Assis Profssor Dprm of Mhmis A Mrs Uivrsiy 9 AYDIN- TURKEY E-mil: filiz@..r Corrspo Ahor Dr. Ali Filiz E-mil: filiz@..r. Co Nmr: (+9)- 8498 x. 4