Matlab Calculate Roots of A Quadratic Equation

Calculating the roots of a quadratic equation is a fundamental mathematical operation that finds applications in various fields including physics, engineering, and economics. This guide explains how to perform this calculation in MATLAB, a powerful numerical computing environment.

Introduction

A quadratic equation is a second-degree polynomial equation in a single variable, typically written in the form:

ax² + bx + c = 0

Where a, b, and c are constants, and a ≠ 0. The roots of the equation are the values of x that satisfy the equation. There are three possible cases for the roots of a quadratic equation:

Two distinct real roots
One real root (a repeated root)
No real roots (complex roots)

MATLAB provides several methods to calculate the roots of a quadratic equation, including the roots() function and direct calculation using the quadratic formula.

Quadratic Equation Formula

The roots of a quadratic equation can be found using the quadratic formula:

x = [-b ± √(b² - 4ac)] / (2a)

Where:

a, b, c are the coefficients of the quadratic equation
√(b² - 4ac) is the discriminant
The discriminant determines the nature of the roots:

If the discriminant is positive, there are two distinct real roots. If it's zero, there's one real root. If it's negative, there are two complex conjugate roots.

MATLAB Method

MATLAB provides several ways to calculate the roots of a quadratic equation. The most straightforward method is to use the roots() function, which returns the roots of a polynomial given its coefficients.

roots([a b c])

For example, to find the roots of x² - 5x + 6 = 0, you would use:

roots([1 -5 6])

Alternatively, you can implement the quadratic formula directly in MATLAB:

a = input('Enter coefficient a: '); b = input('Enter coefficient b: '); c = input('Enter coefficient c: '); discriminant = b^2 - 4*a*c; if discriminant > 0 root1 = (-b + sqrt(discriminant))/(2*a); root2 = (-b - sqrt(discriminant))/(2*a); fprintf('Two distinct real roots: %.4f and %.4f\n', root1, root2); elseif discriminant == 0 root = -b/(2*a); fprintf('One real root: %.4f\n', root); else realPart = -b/(2*a); imagPart = sqrt(abs(discriminant))/(2*a); fprintf('Two complex roots: %.4f + %.4fi and %.4f - %.4fi\n', realPart, imagPart, realPart, imagPart); end

Worked Example

Let's calculate the roots of the quadratic equation x² - 4x + 4 = 0 using MATLAB.

First, identify the coefficients: a = 1, b = -4, c = 4
Calculate the discriminant: (-4)² - 4*1*4 = 16 - 16 = 0
Since the discriminant is zero, there is one real root
Calculate the root: x = -(-4)/(2*1) = 2

In MATLAB, this would be calculated as:

roots([1 -4 4])

The result would be [2], confirming our manual calculation.

Frequently Asked Questions

What is the difference between the roots() function and the quadratic formula in MATLAB?

The roots() function is a general-purpose function that can find roots of any polynomial, while the quadratic formula is specifically designed for quadratic equations. Both methods will give the same results for quadratic equations.

How do I handle complex roots in MATLAB?

MATLAB automatically handles complex roots when the discriminant is negative. The roots will be returned in complex number format, and you can display them using fprintf with appropriate formatting.

Can I use MATLAB to plot the quadratic equation and its roots?

Yes, you can use MATLAB's plotting functions to visualize the quadratic equation and its roots. The fplot() function can be used to plot the equation, and the roots can be marked on the plot using the plot() function.