Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica

Size: px
Start display at page:

Download "Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica"

Transcription

1 Visalisations of Gssian and Mean Cratres by Using Mathematica and webmathematica Vladimir Benić, B. sc., Sonja Gorjanc, Ph. D., Faclty of Ciil Engineering, Kačićea 6, Zagreb, Croatia Abstract. In this paper we hae gien a short oeriew on calclation of the Gassian and mean cratres of a reglar srface and in six examples we hae shown isalisations of the properties of that fnctions by sing Mathematica and color fnction He. We hae described the program webmathematica and presented one web page which is powered by this program.. Introdction Specific conditions concerning the installation of Mathematica in Croatia (the program is aailable on all niersity compters) stimlated some teachers of geometry and mathematics at the Faclties of Ciil Engineering and Geodesy to try to improe teaching and learning process by means of Mathematica and webmathematica. Within the IT project Selected Chapters of Geometry and Mathematics Treated by Means of Mathematica for Ftre Strctral Engineers we designed edcational material which enhances isally standard lectres and stimlates interactie and ttorial way of learning on the Internet. The parts of the edcational material that we hae created so far can be fond, mostly in Croatian langage, at the following address: math/ In this paper we present the part of that edcational material related to the Gassian and mean cratres of a reglar srface, which has been translated into English.. Mathematica isalisations of Gassian and mean cratres For ftre strctral engineers it is important to hae the knowledge of the Gassian and mean cratres. For example: Tensile fabric strctre (e.g. membrane roof) in a niform state of tensile prestress behaes like a soap film stretched oer a wire which is bent in a shape of a closed space cre. Soap film assmes a form which has the minimal area relatie to all other srfaces stretched oer the same wire; this srface is therefore called minimal srface. It can be shown that mean cratre anishes at each point of that srface. The project has been spported by the Ministry of Science and Technology of the Repblic of Croatia since /3.

2 .. Gassian and mean cratre of a reglar srface A reglar srface Φ R 3 is the set of points whose position ectors are the ales of one-to-one continos ector-aled fnction r : U R 3 r(, ) = (x(, ), y(, ), z(, )) () where U R is open and connected, x, y, z : U R are differentiable fnctions and (, ) U, r (, ) r (, ). () At each point of a reglar srface a niqe tangent plane and a normal ector exist. The nit normal ector n at the point with a position ector r(, ) is gien by the following formla: n = r r r r. (3) The qadratic form Ed + Fdd + Gd, where E = r r, F = r r and G = r r, is called the first fndamental form of Φ. The qadratic form Ld + Mdd + Nd, where L = n r, M = n r and N = n r, is called the second fndamental form of Φ. The Gassian cratre K and the mean cratre H of a srface Φ are fnctions K, H : U R gien by the following formlas: K = LN M EG F, EN FM + GL H = (EG F. (4) ).. Visalisations by sing Mathematica The geometrical interpretations of normal, Gassian and mean cratres of a reglar srface can be fond at the following address: math/links/sonja/gasseng/gasseng.html According to Eq. 4 we can define (in the program Mathematica) the fnctions gcratre and mcratre [3, p.394] which calclate the Gassian and mean cratres at each point of a reglar srface. These fnctions enable s to plot the graphs of the Gassian and mean cratres of reglar srfaces and to color srfaces with colors which depend of that cratres. For the following isalizations we sed the periodical Mathematica color fnction He (period ). In black-white print this color-fnction is defined in the following way : He is a color-fnction with ales on a color-spectrm. It is clear that in black-white print in this paper, which for example does not differentiate orange from ble, the distinction of isal data gien by the fnction He is decreased. Een in sch conditions the fnction He can be applied to show some properties of Gassian and mean cratres.

3 Example The parametric eqations of an ellipsoid with a center at the origin O(,, ) are: x(, ) = a cossin, y(, ) = b sin sin, z(, ) = c cos, (, ) [, ] (, ) where a, b, c R +. In Figre we show the nit sphere with radii a = b = c = (Fig. a), the oblate ellipsoid with axes lengths a = b = 4, c = (Fig. b) and the ellipsoid with axes lengths a=.5, b=4, c= (Fig. c) colored by the fnction He[3gcrate] and the graphs of their Gassian cratres (Fig. a, Fig. b, Fig. c ). Figre a Figre b Figre c K(,) K(,) K(,) Figre a Figre b Figre c 3 Example The parametric eqations of the circlar helicoid are: x(, ) = a cos, y(, ) = a cos, z(, ) = b, where a, b R \ {}. In Figre we show the circlar helicoid (a =, b=.5) oer the domain [ 4 7 3, ] [, ] colored by the fnctions He[4gcratre] (Fig. a) and He[mcratre] (Fig. c), as well as the graphs of its Gassian (Fig. b) and mean (Fig. d) cratres. It is clear from Fig. c and Fig. d that it is a minimal srface. K(,) - H(,) Figre a Figre b Figre c Figre d -

4 Example 3 The parametric eqations of the presented hyperbolic paraboloid are: x(, ) =, y(, ) =, z(, ) =, (, ) [, ] [, ]. In Figre 3 we show this paraboloid colored by the fnctions He[gcratre] (Fig. 3a) and He[mcratre] (Fig. 3c), as well as the graphs of its Gassian (Fig. 3b) and mean (Fig. 3d) cratres. K(,) H(,) Figre 3a Figre 3b Figre 3c Figre 3d Example 4 The parametric eqations of the presented 3rd degree parabolic conoid are: x(, ) =, y(, ) =, z(, ) =.5( 3 + 3), (, ) [, 5] [, ]. In Figre 4 we show this conoid colored by the fnctions He[gcratre] (Fig. 4a) and He[mcratre] (Fig. 4c), as well as the graphs of its Gassian (Fig. 4b) and mean (Fig. 4d) cratres. K(,) H(,) Figre 4b Figre 4a Figre 4c Figre 4d Example 5 The parametric eqations of the presented monkey saddle are: x(, ) =, y(, ) =, z(, ) = 3, (, ) [.8,.8] [.8,.8]. In Figre 5 we show this monkey saddle colored by the fnctions He[gcratre] (Fig. 5a) and He[mcratre] (Fig. 5c), as well as the graphs of its Gassian (Fig. 5b) and mean (Fig. 5d) cratres.

5 K(,) H(,) Figre 5a Figre 5b Figre 5c Figre 5d Example 6 The parametric eqations of the presented srface are: 5 3 x(, ) =, y(, ) =, z(, ) = sin sin, (, ) [, ] [, ]. In Figre 6 we show this srface colored by the fnctions He[gcratre] (Fig. 6a) and He[mcratre] (Fig. 6c), as well as the graphs of its Gassian (Fig. 6b) and mean (Fig. 6d) cratres. K(,) H(,) - - Figre 6a 3. Figre 6b Figre 6c Figre 6d webmathematica As it is well known web is based on client/serer architectre. When a ser wants to see some web page on his browser (Internet explorer, Opera, Mozilla etc), the browser (client) sends reqirement to the respectie serer to display that page. Then the serer sends the content of the page to the client which shows the page to the ser. The langage for designing web pages is Hyper Text Markp Langage (HTML). This langage does not spport interactie soltions of mathematically formlated problems and interactie isalisations of reslts. To be able to sole mathematical problems it is necessary to se some of specialised programs, sch as Mathematica. Bt the program Mathematica can not be directly actiated by HTML. Therefore, HTML serer has to be connected with Mathematica and it is done with the program webmathematica. In other words, webmathematica bridges web serer and the program Mathematica which enables interactie calclations and isalisations on web pages. The procedre of interactie commnication is the following: - A ser feeds data for a certain mathematical problem into his client compter. - The display of the page with reslts (nmerical, symbolic or graphical) is reqired from web serer.

6 - Web serer throgh webmathematica actiates the program Mathematica which prodces reslts and forwards them to webmathematica. Then webmathematica sends the reslts to web serer. - Web serer sends the page with the reslts to the client which displays web page to the ser. In order to be able to se webmathematica, it is necessary to install webmathematica and Mathematica on the compter with web serer. The web pages on serer hae to be written in extended HTML which is defined by the rles of webmathematica. In Fig. 7 we show the print-screen of web page powered by webmathematica which we designed within IT project mentioned in the introdction. A ser can write his inpts in white rectangles. Visalize is the command btton to start interactie commnication. The reslts (LieGraphics3D on compter) are shown in Fig. 8. Figre 7

7 Figre 8: The print-screen of web page with reslts.

8 4. Conclsion In teaching geometry Mathematica and can be sed for designing diersified and highqality edcational material. Moreoer, it is an ideal program for connecting the content of geometrical and mathematical sbjects. New technology webmathemaica opens the door to interactie compting and isalization of data directly from the ser s web proider. We hope that the presented edcation material, created for the first year stdents, wold improe stdents nderstanding of the terms related to the normal, Gassian and mean cratres of a reglar srface. References [] Benić V., Gorjanc S., 3, Compting in Geometrical Edcation at the Faclty of Ciil Engineering in Zagreb, pp.77-83, Proceeding of st Symposim on Compting in Engineering, Zagreb [] Goetz A., 97, Introdction to Differential Geometry. Addison Wesley, Reading, Massachsetts [3] Gray A., 998, Modern Differential Geometry of Cres and Srfaces with Mathematica. CRC Press, Boca Raton [4] Wickham-Jones T.,, webmathematica: User Gide. Wolfram Research Inc., pdf format [5] Wolfram S., 993, Mathematica. Addison Wesley

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica 6 th International Conference on Applied Informatics Eger, Hungary, January 27 31, 2004. Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica Vladimir Benič, Sonja Gorjanc

More information

GEOMETRICAL DESCRIPTION OF ONE SURFACE IN ECONOMY

GEOMETRICAL DESCRIPTION OF ONE SURFACE IN ECONOMY GOMTRICAL DSCRIPTION OF ON SURFAC IN CONOMY a Kaňkoá Abstract The principal object of this paper is the reglar parametric srface M in R defined by the formla The geometrical description methods we are

More information

Restricted Three-Body Problem in Different Coordinate Systems

Restricted Three-Body Problem in Different Coordinate Systems Applied Mathematics 3 949-953 http://dx.doi.org/.436/am..394 Pblished Online September (http://www.scirp.org/jornal/am) Restricted Three-Body Problem in Different Coordinate Systems II-In Sidereal Spherical

More information

We automate the bivariate change-of-variables technique for bivariate continuous random variables with

We automate the bivariate change-of-variables technique for bivariate continuous random variables with INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,

More information

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE Dynamic Systems and Applications 24 2015 349-360 TH DIFFRNTIAL OMTRY OF RULAR CURVS ON A RULAR TIM-LIK SURFAC MIN OZYILMAZ AND YUSUF YAYLI Department of Mathematics ge Uniersity Bornoa Izmir 35100 Trkey

More information

Modelling by Differential Equations from Properties of Phenomenon to its Investigation

Modelling by Differential Equations from Properties of Phenomenon to its Investigation Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University

More information

Designing of Virtual Experiments for the Physics Class

Designing of Virtual Experiments for the Physics Class Designing of Virtal Experiments for the Physics Class Marin Oprea, Cristina Miron Faclty of Physics, University of Bcharest, Bcharest-Magrele, Romania E-mail: opreamarin2007@yahoo.com Abstract Physics

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

Approximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method

Approximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method Astralian Jornal of Basic and Applied Sciences, (1): 114-14, 008 ISSN 1991-8178 Approximate Soltion for the System of Non-linear Volterra Integral Eqations of the Second Kind by sing Block-by-block Method

More information

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qalitatie featres of a discrete dynamical system of homogeneos difference

More information

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL 8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

The Open Civil Engineering Journal

The Open Civil Engineering Journal Send Orders for Reprints to reprints@benthamscience.ae 564 The Open Ciil Engineering Jornal, 16, 1, 564-57 The Open Ciil Engineering Jornal Content list aailable at: www.benthamopen.com/tociej/ DOI: 1.174/187414951611564

More information

Nonlinear parametric optimization using cylindrical algebraic decomposition

Nonlinear parametric optimization using cylindrical algebraic decomposition Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic

More information

Digital Image Processing. Lecture 8 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep.

Digital Image Processing. Lecture 8 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Digital Image Processing Lectre 8 Enhancement in the Freqenc domain B-Ali Sina Uniersit Compter Engineering Dep. Fall 009 Image Enhancement In The Freqenc Domain Otline Jean Baptiste Joseph Forier The

More information

Mean Value Formulae for Laplace and Heat Equation

Mean Value Formulae for Laplace and Heat Equation Mean Vale Formlae for Laplace and Heat Eqation Abhinav Parihar December 7, 03 Abstract Here I discss a method to constrct the mean vale theorem for the heat eqation. To constrct sch a formla ab initio,

More information

Konyalioglu, Serpil. Konyalioglu, A.Cihan. Ipek, A.Sabri. Isik, Ahmet

Konyalioglu, Serpil. Konyalioglu, A.Cihan. Ipek, A.Sabri. Isik, Ahmet The Role of Visalization Approach on Stdent s Conceptal Learning Konyaliogl, Serpil Department of Secondary Science and Mathematics Edcation, K.K. Edcation Faclty, Atatürk University, 25240- Erzrm-Trkey;

More information

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

6.4 VECTORS AND DOT PRODUCTS

6.4 VECTORS AND DOT PRODUCTS 458 Chapter 6 Additional Topics in Trigonometry 6.4 VECTORS AND DOT PRODUCTS What yo shold learn ind the dot prodct of two ectors and se the properties of the dot prodct. ind the angle between two ectors

More information

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL BBAA VI International Colloqim on: Blff Bodies Aerodynamics & Applications Milano, Italy, Jly, 2-24 28 FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL

More information

ON THE PERFORMANCE OF LOW

ON THE PERFORMANCE OF LOW Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms

More information

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model . ozdog and W. W. Olson Complex Tire-Grond Interaction Simlation: ecent eelopments Of n danced Shell Theory ased Tire odel EFEECE: ozdog. and Olson W. W. Complex Tire-Grond Interaction Simlation: ecent

More information

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields. Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

University of Bristol

University of Bristol Uniersity of Bristol AEROGUST Workshop 27 th - 28 th April 2017, Uniersity of Lierpool Presented by Robbie Cook and Chris Wales Oeriew Theory Nonlinear strctral soler copled with nsteady aerodynamics Gst

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b)

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b) . Sole Problem.-(c). Sole Problem.-(b). A two dimensional trss shown in the figre is made of alminm with Yong s modls E = 8 GPa and failre stress Y = 5 MPa. Determine the minimm cross-sectional area of

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem Otline Model Predictive Control: Crrent Stats and Ftre Challenges James B. Rawlings Department of Chemical and Biological Engineering University of Wisconsin Madison UCLA Control Symposim May, 6 Overview

More information

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L. .4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre

More information

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007 Graphs and Networks Lectre 5 PageRank Lectrer: Daniel A. Spielman September 20, 2007 5.1 Intro to PageRank PageRank, the algorithm reportedly sed by Google, assigns a nmerical rank to eery web page. More

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

FRTN10 Exercise 12. Synthesis by Convex Optimization

FRTN10 Exercise 12. Synthesis by Convex Optimization FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

III. Demonstration of a seismometer response with amplitude and phase responses at:

III. Demonstration of a seismometer response with amplitude and phase responses at: GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic

More information

Image and Multidimensional Signal Processing

Image and Multidimensional Signal Processing Image and Mltidimensional Signal Processing Professor William Hoff Dept of Electrical Engineering &Compter Science http://inside.mines.ed/~whoff/ Forier Transform Part : D discrete transforms 2 Overview

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

Curves - Foundation of Free-form Surfaces

Curves - Foundation of Free-form Surfaces Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections

More information

DIRECTLY MODELING VOICED AND UNVOICED COMPONENTS IN SPEECH WAVEFORMS BY NEURAL NETWORKS. Keiichi Tokuda Heiga Zen

DIRECTLY MODELING VOICED AND UNVOICED COMPONENTS IN SPEECH WAVEFORMS BY NEURAL NETWORKS. Keiichi Tokuda Heiga Zen DIRECTLY MODELING VOICED AND UNVOICED COMPONENTS IN SPEECH WAVEFORMS BY NEURAL NETWORKS Keiichi Tokda Heiga Zen Google Nagoya Institte of Technology, Japan ABSTRACT This paper proposes a noel acostic model

More information

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory Entrop 3, 5, -3; doi:.339/e54 Article OPEN ACCESS entrop ISSN 99-43 www.mdpi.com/jornal/entrop Deriation of D Power-Law Velocit Distribtion Using Entrop Theor Vija P. Singh,, *, stao Marini 3 and Nicola

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

Lesson 81: The Cross Product of Vectors

Lesson 81: The Cross Product of Vectors Lesson 8: The Cross Prodct of Vectors IBHL - SANTOWSKI In this lesson yo will learn how to find the cross prodct of two ectors how to find an orthogonal ector to a plane defined by two ectors how to find

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics

A State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics A State Space Based Implicit Integration Algorithm for Differential Algebraic Eqations of Mltibody Dynamics E. J. Hag, D. Negrt, M. Ianc Janary 28, 1997 To Appear Mechanics of Strctres and Machines Abstract.

More information

Optimization in Predictive Control Algorithm

Optimization in Predictive Control Algorithm Latest rends in Circits, Systems, Sinal Processin and Atomatic Control Optimization in Predictive Control Alorithm JAN ANOŠ, MAREK KUBALČÍK omas Bata University in Zlín, Faclty of Applied Informatics Nám..

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

Feature extraction: Corners and blobs

Feature extraction: Corners and blobs Featre etraction: Corners and blobs Wh etract featres? Motiation: panorama stitching We hae two images how do we combine them? Wh etract featres? Motiation: panorama stitching We hae two images how do

More information

Efficient quadratic penalization through the partial minimization technique

Efficient quadratic penalization through the partial minimization technique This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC272754474, IEEE Transactions

More information

arxiv: v1 [physics.flu-dyn] 11 Mar 2011

arxiv: v1 [physics.flu-dyn] 11 Mar 2011 arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department

More information

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1 OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps Minimizing Intra-Edge Crossings in Wiring Diagrams and Pblic Transportation Maps Marc Benkert 1, Martin Nöllenbrg 1, Takeaki Uno 2, and Alexander Wolff 1 1 Department of Compter Science, Karlsrhe Uniersity,

More information

Change of Variables. (f T) JT. f = U

Change of Variables. (f T) JT. f = U Change of Variables 4-5-8 The change of ariables formla for mltiple integrals is like -sbstittion for single-ariable integrals. I ll gie the general change of ariables formla first, and consider specific

More information

Mat 241 Homework Set 3 Due Professor David Schultz

Mat 241 Homework Set 3 Due Professor David Schultz Mat 4 Homework Set De Professor David Schltz Directions: Show all algebraic steps neatly and concisely sing proper mathematical symbolism. When graphs and technology are to be implemented, do so appropriately.

More information

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for

More information

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n 12.4 The Cross Prodct 873 12.4 The Cross Prodct In stdying lines in the plane, when we needed to describe how a line was tilting, we sed the notions of slope and angle of inclination. In space, we want

More information

Joint Transfer of Energy and Information in a Two-hop Relay Channel

Joint Transfer of Energy and Information in a Two-hop Relay Channel Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

Fast Algorithms for Restoration of Color Wireless Capsule Endoscopy Images

Fast Algorithms for Restoration of Color Wireless Capsule Endoscopy Images Fast Algorithms for Restoration of Color Wireless Capsle Endoscopy Images Haiying Li, W.-S. L, and Max Q.-H. Meng School of Control Science and Engineering, Shangdon University, Jinan, China Dept. of Electrical

More information

A Proposed Method for Reliability Analysis in Higher Dimension

A Proposed Method for Reliability Analysis in Higher Dimension A Proposed Method or Reliabilit Analsis in Higher Dimension S Kadr Abstract In this paper a new method is proposed to ealate the reliabilit o stochastic mechanical sstems This techniqe is based on the

More information

MATH2715: Statistical Methods

MATH2715: Statistical Methods MATH275: Statistical Methods Exercises III (based on lectres 5-6, work week 4, hand in lectre Mon 23 Oct) ALL qestions cont towards the continos assessment for this modle. Q. If X has a niform distribtion

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

Intelligent Positioning Plate Predictive Control and Concept of Diagnosis System Design

Intelligent Positioning Plate Predictive Control and Concept of Diagnosis System Design Intelligent Positioning Plate Predictive Control and Concept of Diagnosis System Design Matej Oravec - Anna Jadlovská Department of Cybernetics and Artificial Intelligence, Faclty of Electrical Engineering

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 2015 The athors and IOS Press. This article is pblished online with Open Access by IOS Press and distribted nder the terms of the Creative

More information

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use.

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use. MATH 2Q Spring 26 Daid Nichols Change of Variables Change of ariables in mltiple integrals is complicated, bt it can be broken down into steps as follows. The starting point is a doble integral in & y.

More information

EXPT. 5 DETERMINATION OF pk a OF AN INDICATOR USING SPECTROPHOTOMETRY

EXPT. 5 DETERMINATION OF pk a OF AN INDICATOR USING SPECTROPHOTOMETRY EXPT. 5 DETERMITIO OF pk a OF IDICTOR USIG SPECTROPHOTOMETRY Strctre 5.1 Introdction Objectives 5.2 Principle 5.3 Spectrophotometric Determination of pka Vale of Indicator 5.4 Reqirements 5.5 Soltions

More information

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined

More information

Study on the Mathematic Model of Product Modular System Orienting the Modular Design

Study on the Mathematic Model of Product Modular System Orienting the Modular Design Natre and Science, 2(, 2004, Zhong, et al, Stdy on the Mathematic Model Stdy on the Mathematic Model of Prodct Modlar Orienting the Modlar Design Shisheng Zhong 1, Jiang Li 1, Jin Li 2, Lin Lin 1 (1. College

More information

Bäcklund transformation, multiple wave solutions and lump solutions to a (3 + 1)-dimensional nonlinear evolution equation

Bäcklund transformation, multiple wave solutions and lump solutions to a (3 + 1)-dimensional nonlinear evolution equation Nonlinear Dyn 7 89:33 4 DOI.7/s7-7-38-3 ORIGINAL PAPER Bäcklnd transformation, mltiple wave soltions and lmp soltions to a 3 + -dimensional nonlinear evoltion eqation Li-Na Gao Yao-Yao Zi Y-Hang Yin Wen-Xi

More information

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Assignment Fall 2014

Assignment Fall 2014 Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions

More information

Numerical Model for Studying Cloud Formation Processes in the Tropics

Numerical Model for Studying Cloud Formation Processes in the Tropics Astralian Jornal of Basic and Applied Sciences, 5(2): 189-193, 211 ISSN 1991-8178 Nmerical Model for Stdying Clod Formation Processes in the Tropics Chantawan Noisri, Dsadee Skawat Department of Mathematics

More information

), σ is a parameter, is the Euclidean norm in R d.

), σ is a parameter, is the Euclidean norm in R d. WEIGHTED GRAPH LAPLACIAN AND IMAGE INPAINTING ZUOQIANG SHI, STANLEY OSHER, AND WEI ZHU Abstract. Inspired by the graph Laplacian and the point integral method, we introdce a novel weighted graph Laplacian

More information

STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION

STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION 8 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION Čeliovsý S.*, Hospodář P.** *CTU Prage, Faclty

More information

Invariant surfaces in H 2 R with constant (mean or Gauss) curvature

Invariant surfaces in H 2 R with constant (mean or Gauss) curvature Proceedings of XIII Fall Workshop on Geometry and Physics Mrcia, September 20 22, 2004 Pbl. de la RSME, Vol. 9 (2005), 91 103 Inariant srfaces in H 2 R with constant (mean or Gass) cratre Stefano Montaldo

More information

LIGHTWEIGHT STRUCTURES in CIVIL ENGINEERING - CONTEMPORARY PROBLEMS

LIGHTWEIGHT STRUCTURES in CIVIL ENGINEERING - CONTEMPORARY PROBLEMS ITERATIOAL SEMIAR Organized by Polish Chapter o International Association or Shell and Spatial Strctres LIGHTWEIGHT STRUCTURES in CIVIL EGIEERIG - COTEMPORARY PROBLEMS STOCHASTIC CORROSIO EFFECTS O RELIABILITY

More information

SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS

SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS THERMAL SCIENCE, Year 011, Vol. 15, No. 1, pp. 45-55 45 SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS b Khdheer S. MUSHATET Mechanical Engineering Department,

More information

arxiv: v1 [physics.flu-dyn] 4 Sep 2013

arxiv: v1 [physics.flu-dyn] 4 Sep 2013 THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]

More information

STEP Support Programme. STEP III Hyperbolic Functions: Solutions

STEP Support Programme. STEP III Hyperbolic Functions: Solutions STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

Direct linearization method for nonlinear PDE s and the related kernel RBFs

Direct linearization method for nonlinear PDE s and the related kernel RBFs Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract

More information

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper

More information