// Composite Trapezoidal Rule ---- Assignemnt 1 --- Question A #define _USE_MATH_DEFINES #include<stdio.h> #include <cmath> #include <iostream> #include <fstream> #include <cmath> #include <tuple> using namespace std; // Function --- f(x) = 4/(1+x^{2}) float y(float x) { // Declaring the computational function f(x) = 4/(1+x^{2}) return 4/(1+x*x); } // Computational function to evaluate the definite Integral float trapezoidal(float lowlim, float uplim, float n) { // Grid Spacing or Width of the Trapezoid along x-Axis float h = (uplim-lowlim)/n; // Calculating the f(lowlim) and f(uplim) summation float s = y(lowlim)+y(uplim); // Calculating the Middle Terms of the Composite Trapezoid Formula for (int i = 1; i < n; i++) s += 2*y(lowlim+i*h); // Multiplying with Grid Spacing by 2 (Sum of Parallel Sides * (Height/2)) return (h/2)*s; } int main() { // Range of Definite Integral float x0 = 0; float xn = 1; printf("PI %6.15f\n", M_PI ); printf("\n"); long long int n = 5; //Loop to Evaluate the Answer for Multiple Values of n while( n < 1000000 ){ // Printing the Answer to a 20 Digit Precision printf("Value of integral is = %6.15f\n", trapezoidal(x0, xn, n)); double e = abs(M_PI - trapezoidal(x0, xn, n)); double loge = log(abs(M_PI - trapezoidal(x0, xn, n))); double logn = log(n); printf("log(n) = %6.15f\n", logn); printf("Error in the Intgral Computation is = %6.18f\n", e); printf("log(abs(error)) = %6.15f\n", loge); printf("\n"); tuple <double, double> var; var = make_tuple(logn, loge); n = n*5; } return 0; }
Write, Run & Share C++ code online using OneCompiler's C++ online compiler for free. It's one of the robust, feature-rich online compilers for C++ language, running on the latest version 17. Getting started with the OneCompiler's C++ compiler is simple and pretty fast. The editor shows sample boilerplate code when you choose language as C++
and start coding!
OneCompiler's C++ online compiler supports stdin and users can give inputs to programs using the STDIN textbox under the I/O tab. Following is a sample program which takes name as input and print your name with hello.
#include <iostream>
#include <string>
using namespace std;
int main()
{
string name;
cout << "Enter name:";
getline (cin, name);
cout << "Hello " << name;
return 0;
}
C++ is a widely used middle-level programming language.
When ever you want to perform a set of operations based on a condition If-Else is used.
if(conditional-expression) {
//code
}
else {
//code
}
You can also use if-else for nested Ifs and If-Else-If ladder when multiple conditions are to be performed on a single variable.
Switch is an alternative to If-Else-If ladder.
switch(conditional-expression){
case value1:
// code
break; // optional
case value2:
// code
break; // optional
......
default:
code to be executed when all the above cases are not matched;
}
For loop is used to iterate a set of statements based on a condition.
for(Initialization; Condition; Increment/decrement){
//code
}
While is also used to iterate a set of statements based on a condition. Usually while is preferred when number of iterations are not known in advance.
while (condition) {
// code
}
Do-while is also used to iterate a set of statements based on a condition. It is mostly used when you need to execute the statements atleast once.
do {
// code
} while (condition);
Function is a sub-routine which contains set of statements. Usually functions are written when multiple calls are required to same set of statements which increases re-usuability and modularity. Function gets run only when it is called.
return_type function_name(parameters);
function_name (parameters)
return_type function_name(parameters) {
// code
}