If only one fertilizer x is used, the dependence of yield z(x) on x first was given by Mitscherlich (1909) in form of the differential equation

Similar documents
Lecture 1 - Introduction and Basic Facts about PDEs

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx,

ANALYSIS AND MODELLING OF RAINFALL EVENTS

MAC-solutions of the nonexistent solutions of mathematical physics

University of Sioux Falls. MAT204/205 Calculus I/II

Solutions to Assignment 1

ON THE INEQUALITY OF THE DIFFERENCE OF TWO INTEGRAL MEANS AND APPLICATIONS FOR PDFs

Section 4.4. Green s Theorem

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

Gauss Quadrature Rule of Integration

Integrals Depending on a Parameter

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

K 7. Quadratic Equations. 1. Rewrite these polynomials in the form ax 2 + bx + c = 0. Identify the values of a, b and c:

AP Calculus AB Unit 4 Assessment

( ) { } [ ] { } [ ) { } ( ] { }

f (z) dz = 0 f(z) dz = 2πj f(z 0 ) Generalized Cauchy Integral Formula (For pole with any order) (n 1)! f (n 1) (z 0 ) f (n) (z 0 ) M n!

Core 2 Logarithms and exponentials. Section 1: Introduction to logarithms

1B40 Practical Skills

ILLUSTRATING THE EXTENSION OF A SPECIAL PROPERTY OF CUBIC POLYNOMIALS TO NTH DEGREE POLYNOMIALS

MATH Final Review

2008 Mathematical Methods (CAS) GA 3: Examination 2

Population Dynamics Definition Model A model is defined as a physical representation of any natural phenomena Example: 1. A miniature building model.

Linear Algebra Introduction

Co-ordinated s-convex Function in the First Sense with Some Hadamard-Type Inequalities

A STUDY OF Q-CONTIGUOUS FUNCTION RELATIONS

Lecture 2e Orthogonal Complement (pages )

#A42 INTEGERS 11 (2011) ON THE CONDITIONED BINOMIAL COEFFICIENTS

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016

The study of dual integral equations with generalized Legendre functions

1. Extend QR downwards to meet the x-axis at U(6, 0). y

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Lecture 21: Order statistics

CHM Physical Chemistry I Chapter 1 - Supplementary Material

Jim Lambers MAT 280 Spring Semester Lecture 26 and 27 Notes

Precalculus Spring 2017

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Chapter Gauss Quadrature Rule of Integration

Part 4. Integration (with Proofs)

Gauss Quadrature Rule of Integration

A Non-parametric Approach in Testing Higher Order Interactions

I. INTEGRAL THEOREMS. A. Introduction

10 Vector Integral Calculus

We divide the interval [a, b] into subintervals of equal length x = b a n

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Discrete Least-squares Approximations

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Best Approximation. Chapter The General Case

MATRICES AND VECTORS SPACE

1 Linear Least Squares

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

than 1. It means in particular that the function is decreasing and approaching the x-

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

Chapter 9 Definite Integrals

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

Kinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums

Polynomials. Polynomials. Curriculum Ready ACMNA:

Review of Gaussian Quadrature method

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Handout 4. Inverse and Implicit Function Theorems.

Math 259 Winter Solutions to Homework #9

Predict Global Earth Temperature using Linier Regression

Learning Objectives of Module 2 (Algebra and Calculus) Notes:

Hyers-Ulam stability of Pielou logistic difference equation

On the Scale factor of the Universe and Redshift.

AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION

Polynomials and Division Theory

Journal of Inequalities in Pure and Applied Mathematics

Introduction to Olympiad Inequalities

On Error Sum Functions Formed by Convergents of Real Numbers

Absolute values of real numbers. Rational Numbers vs Real Numbers. 1. Definition. Absolute value α of a real

Orthogonal Polynomials and Least-Squares Approximations to Functions

ON THE EXCEPTIONAL SET IN THE PROBLEM OF DIOPHANTUS AND DAVENPORT

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

Lecture Solution of a System of Linear Equation

Math 100 Review Sheet

Chapter 8 Roots and Radicals

Propositional models. Historical models of computation. Application: binary addition. Boolean functions. Implementation using switches.

Linear Inequalities. Work Sheet 1

The Periodically Forced Harmonic Oscillator

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Section 6: Area, Volume, and Average Value

AT100 - Introductory Algebra. Section 2.7: Inequalities. x a. x a. x < a

TOPIC: LINEAR ALGEBRA MATRICES

A Study on the Properties of Rational Triangles

A Matrix Algebra Primer

1 ELEMENTARY ALGEBRA and GEOMETRY READINESS DIAGNOSTIC TEST PRACTICE

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

For a, b, c, d positive if a b and. ac bd. Reciprocal relations for a and b positive. If a > b then a ab > b. then

Research Article Harmonic Deformation of Planar Curves

ad = cb (1) cf = ed (2) adf = cbf (3) cf b = edb (4)

Homework Assignment 6 Solution Set

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

The mth Ratio Convergence Test and Other Unconventional Convergence Tests

Necessary and sucient conditions for some two. Abstract. Further we show that the necessary conditions for the existence of an OD(44 s 1 s 2 )

Transcription:

Mitsherlih s Lw: Generliztion with severl Fertilizers Hns Shneeerger Institute of Sttistis, University of Erlngen-Nürnerg, Germny 00, 5 th August Astrt: It is shown, tht the rop-yield z in dependene on two fertilizers x nd y is the produt of two omponents: z in dependene on x lone nd z in dependene on y lone, divided y, the yield without extern fertilizers, i.e. with xy0. For n fertilizers, we hve the n produt of n omponents, divided y. Introdution If only one fertilizer x is used, the dependene of yield z(x) on x first ws given y Mitsherlih (909) in form of the differentil eqution dẑ dx ( ẑ) () where is the symptoti vlue of z, the ftor of proportionlity. As usul in sttistis, ẑ is the hypothetil vlue, z the experimentl vlue of the rop-yield. For eqution () it is ssumed: No over-fertiliztion. For the se of overfertiliztion see Shneeerger (009). Solution of formul () with oundry ondition ẑ (x 0) is Mitsherlih s urve in the espeil instrutive form x ẑ(x) + ( )( e ) () demonstrted in figure with 0.5,., 0,75. For estimtion of the prmeters in prtie see Shneeerger (009). Figure : Crop-yield ẑ(x) s funtion of one fertilizer x

An other form of formul () is: used in the following. ẑ(x) x ( )e () Bule (98) gve the solution of eqution () in the form (x d) ẑ ( e ) (3) with d ln (<0) (3) in figure we hve d-0.79. Generliztion Now we ssume two vriles (fertilizers) x nd y (see figure ) Figure : Crop-yield ẑ (x, y) s funtion of two fertilizers x nd y

Then we hve for the rop-yield ẑ (x, y) with formul () ẑ(0, y) nd ẑ(x,0) y ( )e (4) x (x) ( )e (4) Note: For short we write (y 0), (x 0), (y 0) (x 0), (y 0), (x 0), d (y 0) d, d (x 0) d. Generlizing formule (4) we hve ẑ(x, y) ẑ(x, y) (x)y (x) ( (x) (x))e (5) x ( )e (5) nd herewith (x) x (6) (x) (x) (x)y e (x) e Now I mke use of result of Mitsherlih (947): given different fertilizers, the prmeter (in Mitsherlih s nottion) is onstnt for fixed fertilizer, s I ould prove in tens of yers of work. Mitsherlih s Wirkungsgrd is our prmeter exept for the onstnt ln0. This mens: Herewith we get from formul (6): nd (x) independent of y resp. x (7) is independent of y, i.e. (y 0) (8) (y 0) nd (x) independent of x (x) (8) nd with this d independent of y (9) ln ln d (x) d nd d independent of x. (9) Finlly we hve from formul (5) (or (5)): y y ẑ(x, y) (x) (x)e (x)( ( )e ) (x) (0) or in the most instrutive form x y ẑ(x, y) ( + ( )( e ))( + ( )( e )) ()

RESULT: The generlized Mitsherlih formul in two vriles is the produt of the onedimensionl formule, multiplied y /. With formule (9) one n show, tht formul () is identil with the formul of Bule (98). ẑ(x, y) (x d ) (y d ) A( e )( e ) with () A Equtions (0) nd () n esily e generlized for n fertilizers: Applition ẑ(x,...x ) ẑ(x,0,...0)ẑ(0, x,...0)...ẑ(0,0,...x n ) (3) n n The following dt re from n exmple of Steinhuser, Lngehn nd Peters (99) with x (in 00 kg/h of P O 5 ), y (in 00 kg/h of O), z (in 000kg/h of rye). K Tle : Crop-yield ẑ (x, y) in 000 kg/h of rye, x in 00 kg/h of P O 5, y in 00 kg/h of K O x 0.5 x 0. 50 x 3 0. 75 x 4. 00 x 5. 5 x 6. 50 y 0.5.00..4.58.73.87 y 0.50.4.79.09.34.55.73 y 3 0.75.7..59.90 3.5 3.35 y 4.00.00.55.98 3.3 3.59 3.8 y 5.5..8 3.9 3.65 3.94 4.8 y 6.50.4 3.05 3.55 3.93 4.4 4.50 The 5 prmeters,,,, were determined with the method of Lest Squres of Guss f (,,,, ) (z(x, y) ẑ(x, y)) Min (4) x y summing over ll 36 dt-points. The minimum ws gined itertively with the non-liner Simplex-Method of Nelder nd Med (965). The result is 0.75, 0.947, 0.899,.9438,. 07 nd ẑ 3.683(0.75 + 0.6757( e 0.899x ))(0.75 +.07y.673(- e )) In figure 3 the ontour-lines ẑ (x, y) 0., 0.4, 4.0 re drwn, in figure 3 interseting urves ẑ (x, y onst.) for y d,0, 0.5,.0,.5 nd re plotted. It is ovious tht they re Mitsherlih-urves. We hve d 0.376, d --0.46. The symptotes of the urves ẑ (x, y ) re horizontl dotted stright lines in figure 3. Espeilly for y we get

ẑ(, ) 6.78, Bule s prmeter A in formul (). ẑ(x, y ) is urve (x) of figure. In nlogy interseting urves for xonst. ould e plotted. Figure 3: Contour-lines ẑ (x, y) onst. - in 000 kg/h of rye, x in 00 kg/h of P O 5, y in 00 kg/h of K O

Figure 3: Mitsherlih-urves ẑ (x, y onst.) Exmple: With fertilizer 00 kg/h of P O 5 (x) nd 50 kg/h of K O (y0.5) we get the rop-yield 335 kg/h of rye ( ẑ.335).. With xy0 (i.e. without externl fertilizers) we would get 7 kg/h of rye ( ẑ 0.7). Aknowledgement I hve to thnk Dr. Emher, Munih, who helped me to pulish these ppers in the internet. A generliztion with overfertiliztion is given in pper 5 (Pper 5: Mitsherlih's Lw: Generliztion with severl Fertilizers nd Overfertiliztion) Referenes Bule B. (98). Zu Mitsherlihs Gesetz der physiologishen Beziehungen, Lndwirtshftlihe Jhrüher 5, 363-385 Mitsherlih E.A. (909). Ds Gesetz des Minimums und ds Gesetz des nehmenden Bodenertrgs, Lndwirtshftlihe Jhrüher 38, 537-55 Mitsherlih E.A. (947). Ds Ergenis von üer 7000 Feld-Düngungsversuhen, Zeitshrift für Pflnzenernährung, Düngung, Bodenkunde 38. Bnd, 947, Verlg Chemie Nelder J.R. nd Med R. (965). A Simplex Method for funtion minimiztion. The Computer Journl 7, 303-33 Shneeerger H. (009). Mitsherlih s Lw: Sum of two Exponentil Proesses. Conlusions. Internet: www.soil-sttisti.de, pper nd pper Shneeerger H. (009). Over-Fertiliztion: An Inverse Mitsherlih Proess. Internet: www.soil-sttisti.de, pper 3 Steinhuser H., Lngehn C. nd Peters U. (99). Einführung in die lndwirtshftlihe Betrieslehre, Bnd, 5.Auflge