Rensselaer. Working Papers in Economics

Size: px
Start display at page:

Download "Rensselaer. Working Papers in Economics"

Transcription

1 Rensselaer Working apers in Economics Department of Economics, Rensselaer olytechnic Institute, 8th Street, Troy, NY, 8-59, US. Tel: ; Fax: ; URL: Tracking Global Factor Inputs, Factor Earnings, an Emissions ssociate with Consumption in a Worl Moeling Framework Faye Duchin Rensselaer olytechnic Institute Stephen H. Levine Tufts University Number 74 November 7

2 Tracking Global Factor Inputs, Factor Earnings, an Emissions ssociate with Consumption in a Worl Moeling Framework by Faye Duchin Department of Economics Rensselaer olytechnic Institute Stephen H. Levine Department of Civil an Environmental Engineering Tufts University bstract. This paper presents a new approach for estimating the amount of carbon emboie in a prouct consume in a given economy, taking account of where the inputs to that prouct were extracte an processe all along the supply chain. The metho is generalie to apply to all factor inputs, incluing materials an energy, as well as pollutant emissions an can track not only the flows of factors an goos as imports an exports along the global supply chain but also the payments for these inputs mae by ultimate consumers along the global value chain. The new metho makes use of absorbing Markov chains that track ownstream an upstream flows. These chains are first escribe in terms of the mathematics of a one-region input-output moel an then generalie to the global framework of a multiregional worl economy. The paper also escribes the stanar way of solving this problem, which we call the Big metho, an inicates the main avantages of the Markov chain approach, namely that it is implemente without loss of information using a more compact atabase an can aress a wier range of questions, especially ones relate to the recycling of materials. Finally, the paper iscusses the parameter requirements istinguishing this type of ex-post analysis from moel-base exploration of alternative scenarios about the future an makes the case for combining the two. C6, C67, F8, Q56, Q57

3 Tracking Global Factor Inputs, Factor Earnings, an Emissions ssociate with Consumption in a Worl Moeling Framework by Faye Duchin Department of Economics Rensselaer olytechnic Institute Stephen H. Levine Department of Civil an Environmental Engineering Tufts University. Introuction Economists often portray an economy as a circular money flow, where proucers pay salaries to workers who in turn use their incomes to purchase consumption goos an services from the proucers. Goos are assume to move aroun the circle in the opposite irection from the money flows. While the circular flow concept has been use to escribe iniviual economies, it is equally applicable to the more complex system of the worl economy. This image accurately conveys the uality between the physical flows an the money flows that constitute any economy. However, while the money flows are entirely containe within the economic system, this is not true of the physical flows. rouction requires inputs of resources that are extracte from the environment; an uring both prouction an consumption, wastes are ischarge into the environment. While some ischarges o re-enter the prouction process, reuse of proucts an recycling of materials are outsie the scope of this paper. Thus, a typical ton of iron ore will follow a path through the system, ening up emboie in a range of consumer items, incluing househol appliances an vehicles. Starting from the other en, a typical family car traces its ultimate roots to not only the iron require in the course of its fabrication but also other materials an energy as well as labor an capital an is associate also with the carbon ischarge at each stage along the way. This paper presents a metho for quantifying where that ton of ore ens up, where that car came from, an where the carbon was emitte. It also tracks the associate money flows from consumption outlays to factor payments. The approach taken escribes these paths through an economy as absorbing Markov chains an shows, both symbolically an using a numerical example, how they are implemente using inputoutput ata. This framework supplements the mathematics of input-output moeling an is a fully general treatment for associating all factor inputs an all wastes generate with specific components of final consumption, in both physical quantities an money values.

4 4 The ability to track the supply chain for consumer goos has been recognie as vital for climate change policy. party to the Kyoto rotocol now commits to targets for the carbon emitte on its territory. mechanism consiere both more effective for reucing global emissions an fairer than imposing proucer responsibility is to target emissions emboie in a country s consumption, or its consumer responsibility. To calculate the latter quantity, it must be possible to compute the amount of carbon, an where in the worl it is emitte, associate with a given final eman. Over the past several years a literature has accumulate on the use of input-output moels to calculate a country s carbon emissions uner both consumer responsibility an proucer responsibility. The carbon associate with prouction is relatively straightforwar to estimate, but quantifying the carbon emboie in consumption is more problematic because it involves ientifying the countries of origin for imports an the arrangements surrouning prouction in these countries, incluing their own imports, then the arrangements surrouning prouction in the countries from which they in turn import, an so on. Techniques have been evelope for estimating the total carbon emboiment for given countries using multi-regional input-output analyses (see Lenen, ae, an Munksgaar 4; eters an Hertwich, forthcoming. Further simplifications permit estimates of emboie carbon, incluing for imports, base on ata for only the one or several given countries (see, for example, eters an Hertwich, 6. The main objective of these stuies has been to compare, for iniviual countries, the carbon emboie in their consumption an in their prouction. This paper evelops path-base methos base on Markov chain analysis to aress a relate family of questions for all countries simultaneously an for all factor inputs an pollutant ischarges, while also accounting for irect an inirect imports. Other kins of path analysis have been applie to ecomposing input-output matrices an multipliers (Defourny an Thorbecke 984; Khan an Thorbecke 989; Sonis an Hewings, 998; eters an Hertwich 6; Lenen 7, but none has aresse the question pose here or mae use of the absorbing Markov chains. Bailey et al. (4 a,b use environ analysis, a metho rawing on input-output an to a lesser egree Markov chain analysis an previously applie to ecosystems (atten an Finn, 979, to perform path analysis, ultimately focusing on system wie measures such as average path length in the entire system, system throughput, an measures of recycling. Markov chains have been use in the analysis of energy an biomass flows in ecosystems, incluing etermining a measure of recycling (Barber, 976. We know of one pioneering effort to formally use Markov chains to stuy the recycling of materials in inustrial systems (Yamaa et al. 6; Matsuno et al. 7. The Bailey, Yamaa, an Matsuno papers are narrowly focuse on the flows of specific factors, such as aluminum, an not on their emboiment in proucts an consumption goos. They o not aress the broaer question of relating prouction to consumption as o many input-output moels. This paper, by contrast, aresses just these broaer issues an in oing so provies a brige between an input-output analysis of the entire economy an a promising approach to stuying recycling.

5 5 Within the framework of input-output moels, the Leontief inverse matrix plays a privilege role. This is the matrix (I - that is the efining feature of the basic static input-output moel, (I x = y or x = (I - y, (- where y is a vector of final eman, is a matrix of input-output coefficients, an x is the vector of output require to satisfy y. The power of the Leontief inverse is that it captures not only the irect but also the inirect input requirements. Thus, if F is the matrix of factor requirements per unit of output, the Leontief inverse makes it possible to etermine the vector of factor inputs,, require irectly an inirectly to satisfy y: = Fx = F(I - y. (- If we substitute for F the matrix C measuring pollution generate, or emissions, per unit of output, then c = Cx quantifies emissions. (Note that this is the stanar representation for emissions in an input-output framework. While in fact F an C shoul be conceptually relate, this is a workable first approximation. To simplify the notation, we use F to represent either factor inputs or emissions. The paths joining final eman an factor inputs (or emissions, represente both in physical units an in money values, are the subject of this paper. If goos prices are applie to y an factor prices to, the associate paths of money flows implie by Eq. (-, i.e., those between factor incomes (or value ae an outlays for final eman, can also be tracke. These paths can be sai to escribe the supply chain an the value chain, respectively. Some specific instances of the general question to be aresse are: How much carbon is emitte in other economies to satisfy consumer eman in the US? Which countries consumption accounts for the carbon emitte in China? How much carbon is associate with the transport of goos? How much Mile Eastern oil is associate irectly an inirectly with consumption in the US? How much of the money outlays for total consumption in the US goes to paying royalties on this oil? What are the factor input requirements for refrigerators sol in the EU, an where o the factors originate? Of the money pai for these refrigerators in the EU, where oes it en up in factor payments? How much Chinese labor is emboie in consumption of specific other countries? What portion of labor income in China is pai by consumers in these other countries? If one writes the equivalent of Eq. (- for the worl economy an analyes the chain beginning with each factor in each region an the chain ening in an element in each region s final eman, one can solve this problem for all regions simultaneously. It is stanar to o this using what we will efine as the Big metho. However, it will be seen that there are avantages for formulating the problem instea as an absorbing Markov chain.

6 6 The remainer of this paper is organie as follows. In Section, we efine ownstream flows an upstream flows in the supply chain an the value chain for each prouct, taking account of the webs of interepenence among sectors; these concepts are formalie an quantifie in subsequent sections. bsorbing Markov chains that track ownstream an upstream flows in a single economy are introuce in Section, where they are escribe in relation to the mathematics of the one-region input-output moel. Section 4 moves on to a global framework an presents the algorithms for a similar analysis of the ownstream an upstream flows in a multiregional worl economy. numerical example for tracking the supply chain an value chain, both ownstream an upstream, in the global framework is provie in Section 5. Section 6 escribes the stanar way of solving the problem, namely the Big metho, an inicates the main avantages of the Markov chain approach. The funamental istinction between this analysis of one given outcome an a moel for analying alternative scenarios is iscusse in Section 7, which escribes a global multiregional moel that can provie the inputs for the kin of analysis aresse in this paper. The final section provies conclusions, an an appenix elaborates on some of the mathematical analysis.. Downstream Flows an Upstream Flows in rouction Chains The supply chain for a particular goo can be sai to begin with the extraction of raw materials an mobiliation of labor an other factors of prouction, continues through various stages of processing an fabrication, an ens with elivery of the finishe goo to consumers. From any point mistream in the chain, prouct output flows ownstream in the irection of the consumer, while the activities consiere upstream are those proucing intermeiate inputs an ultimately the factor inputs. proucer is concerne with securing upstream suppliers an ownstream customers. From the point of view not of iniviual proucers but of the economy as a whole, the terms have a similar significance: we are concerne with tracking consumer goos upstream to the factor inputs require for their prouction an with tracking factors of prouction ownstream to the consumer goos in which they are eventually emboie. While iniviual proucers are concerne mainly with securing their own inputs, in fact factors of prouction an other inputs are also require at every intermeiate stage in the upstream prouction chain. n input-output framework is neee to escribe an quantify all of the inirect relationships that join the chains in a web-like structure. In the global context we will be intereste in tracking a goo consume in one economy upstream to the factor inputs that were utilie in proucing it in the same an other economies. We will also wish to track a factor input extracte an use in a given economy ownstream to the consumption goos in the same an other economies in which it is embee, both irectly an inirectly. Money flows in the opposite irection from the material goos. Outlays by consumers for a particular goo in one country flow upstream to the owners of the constituent

7 7 factors of prouction in the same an other economies, such that the total outlay of all consumers purchasing a given prouct is equal to the incomes receive by the owners of the factors emboie in the prouct. Moving in the other irection, the earnings of a particular factor in a given economy are pai by all ownstream users of the factor an ultimately by consumers of the goos in which the factor is embee. In the case where we are tracking not factors of prouction but wastes like carbon emissions, this waste ischarge may or may not be price (e.g., taxe. In the event, for example, of a carbon tax, the logic of the last paragraph hols. This means that the outlays of consumers can be trace upstream to factor payments an carbon taxes associate with all or part of a final bill of goos, an the earnings from a carbon tax can be trace ownstream to all the ultimate consumers whose outlays contribute to paying it. If the carbon emissions are not taxe, they simply have a ero price an can be hanle formally like any price factor or waste.. The Input-Output Moel an bsorbing Markov Chains Before consiering the multiregional global framework, we escribe the simpler context of a one-region economy with n sectors, each proucing a characteristic goo, an k factors of prouction. ll variables an coefficients are measure in relevant physical units, such as tons of steel or kwh of electricity per ton of steel. Money values are, of course, a special case of physical unit. Our examples are in general physical units, an unit prices are subsequently introuce explicitly into the analysis. In this section factors are associate with final eman first using the familiar inputoutput equations. Then we introuce an absorbing Markov chain moel an solve the same problem. Markov chains can be use to moel systems where the selection of the next state on a path is a function only of the present state an the transition probabilities associate with the branches available at the present state. In an economic system ifferent factors, an ifferent units of a given factor, take ifferent paths through the system. Thus a proportion of each factor is use irectly in sector, another proportion in sector, an so on. These proportions, of course, a to., an it is this property that a factor (or an intermeiate goo is entirely istribute among a clearly efine set of irect uses -- that allows us to use a Markov chain moel for this system. We interpret the Markovian transition probabilities as proportions of factors an goos, an it is in terms of these proportions that we will work. The section ens with a formal escription of the Markov moel an its relationship to the input-output moel in the case of a oneregion system. The Input-Output Moel Starting from Eqs. (- an (-, we istinguish the requirements for a particular factor, say i, corresponing to iniviual components of final eman by replacing the vector y by the iagonal matrix ŷ. If we call the resulting matrix Φ, then

8 8 Φ = F(I - ŷ, (- where each row of the k x n matrix Φ quantifies the istribution of one factor among the n components of final eman, while the columns escribe the requirements of all factors to satisfy each component of final eman. For a numerical example, consier an economy escribe by, F, an y: F y 8 9 Then one can etermine x, = , an Φ Reaing across the first row of the Φ matrix shows the ways in which the require 8.5 units of the first factor are istribute:. units are require to eliver the first goo to consumers, 6.7 to eliver the secon goo, an 8.5 to eliver the thir goo. The first column of the same matrix shows that elivering the first goo to consumers requires. units of the first factor an.6 units of the secon factor. Thus the Φ matrix provies the solution to our problem. The bsorbing Markov Chain Moel Next consier an analysis base on representing this economic system as a Markov chain. The system structure of our example economy is shown in Figure -, where the network noes represent the states of the system. The noes labele F refer to factors of prouction, X to outputs, an Y to final emans.

9 9 x ij = a ij x j - flow from sector i to sector j x Y y i - final eman for sector i rj = f rj x j - flow of factor r to sector j X y x x x x F F x x X X x x Y y y Figure -: Flows of Goos an Factors in a Three-Sector, Two-Factor Economy Y The objective is to follow each unit of a factor of prouction as it progresses through the system shown in Figure -, eventually ening by satisfying final eman. Because all utilie factors eventually en up embee in one or more components of final eman, the state of being part of final eman is referre to as an absorbing state an the Markov chain is referre to as an absorbing chain. (We remin the reaer that wastes an recycling are not consiere in this paper. The other states, corresponing to factors an goos, are referre to as transient states. Since the Markov chain moel quantifies each flow out of a state as a proportion of the total flow out of that state, we nee to etermine exactly what proportion of factor r is use irectly by sector j. The number of units of factor r use by sector j is rj, an the total use of factor r is r = r + r + r, so the proportion in question will be: f rj = rj / r (- This is the rj th term in the k x n matrix F. We now efine so-calle irect transition matrices of the form (a,b, where a an b are vectors representing two states, an the ij th entry of (a,b is the portion of flow out of state a i going irectly to state b j. By this efinition, (,x = F, the matrix of proportions of irect flows of factors to sectoral outputs x. Next we consier the proportional flows from the total supply of proucts, x, to final eman, y, taking into account that an intermeiate goo can flow from one sector to another any number of times before being absorbe in final eman. In the absence of

10 imports, the omestic supply is simply total omestic output, x. The proportion of supply that flows from the i th sector irectly to the i th component of final eman is simply the ratio y i /x i, an no output flows irectly from the i th sector to any other component of final eman; that is, the latter proportion is. In matrix form these proportions are elements of the n x n iagonal matrix, (x,y = x - ŷ ( = ŷ x -. Next we examine the intermeiate use of the total supply of output of each goo. The proportion of total supply, or output, flowing irectly from sector i to sector j is a ij = x ij /x i. This is the ij th term in the n x n matrix, where of course = (x,x. fter this first roun of irect flows among sectors, the intermeiate goos may be elivere either to final eman or to other sectors, an so on for all subsequent rouns. The matrix of proportional flows from x to y, incluing both irect an inirect flows, is therefore (x,y = (x,y + (x,x (x,y + (x,x (x,x (x,y + = ŷ x - + ŷ x - + ŷ x - + = (I ŷ x - = (I - ŷ x -. (- The evelopment of Eq. (- is base on the same logic as the representation of the Leontief inverse (I - as the matrix power series (I + + +, an the relation between an is spelle out below. Combining Eqs. (- an (- leas to the matrix of overall proportional flows from factors to final eman (,y = (, x(x,y = F (I - ŷ x -. (-4 The rows of this matrix istribute proportions of each factor over all components of final eman; it is reaily seen that each row total equals. since (,y is the prouct of two matrices that both have this property. Finally, to calculate the actual amounts of factors, we calculate the matrix Θ = (,y = F (I - ŷ x -. (-5 Evaluating Eq. (-5, we get.4 (,x = ,.455

11 (x, y , (,y = ,.46 an finally Θ It shoul come as no surprise that this Θ = Φ, the result we obtaine by the familiar input-output metho. It shoul be pointe out that we can express Eq. (-5 in terms of the input-output matrices an F by recogniing that f rj = f rj (x j / r an a ij = a ij (x j /x i. Therefore where X is the intersectoral flow matrix. F = - F x, an (-6a = x - x = x - X (-6b The Formal Markov Chain Moel The familiar input-output manipulations on a atabase measure in a common unit (namely money values or mass are a special case of the more general Markov chain analysis. When the flows are measure in mixe units, traitional input-output moels are not Markov chains. In this case the input-output moel can be solve provie that the matrix satisfies the Simon-Hawkins conitions, meaning that the economy oes not operate at a loss. The Markov chain moel in mixe units can be analye even for an economy that epens upon subsiies. We conclue this section by presenting a formal absorbing Markov chain analysis of this system to emonstrate its aitional properties. In Section 4, we will exploit the greater generality of the Markov metho to link factors to final eman in the more complex system of the global multiregional economy. In a Markov chain the system is represente by a number of states an the one-step transition probabilities, what we are referring to as the irect proportional flows, from

12 one state to another. These irect transition probabilities are represente for the entire system in what is calle the systemwie transition matrix, enote as M. Thus m ij is the probability that state i will transition in the next step to state j. In an absorbing Markov chain analysis, the states are partitione into absorbing states an transient states, an the transition matrix is put into canonical form by orering the states starting with the absorbing states first. Then M can be partitione as follows, I M. (-7 R Q The absorbing states have the property that once entere they can never be exite. Thus if state i is an absorbing state m ii =. (since % of the flow remains in the same state an m ij = for j i. The ientity matrix I an the matrix, the latter being rectangular, reflect these properties. The irect proportional flows among the transient states are represente by the matrix esignate Q an the irect proportional flows from transient to absorbing states by the matrix esignate R. The absorbing chain is sai to have a funamental matrix (Kemeny an Snell, 976, efine as N = I + Q + Q + = (I Q -. (-8 This matrix contains information about the average number of times a unit of a factor of prouction or a sectoral output passes through each transient state before reaching an absorbing state. Since R contains the irect proportional flows from those transient states to the absorbing states, the matrix B = NR (-9 contains the proportion of the output of each transient state that ultimately reaches each absorbing state, either irectly or inirectly. The B matrix is therefore what we sought in the analysis an the numerical example above. In our example the final emans (y correspon to the absorbing states while the transient states are the factors ( an the sectoral outputs (x. Noting that, an making use of the irect transition matrix notation, (y,y = I, the partitione M matrix is

13 y x F I (- (x,x (x,y x ( (y,y M, Therefore, F Q (- an the funamental matrix is (I (I F I N. (- The matrices in the first row an column of N correspon to the factors while the matrices in the secon row an column correspon to the sectoral outputs. Finally, y x R (- an y x (I y x (I F B. (-4

14 4 The matrix constituting the first block of B correspons to the proportional flows from the factors to final eman. We will efine this as B = F (I - x - ŷ. (-5 The secon block matrix of B correspons to the proportional flows from the sectors to the final emans an is efine as B x = (I - x - ŷ. (-6 B is the matrix of interest, an we see that in fact it is exactly what we previously labele (,y. The other component of B, B x, correspons to (x,y. In this section we have shown how the absorbing Markov chain analysis replicates the results of the more traitional input-output analysis. However, by virtue of Eq. (- being a funamental matrix of an absorbing chain, we can erive aitional results escribing the paths taken by the factors as they flow through the system. s one example, the row sums in N corresponing to the factors are the average path lengths taken by those factors before being absorbe as final eman. Similarly, the row sums corresponing to the sector outputs are the average path lengths taken by those outputs. In the form presente in Eq. (- this is an average across all the final emans, but the funamental matrix can also be moifie to escribe the average path length until absorption by a specific component of final eman, that is, a specific absorbing state. Note that while (I - has the property that its row sums are average path lengths, the same is not true for the row sums (or the column sums when flows are measure in physical units of (I The Global Supply an Value Chains as bsorbing Markov Chains In this section we generalie the concepts presente in the one-region moel of Section in orer to apply them to a worl economy consisting of m regions, each escribe in terms of n sectors an k factors of prouction. In this global framework, we explicitly escribe the role of imports an exports in each region s economy. The inter-regional transport of trae flows itself calls for factor inputs an generates emissions, an these will be explicitly aresse. Note that this transport is require per unit of trae flow between specific regions, not per unit of one region s output, so its representation requires special treatment in an input-output framework. We treat the inustries proviing international transport as sectors in each regional economy an for now represent the eman for their output as part of the final eman of importing regions. Later the treatment of this istinctive sector will be further elaborate. Figure 4- below is a network moel illustrating the flows of goos an services, measure in physical units, in a -region system.

15 5 Region Final Deman Region Factors of rouction u = y + t Region Final Deman Region Factors of rouction u = y + t Region,F e e Region,F e e e e Region,F Region Factors of rouction u = y + t Region Final Deman y consumer eman t inter-regional transport eman e import/export technical matrix T inter-regional transport matrix F factor matrix t = T e + T e t = T e + T e t = T e + T e Figure 4-: Flows of Goos an Factors in -Region Economic System Consumer eman an inter-regional transport eman are combine into one vector, u, an referre as final eman. The u s an e s are n-element vectors, one element for each sector, an the s are k-element vectors, one for each factor. The noes representing the regions are more complex than the others because they involve the economic structure characterie by the matrices of intermeiate inputs an of factors of prouction, an F, consiere in Section. Figure 4- illustrates the internal structure of the noe for region ; the noes for regions an are similarly organie. The intermeiate an total output vectors are w an x, respectively, while is the vector of total omestic supply incluing both omestic prouction an imports.

16 6, F x w + u e e e e Figure 4-: Details of Region Noe for Flows of Goos an Factors s in the previous section we will analye the system in Figure 4- as an absorbing Markov chain, with the final eman noes treate as absorbing states. Since each noe in Figure 4- represents a set of noes, one for each economic sector or factor of prouction, the flows associate with the branches are all vectors. We first evelop the irect transition matrices, (a,b, an then incorporate them into a formal absorbing chain moel. aralleling the escription in Section, we begin with the flow in region g from g to x g. The proportion of the r th factor flowing to the j th sector is grj f grj n. (4- s This is the rj th term in the k x n matrix F g = ( g,x g. The proportion of the i th sector s intermeiate output flowing irectly to the j th sector is x gij /w gi, where w gi is the total intermeiate use of the output of sector i. This is the ij th term of the n x n matrix (w g,x g = where X g is the intersectoral flow matrix for region g. grs ŵ - X g g, (4-

17 7 The proportion of x g flowing irectly to g, total omestic supply, is clearly (x g, g = I, (4- the n x n ientity matrix, since each component of x g is irectly an completely incorporate in the corresponing component of g. The same is true for (e hg, g for h g. The proportion of g flowing irectly to w g is the proportion of each unit of supply of g that is incorporate into w g, an is therefore ( g,w g = g ẑ - ŵ. (4-4 g The same relationship hols for g to e gh, h g, as well as g to u g (recalling that u g = y g + t g. Flows from region g to region h are represente as Finally, ( g, h = ( g,e gh (e gh, h ẑ - ê I = g gh = ẑ - g gh ( g,x g = ( g,w g (w g,x g ẑ - ŵ ŵ - X g = g ê. (4-5 g g = ẑ - X g g = g, (4-6 where x gij a gik. (4-7 gi is, as in Section (see Eq. (-6b, the irect istribution of output to ifferent sectors as a proportion of total omestic supply. In the multiregional case, omestic supply inclues imports as well as omestic output: x e. (4-8 g g h g hg Noes, the final emans, contain the n absorbing states. For each of these noes there is only one branch leaving the noe an it returns into the same noe. (These branches, not shown in Figure 4-, are shown in Figure 4-. The matrix of proportions of flow associate with this single branch, (u,u, is the n x n ientity matrix I.

18 8 We now have etermine all the matrices require for the absorbing Markov chain. To put the system in canonical form, the absorbing states (i.e., the final emans, are labele noes -; an the transient states are noes 4. We recast the network moel in Figure 4- as a Markov chain, shown in Figure 4-. Each region is now represente by two noes corresponing to the vectors x an, the latter incluing imports. Vector States u in Region u in Region - u in Region 4 - in Region 5 - in Region 6 - in Region (,u (u,u (x, 8 (,x 7 x in Region 8 - x in Region 9 - x in Region in Region in Region in Region (,x 5 (, (, (, ( (u,u (,, (,u (,u (u,u (, (,x (x, (,x 7 (x, 4 (,x (,x 9 6 Figure 4- bsorbing Markov Chain Moel for -Region System Showing the Flows of Goos an Factors. Each (, is a Direct roportional Flow Matrix The canonical form for the matrix of proportions of irect flows of an absorbing chain, as escribe in Section, is shown in Figure 4-4. For the multiregional system, the imension of the I matrix is (nm x nm, the matrix is (nm x m(k+n, the R matrix is (m(k+n x nm, an the Q matrix is (m(k+n x m(k+n. Recall that the R matrix escribes the irect proportions of flows from the transient states (the s, s an x s to the absorbing states (the u s, an the Q matrix escribes the irect proportions of flows between the transient states. ll row sums of M equal.: this is evient for the first nm rows, an also hols for the row sums of R plus Q, since R contains the portions of omestic supply of goos ( istribute irectly to final eman an Q contains the portions use as intermeiate inputs an the portions exporte.

19 9 M = I O R Q I irect proportional flows from between absorbing states (ientity matrix irect proportional flows from absorbing states to transient states (ero matrix R irect proportional flows from transient states to absorbing states Q irect proportional flows from transient states to transient states Figure 4-4: Canonical Form of bsorbing Chain Direct roportional Flow Matrix For our -region system the M matrix is shown in Figure 4-5. The first of these matrices is the irect representation of Figure 4- an is partitione in accorance with its canonical form into the I,, R an Q matrices, with the last three further partitione. The matrix M is sparse: a number of the partitions are matrices, an others are iagonal matrices.

20 M =, (, ( x, (,u (, (, ( x,,u (, (, ( x, (,u (, (x, (x, (x x, ( x, ( x, (,u (u,u (u,u (u = e e u e e u e e u I I I F F F I I I Figure 4 5: -Region roportional Flow Matrix in Canonical Form for the Flow of Goos an Factors

21 s in Section, we wish to etermine the absorption matrix B = NR. To exploit the wellefine structure an sparsity of M, we write: R u, Q F I e an B B B x. B Then we are able to express B as follows: B = (I Q - R (I QB = R B QB = R. Solving for B we obtain: This leas to B F B x =, B x - B =, an B B x ẑ - ê B = ẑ - û. B = F (I - ẑ - ê - ẑ - û, an B x = B = (I - ẑ - ê - ẑ - û. (4-9a (4-9b (4-9c (4-a (4-b We note in Figure 4-5 that is a block-iagonal matrix while the submatrix we have calle ẑ - ê has eroes in the iagonal blocks. If we efine à = + ẑ - ê, then à is the mn x mn matrix with the proportional istributions for regional prouction matrices own the iagonal an the information on exports in off-iagonal blocks. We can rewrite the key equations as follows: B = F (I à - ẑ - û, an B x = B = (I à - ẑ - û. (4-a (4-b These equations generalie Eq. (-5 an (-6 to the case of multiple regions that trae among themselves. To calculate the amounts, an not just the proportions, of factors associate with each component of final eman, we construct the augmente vector of multiregional factor use

22 an then, parallel to Eq. (-5, Φ = B = F (I Ã - ẑ - û. (4- Finally we return to the important fact that factors of prouction are require, an carbon an other pollutants emitte, in carrying internationally trae goos. These requirements have been accounte for but combine with consumer eman in the u vector. We can o this without istortion or loss of information since there is ero consumer eman for the output of the inter-regional transport sectors. However, as our goal is to associate all factor inputs an pollutant emissions with consumer eman for the goos-proucing sectors, it will be necessary to appropriately reallocate them to final eman for the outputs of these other sectors. This can be one by appropriately moifying the Φ matrix using a proceure that will be escribe an illustrate in the numerical example of Section Region Numerical Example In this section we emonstrate the absorbing chain approach through a -region example. Each region s economy is escribe in terms of four sectors: ( agriculture, ( manufacturing, ( oil extraction, an (4 inter-regional transportation. There are three factors of prouction: ( labor, ( crue oil, an ( lan. The illustrative ata are such that region is a stylie representation of a inustrialie economy, region of an agricultural economy, an region of a non-inustrialie economy that is well-enowe with oil. It is assume that not only the parameters but also the values of all variables are known: the implications of this assumption are the subject of Section 6. The specifics of the example follow, starting with the matrices of intermeiate an factor inputs per unit of output with all quantities, unless otherwise specifie, measure in physical units: F F F Final eman an output in each region are:

23 y 5 y 6 8 y an x 75 x 98.7 x. The vectors of factor use are: , 6, an Trae goos are transporte between regions with the requirements for transport of each unit of imports to region j from region i escribe by entries in the rows of the matrix T ij (shown below corresponing to ifferent transport sectors a single such sector in the example below. Each entry in that row (the 4 th row, measure in ton-kilometers, is the prouct of the mass of the goo transporte multiplie by the istance between the origin region an the estination region. Thus for each unit of agricultural goos importe from region to, or region to,. ton-km of inter-regional transport is require. For heavier manufacture goos, the corresponing figure is.5 ton-km T T.5..4 T T...48 T T The trae vectors quantify the flows of goos an services among the three regions: 5.5 e e e 67.8 e 8 e 7.5 e. The unit prices of goos, which reflect ifferences in transport costs, are:

24 4 p p p, an the per-unit factor prices an scarcity rents are π 5 π.5 5 π an r r. 5 r, respectively. In this case, there is a scarcity rent only on lan in region. When there are no scarcity rents, or when they are presume to be inclue in the factor prices, the π vectors alone suffice. (The reason for incluing explicit scarcity rents becomes clearer in Section 7, where the Worl Trae Moel is iscusse. The matrices B, B x, an B are compute as escribe in Section 4. The row sums of B are equal to., except for those rows corresponing to factors that were not utilie, since each row shows the proportions of flows, both irect an inirect, from a given factor in one region to each component of final eman in all regions. B The columns of B correspon to the components of the final eman vectors u, u, an u, each with four components: agriculture, manufacturing, oil extraction, an international transportation. The nine rows correspon to the components of the vectors,, an, each with three components for labor, crue oil, an lan. Thus, the first row shows that 44% of the labor in region (the sum of the first 4 entries is require to satisfy its own final eman both irectly an inirectly. That means that 56% of the region s labor is emboie in the consumption goos of the other regions. Note that in this very aggregate example, regions an extract no crue oil an the oil-rich region uses no lan for agricultural prouction, as reflecte by the rows of eroes in B. Note

25 5 also that rows 4 an 6, corresponing to an, the istribution of its utilie lan an labor, are ientical because region prouces only one output. The same is true for rows 7 an 8. Utiliing Eq. (4-, we can calculate the Φ matrix in physical quantities. Note that the row sums in B that are correspon to the failure to use specific factors an therefore get multiplie by in this operation: Φ Each row of Φ shows how a particular factor in a given region is istribute among all components of final eman in all regions, while each column shows how many units of each factor originating in all regions are require to meet final eman for a specific goo in a particular region. Naturally, the vector of row sums is exactly equal to the vector. Columns can in general not be totale since factors may be measure in ifferent units. Now we examine the istribution of income in this multiregional system. The relevant matrices are etermine by multiplying each row of B, an each row of Φ, by the region-specific factor price, incluing any scarcity rents, that is, by the relevant component of (π i + r i for region i; we call the resulting matrices B W an Φ W : B W

26 6 Each row sum of B W is equal to the earnings of one unit of that factor, an the row shows the portion of the unit earnings pai by consumers of each goo in each region. Φ w Each row total of Φ W shows the total earnings of one factor in each region, an iniviual entries in the row show the source of the payment in the consumption of a particular final goo in some region. For example, the total of row is 7.7 money units, say thousans of ollars; this means that workers in region earn $7,7, of which $75,8 comes out of consumption outlays in the same region (the sum of the first 4 figures in the row. Most worker income in region is associate with manufacturing, amounting to over $, of the total. The largest factor earnings (i.e., the highest row total are for lan in region, which takes in $75,; this is not surprisingly the region that specialies in agriculture. Total income in region is $8,6, the sum of the first row sums. Now assume that the thir factor correspons to water pollution instea of lan inputs an that the ischarge of a unit of pollutant is taxe at ifferent rates in the ifferent regions. Then we see that total worl taxes amount to $76,5 (the sum of row-sums, 6 an 9, most of which are pai by region, while these ischarges are either insignificant in quantity, or not taxe, in region. Most of the tax is pai by consumers of agricultural goos, mainly the large number of them in region. Each column total of Φ W measures the outlays for all factors associate with the consumption of a particular goo in a given region, an the sum of all column totals for a given region is the total value of consumption outlays, p i T y i. Thus consumption of agricultural goos in region amounts to $, (column sum, most of which goes to pay lan owners in region ($6, an workers in region ($5,. The sum of all row totals for a given region is the total factor income, or the (supply-sie GD for that region. The sum of all column totals for each region is the value of total consumption outlays. The latter iffers from the former by the value of net exports. For example, region has total factor earnings of $84, but nees to lay out $475, for consumption; this eficit is explaine by the fact that it is a net importer, an the value of its trae eficit accounts for the ifference. By contrast region is a net exporter. Its

27 7 factor earnings are $,5,, but its outlays for consumption are only $9,7, with its trae surplus accounting for the ifference. Finally, we return to the question of the factor resources require by inter-regional transport of trae goos an reallocate these factors to the importe goos. Reallocation of Transport to Final Deman for Goos We wish to reallocate the factors associate in each region with the international transport of its imports to the final eman in the same an other regions where those imports are absorbe. To o this, the information neee, an its sources, are as follows: (a The amounts of factors to be reallocate: in the column of Φ for each region corresponing to the transportation sector(s. (b The proportion of each region s transport eman (in ton-km associate with the import of each goo: obtaine from T ij e ij. The factors are allocate among imports in these proportions. (c The proportions among components of final eman, by region an by sector, where each region s imports en up, which is the same as where each region s total omestic supply, (output plus imports, ens up: obtaine from the rows of B. The factors associate with each import are allocate to final eman in these proportions. See the appenix for a more etaile escription of the algorithm. The process may nee to be iterate but is assure to converge. In the -region example, this yiels (after iterations the matrix we call B R, to be compare with B : B R The bsorbing Markov Chain pproach Compare to the Big pproach s shown in Section, the upstream an ownstream chains in the one-region case can be etermine using either the stanar input-output approach or the absorbing Markov

28 8 chain approach. In fact, the same is true in the case of the multiregional worl economy, for which the absorbing Markov chain approach was escribe in Sections 4 an 5. Here we escribe the stanar input-output approach to the multiregional problem. Following the input-output notation of Section, the solution to the -region, -factor, 4- goo problem coul be written as Φ = F(I Å - ŷ, where ŷ is the x iagonal matrix of final eman for the regions, F is the 9 x block-iagonal matrix with all regions F matrices (each of imension x 4 arrange own the iagonal, an, which we will call the Big matrix, is efine as follows:. (6- In the Big matrix, each region s matrix, i, is ivie column-wise into parts, one corresponing to omestically supplie inputs per unit of its output ( ii an the others to imports from each other region, j, per unit of region i s output ( ji. Thus, taking region as an example, = + +. It is reaily seen that the Big matrix satisfies the Hawkins-Simon conitions for an acceptable input-output matrix if the matrices i o. This Big matrix has become familiar because the multiregional input-output moel use by regional economists generally takes the following form (see, for example, (eters an Hertwich forthcoming: (I x = y + e, or I I I x x x = y y y + e e e, (6- which can be solve for x, given y (here efine as final eman for only those goos prouce an consume in a region an the vectors of exports to final eman in other regions, e i. Using the matrices F, ij, an, Φ can be calculate from Eq. (-. Thus we have expressions for Φ: Eq. (6- shows the comparable Big an path-base expressions, an Eqs. (6-4 an (6-5 evelop them to highlight the similarities an ifferences. In Eq. (6- an (6-5, is a block-iagonal matrix, an each block escribes factor flows ( ij = f ij x j in a region. Φ F(I ~ y (I y (6-

29 9 y y y I I I F F F Φ (6-4 y y y I e e e I e e e I Φ (6-5 Formulating the problem in terms of the Big matrix has the avantage that manipulating the Big matrix is alreay familiar to regional input-output economists: it is hanle like a one-region input-output matrix, an its Leontief inverse is compute. The properties of the Leontief inverse are sufficiently attractive that analysts are incline to put all extensions to the basic static input-output moel in the form of Eq. (-a. Eq. (6- is a multi-regional moel written in that form, an the ynamic input-output moel an the moel close for househols have also been written an solve in a Big form, among other formulations (but naturally with a ifferent Big matrix, for this reason. From a practical point of view, statistical offices have begun to publish the breakown of a country s input-output table into a omestic input-output table an an import table -- although the latter are not now further broken own to istinguish iniviual country origins of imports in response to the eman of analysts. However, the solution by the Big metho of Eq. (6-4 has several rawbacks. First, the numerous off-iagonal components of the Big matrix are not available in publishe ata. Even when a regionally-aggregate import matrix is available, it oes not contribute aitional information content over the smaller quantity of ata require for the Markov metho, as import matrices are generally erive on the assumption that the omestic an import shares for a given goo are the same in all uses of that goo. That is, ij = i r ji ŝ i, where i r is the iagonal matrix showing the share of imports in omestic supply in region i an ji ŝ shows the share of all imports to region i from region j. There are many categories of ata that analysts woul like from statistical offices, but arguably better ij matrices for benchmark years base on the collection of more etaile ata an final eman vectors istinguishing omestically prouce from importe goos - shoul not be accore high priority. While matrices are a goo choice of parameters because they change slowly an in ways that can be unerstoo in terms of technological change, this is not the case for ij matrices where a proucer s choice between a omestically prouce input an its similarly-price importe counterpart, or between an import of a particular goo from one region rather than another region at similar prices is circumstantial an volatile. By contrast, the off-iagonal blocks in the path-base approach are iagonal matrices with exports as a share of omestic supply own the iagonal.

30 secon rawback of the Big approach is conceptual. It is natural, an generally effective, to approach new research problems by making use of concepts an techniques that worke well for relate problems in the past; this is the reason the Leontief inverse an its constituent multipliers are put in the service of new problems. However, this practice also has the effect of slowing the aoption of new approaches able to aress a broaer set of questions by generaliing existing concepts an methos. Duchin an Steenge (7 provie a critique of the Big strategy. The path-base approach set out to follow the physical flows rather than to erive an extene Leontief inverse. The analysis in Sections 4 an 5 provies a metho to quantify the factors (an emissions emboie in a region s consumption, taking full account of its imports, both irect an inirect, an their transport. Resources use in transporting imports are reallocate to the importe goos. The T matrices require to internalie these inputs o not fit into a Big matrix because they are multiplie not by an output vector but by a region s imports from each istinct origin. Distinguishing the transport associate with trae flows requires a more general representation. Finally, the absorbing Markov chain analysis offers aitional information about the input-output problem that is not available using stanar input-output methos, that has not been exploite in the present paper. In particular, the N matrix, erive from the canonical Markov matrix, measures the number of times a noe is visite. This property is use by Yamaa et al. (6 an Matsuno et al. (7 to examine the paths of recycle a material through an economy an count the average number of times it is recycle before ening up in a lanfill. Similarly, it was utilie by Levine (98 to erive a measure of trophic position within an ecosystem. The approach can profitably be generalie to a global analysis such as that reporte here. The Big metho has the avantage of familiarity among input-output analysts, by contrast with the absorbing Markov chain analysis propose in the paper. Our claim is that the approaches are complementary, enjoying areas of overlap where they provie results that are ientical, an areas of non-overlap where the Markov results require minimal information to answer familiar questions an provie consistent answers to a larger set of emerging questions. The basic logic of the Markov approach is natural an intuitive, as illustrate by the iagrams of physical flows shown in Figures., 4., an 4.. The efining characteristic of the metho is to associate each flow of a factor or goo with a set of shares aing to. that inicate the estination (or origin in the case of pollution of that flow. Once this simple logic is internalie, the mathematics of the Markov metho is as straightforwar as that for manipulating an input-output moel. However, both the Big metho an the Markov chain analysis escribe in this paper are ex-post approaches to analying historical ata or ata generate by a worl moel; neither is a substitute for a moel capable of analying alternative scenarios.

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

Chapter 3 Notes, Applied Calculus, Tan

Chapter 3 Notes, Applied Calculus, Tan Contents 3.1 Basic Rules of Differentiation.............................. 2 3.2 The Prouct an Quotient Rules............................ 6 3.3 The Chain Rule...................................... 9 3.4

More information

A Review of Multiple Try MCMC algorithms for Signal Processing

A Review of Multiple Try MCMC algorithms for Signal Processing A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization Web-Base Technical Appeni: Multi-Prouct Firms an Trae Liberalization Anrew B. Bernar Tuck School of Business at Dartmouth & NBER Stephen J. Reing LSE, Yale School of Management & CEPR Peter K. Schott Yale

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize

More information

Generalizing Kronecker Graphs in order to Model Searchable Networks

Generalizing Kronecker Graphs in order to Model Searchable Networks Generalizing Kronecker Graphs in orer to Moel Searchable Networks Elizabeth Boine, Babak Hassibi, Aam Wierman California Institute of Technology Pasaena, CA 925 Email: {eaboine, hassibi, aamw}@caltecheu

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

arxiv:hep-th/ v1 3 Feb 1993

arxiv:hep-th/ v1 3 Feb 1993 NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

lim Prime notation can either be directly applied to a function as previously seen with f x 4.1 Basic Techniques for Finding Derivatives

lim Prime notation can either be directly applied to a function as previously seen with f x 4.1 Basic Techniques for Finding Derivatives MATH 040 Notes: Unit Page 4. Basic Techniques for Fining Derivatives In the previous unit we introuce the mathematical concept of the erivative: f f ( h) f ( ) lim h0 h (assuming the limit eists) In this

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Dot trajectories in the superposition of random screens: analysis and synthesis

Dot trajectories in the superposition of random screens: analysis and synthesis 1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Fluid Models with Jumps

Fluid Models with Jumps Flui Moels with Jumps Elena I. Tzenova, Ivo J.B.F. Aan, Viyahar G. Kulkarni August 2, 2004 Abstract In this paper we stuy a general stochastic flui moel with a single infinite capacity buffer, where the

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10 Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Quality competition versus price competition goods: An empirical classification

Quality competition versus price competition goods: An empirical classification HEID Working Paper No7/2008 Quality competition versus price competition goos: An empirical classification Richar E. Balwin an Taashi Ito Grauate Institute of International an Development Stuies Abstract

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

Monopoly Part III: Multiple Local Equilibria and Adaptive Search

Monopoly Part III: Multiple Local Equilibria and Adaptive Search FH-Kiel University of Applie Sciences Prof Dr Anreas Thiemer, 00 e-mail: anreasthiemer@fh-kiele Monopoly Part III: Multiple Local Equilibria an Aaptive Search Summary: The information about market eman

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

2Algebraic ONLINE PAGE PROOFS. foundations

2Algebraic ONLINE PAGE PROOFS. foundations Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using

More information

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3 RAIRO Operations Research RAIRO Oper. Res. 41 (27) 455 464 DOI: 1.151/ro:2731 NEWSVENDOR SOLUTIONS WITH GENERAL RANDOM YIELD DISTRIBUTIONS Scott E. Grasman 1, Zaki Sari 2 an Tewfik Sari 3 Abstract. Most

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

NBER WORKING PAPER SERIES FACTOR PRICES AND INTERNATIONAL TRADE: A UNIFYING PERSPECTIVE. Ariel Burstein Jonathan Vogel

NBER WORKING PAPER SERIES FACTOR PRICES AND INTERNATIONAL TRADE: A UNIFYING PERSPECTIVE. Ariel Burstein Jonathan Vogel NBER WORKING PAPER SERIES FACTOR PRICES AND INTERNATIONAL TRADE: A UNIFYING PERSPECTIVE Ariel Burstein Jonathan Vogel Working Paper 6904 http://www.nber.org/papers/w6904 NATIONAL BUREAU OF ECONOMIC RESEARCH

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection Online Appenix for Trae Policy uner Monopolistic Competition with Firm Selection Kyle Bagwell Stanfor University an NBER Seung Hoon Lee Georgia Institute of Technology September 6, 2018 In this Online

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Web Appendix to Firm Heterogeneity and Aggregate Welfare (Not for Publication)

Web Appendix to Firm Heterogeneity and Aggregate Welfare (Not for Publication) Web ppeni to Firm Heterogeneity an ggregate Welfare Not for Publication Marc J. Melitz Harvar University, NBER, an CEPR Stephen J. Reing Princeton University, NBER, an CEPR March 6, 203 Introuction his

More information

Consider for simplicity a 3rd-order IIR filter with a transfer function. where

Consider for simplicity a 3rd-order IIR filter with a transfer function. where Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation

More information

A Path Planning Method Using Cubic Spiral with Curvature Constraint

A Path Planning Method Using Cubic Spiral with Curvature Constraint A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw

More information

Two formulas for the Euler ϕ-function

Two formulas for the Euler ϕ-function Two formulas for the Euler ϕ-function Robert Frieman A multiplication formula for ϕ(n) The first formula we want to prove is the following: Theorem 1. If n 1 an n 2 are relatively prime positive integers,

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences. S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial

More information

Basic IIR Digital Filter Structures

Basic IIR Digital Filter Structures Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

What s in an Attribute? Consequences for the Least Common Subsumer

What s in an Attribute? Consequences for the Least Common Subsumer What s in an Attribute? Consequences for the Least Common Subsumer Ralf Küsters LuFG Theoretical Computer Science RWTH Aachen Ahornstraße 55 52074 Aachen Germany kuesters@informatik.rwth-aachen.e Alex

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Text S1: Simulation models and detailed method for early warning signal calculation

Text S1: Simulation models and detailed method for early warning signal calculation 1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009 Spurious Significance of reatment Effects in Overfitte Fixe Effect Moels Albrecht Ritschl LSE an CEPR March 2009 Introuction Evaluating subsample means across groups an time perios is common in panel stuies

More information

4.2 First Differentiation Rules; Leibniz Notation

4.2 First Differentiation Rules; Leibniz Notation .. FIRST DIFFERENTIATION RULES; LEIBNIZ NOTATION 307. First Differentiation Rules; Leibniz Notation In this section we erive rules which let us quickly compute the erivative function f (x) for any polynomial

More information

arxiv: v1 [cs.lg] 22 Mar 2014

arxiv: v1 [cs.lg] 22 Mar 2014 CUR lgorithm with Incomplete Matrix Observation Rong Jin an Shenghuo Zhu Dept. of Computer Science an Engineering, Michigan State University, rongjin@msu.eu NEC Laboratories merica, Inc., zsh@nec-labs.com

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE ELECTRIC CHARGE Introuction: Orinary matter consists of atoms. Each atom consists of a nucleus, consisting of protons an neutrons, surroune by a number of electrons. In electricity, the electric charge

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

New Technology, Human Capital and Growth for Developing Countries.

New Technology, Human Capital and Growth for Developing Countries. New Technology, Human Capital an Growth for Developing Countries Cuong Le Van, Manh-Hung Nguyen an Thai Bao Luong Centre Economie e la Sorbonne, Université Paris-1, Pantheon-Sorbonne 106-112 B e l Hôpital,

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca

More information

u!i = a T u = 0. Then S satisfies

u!i = a T u = 0. Then S satisfies Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Computing Derivatives

Computing Derivatives Chapter 2 Computing Derivatives 2.1 Elementary erivative rules Motivating Questions In this section, we strive to unerstan the ieas generate by the following important questions: What are alternate notations

More information

LeChatelier Dynamics

LeChatelier Dynamics LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a

More information

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification Why Bernstein Polynomials Are Better: Fuzzy-Inspire Justification Jaime Nava 1, Olga Kosheleva 2, an Vlaik Kreinovich 3 1,3 Department of Computer Science 2 Department of Teacher Eucation University of

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

More information

Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Method

Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Method Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Metho Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Computing Inverse Trig Derivatives. Starting with the inverse

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

A tale of two cities:

A tale of two cities: Highly preliminary. Do not quote or refer to. A tale of two cities: Urban spatial structure an moe of transportation Takaaki Takahashi Center for Spatial Information Science, the University of Tokyo, 5-1-5,

More information