Initial value problems, math, Mathematics

Assignment Help:
Write a Matlab function MyIVP that solves an initial-value problem (IVP) for a system of ordinary
differential equations (ODEs) of the form
x ?(t) = f (t, x(t)),
where f : R × Rn ? Rn is an arbitrary function with one one-dimensional input (for time t) and one n-dimensional input x, and n-dimensional output. The function should implement a Runge-Kutta formula (for example, the rk34 formula or the Dormand & Prince formula).
The initial value x0 is provided by the user of MyIVP. The first line of MyIVP (saved in a file MyIVP.m) should look like this
function [xend,t,xt]=MyIVP(f,x0,tspan,N) Inputs
• f: function defining the right-hand side of the ODE. f should accept two arguments: t (a number) and x (an n-dimensional vector). The function f should return an n-dimensional vector y (the time derivative). Typical calling sequence: y=f(t,x), returning the value of f at time t in position x.
• x0: initial value where integration starts from (n-dimensional vector).
• tspan: Starting time and end time for integration. Integration has to run from time t =tspan(1)
to time t =tspan(2).
• N: number of steps for integration. The integration stepsize h=(tspan(2)-tspan(1))/N should
be small.
Outputs
• xend: result of integration at t =tspan(2).
• t: vector of times at which intermediate values have been computed (this should have N + 1
entries).
• xt: intermediate values (n × (N + 1)-array). xt(:,k) should be the solution at t(k).
You can check the built-in variable nargout inside your function to see if the user wants to get three outputs or only the end value xend. If nargout==1 you don’t need to store the intermediate values.

Related Discussions:- Initial value problems, math

Opening Account, I am expert in mathematics. How i open my expert account?

I am expert in mathematics. How i open my expert account?

Solution by factorization, Solution by Factorization, please solve quadrati...

Solution by Factorization, please solve quadratic equations by Factorization.

Right- and left-handed limits , Right- and left-handed limits : Next, let'...

Right- and left-handed limits : Next, let's see precise definitions for the right- & left-handed limits. Definition   For the right-hand limit we say that, if for eve

Show that tan = 1/v3 , If 7sin 2 ?+3cos 2 ? = 4, show that tan? =   1/√3  ...

If 7sin 2 ?+3cos 2 ? = 4, show that tan? =   1/√3                      . Ans:    If 7 Sin 2 ? + 3 Cos 2 ? = 4 S.T. Tan?  1/√3 7 Sin 2 ? + 3 Cos 2 ? = 4 (Sin 2 ? + Cos 2 ?)

Determine the widest piece of mail, A mailbox opening is 4.5 inches high an...

A mailbox opening is 4.5 inches high and 5 inches wide. Determine the widest piece of mail able to ?t in the mailbox without bending? a. 9.5 inches b. 2.2 inches c. 6.7 in

Calculate expected average time , A car was machine washes every car in 5 m...

A car was machine washes every car in 5 minutes accurately. It has been calculated that customers will arrive as per to a Poisson distribution at an average of 8 per hour. Calculat

Prove that the poset has a unique least element, Prove that the Poset has a...

Prove that the Poset has a unique least element Prove that if (A, ) has a least element, then (A,≤)  has a unique least element. Ans: Let (A, ≤) be a poset. Suppose the po

Solve the initial value by laplace transform method, Question: Solve the i...

Question: Solve the initial value problem 2x'' +x'-x =27 Cos2t +6 Sin 2t, x(0)=2 , x'(0)= -2 by using Laplace transform method.

Write Your Message!

Captcha
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd