New User Special Price Expires in

Let's log you in.

Sign in with Facebook


Don't have a StudySoup account? Create one here!


Create a StudySoup account

Be part of our community, it's free to join!

Sign up with Facebook


Create your account
By creating an account you agree to StudySoup's terms and conditions and privacy policy

Already have a StudySoup account? Login here


by: Vivien Bradtke V


Vivien Bradtke V
Texas A&M
GPA 3.6

Joseph Pasciak

Almost Ready


These notes were just uploaded, and will be ready to view shortly.

Purchase these notes here, or revisit this page.

Either way, we'll remind you when they're ready :)

Preview These Notes for FREE

Get a free preview of these Notes, just enter your email below.

Unlock Preview
Unlock Preview

Preview these materials now for free

Why put in your email? Get access to more of this material and other relevant free materials for your school

View Preview

About this Document

Joseph Pasciak
Class Notes
25 ?




Popular in Course

Popular in Mathematics (M)

This 2 page Class Notes was uploaded by Vivien Bradtke V on Wednesday October 21, 2015. The Class Notes belongs to MATH 609 at Texas A&M University taught by Joseph Pasciak in Fall. Since its upload, it has received 14 views. For similar materials see /class/226031/math-609-texas-a-m-university in Mathematics (M) at Texas A&M University.

Similar to MATH 609 at Texas A&M

Popular in Mathematics (M)




Report this Material


What is Karma?


Karma is the currency of StudySoup.

You can buy or earn more Karma at anytime and redeem it for class notes, study guides, flashcards, and more!

Date Created: 10/21/15
Programming assignment3 Due Nov 21 Math 609 We consider applying explicit and implicit methods for solving a mildly stiff system of ODE s in this assignment Speci cally7 we consider the system Y Yt 01 a R where Y is the solution to n 12 Y iAY 0 01 t 5 Y0 Y0 E 111 Here A is the n x 71 matrix given by 2 ifz j A 71 iflz ijl17 0 otherwise Problem 1 Show that the set of vectors 117 7 1 C R with compo nents given by i ij gtgtjsmn1 241m are eigenvectors for A and compute the corresponding eigenvalues Write down the solution to 01 given the expansion of the initial data7 V L Y0 Zimm i1 You need not compute the coef cients above and they will appear in your expansion of the solution Note that every component of the solution decays exponentially as 25 increases Problem 2 Given a time step size k let Y denote the sequence obtained by applying the forward Euler method to 01 Give an expression for the coef cients of Y expanded in the basis of eigenvectors The forward Euler method is considered stable for this problem provided that every component of the solution decays Find a resonably sharp see the next problem77 condition on the time step size as a function of n which guarantees stability Problem 3 Code the forward Euler method applied to the ODE 01 You can use CRS to store the matrix or simply store the diagonals in three vectors The forward Euler method only requires matrix evaluation and is simple to code using either storage mode note that the generation of a full 71 x 71 matrix is unacceptable Apply your method to the cases of 1 2 n 16 32 64 128 First use time step sizes which are twice as large as your condition in the previous problem If your condition above is sharp enough you will see blowup of the ODE approximation Next use time step sizes which are smaller than you condition by a factor of two If your condition above is sharp enough you will now see convergence at t 1 Problem 4 Code the backward Euler method applied to the ODE 01 As the matrix A is tri diagonal the matrix systems appearing are also tri diagonal These can be solved using a tri diagonal solver as illustrated in the c code httpwwwmathtamuedu pasciakclasses609trisolc This code is simply Gaussian elimination in a storage mode where calculations involving the zeroes of the matrix are avoided As above the construction of the full 71 x 71 matrix is not allowed This method is stable for any time step size Run the method for n 163264128 and in each case try to nd a time step size where you observe two decimal places of accuracy in each component This you can do by comparing the results of different decreasing time step sizes while observing the resulting changes in the solution at t 1 You should be able to get the desired accuracy with fewer time steps than necessary for stability of the forward Euler method above


Buy Material

Are you sure you want to buy this material for

25 Karma

Buy Material

BOOM! Enjoy Your Free Notes!

We've added these Notes to your profile, click here to view them now.


You're already Subscribed!

Looks like you've already subscribed to StudySoup, you won't need to purchase another subscription to get this material. To access this material simply click 'View Full Document'

Why people love StudySoup

Steve Martinelli UC Los Angeles

"There's no way I would have passed my Organic Chemistry class this semester without the notes and study guides I got from StudySoup."

Kyle Maynard Purdue

"When you're taking detailed notes and trying to help everyone else out in the class, it really helps you learn and understand the I made $280 on my first study guide!"

Jim McGreen Ohio University

"Knowing I can count on the Elite Notetaker in my class allows me to focus on what the professor is saying instead of just scribbling notes the whole time and falling behind."

Parker Thompson 500 Startups

"It's a great way for students to improve their educational experience and it seemed like a product that everybody wants, so all the people participating are winning."

Become an Elite Notetaker and start selling your notes online!

Refund Policy


All subscriptions to StudySoup are paid in full at the time of subscribing. To change your credit card information or to cancel your subscription, go to "Edit Settings". All credit card information will be available there. If you should decide to cancel your subscription, it will continue to be valid until the next payment period, as all payments for the current period were made in advance. For special circumstances, please email


StudySoup has more than 1 million course-specific study resources to help students study smarter. If you’re having trouble finding what you’re looking for, our customer support team can help you find what you need! Feel free to contact them here:

Recurring Subscriptions: If you have canceled your recurring subscription on the day of renewal and have not downloaded any documents, you may request a refund by submitting an email to

Satisfaction Guarantee: If you’re not satisfied with your subscription, you can contact us for further help. Contact must be made within 3 business days of your subscription purchase and your refund request will be subject for review.

Please Note: Refunds can never be provided more than 30 days after the initial purchase date regardless of your activity on the site.