https://drive.google.com/drive/folders/1cgFnguKAojQSi_nITflq-E1BRFg4iTj8 clear all close all clc b = 0.05; g = 9.81; l = 1; m = 1; % initial condition theta_0 = [0 3]; % time points t_span = linspace(0,20,500); % solve ODE [t, results] = ode45(@(t,theta) ode_func(t, theta, b, g, l, m), t_span, theta_0); plot(t,results(:,1))…

Main code clear all close all clc % Input variables x = linspace(0, 0.6, 150); y = linspace(0, 0.6, 150); num_case = 250; % creating an arraay [xx, yy] = meshgrid(x,y); % loop for passing the stalagmite function for i = 1 : length(xx) for j = 1 : length(yy) input_vector(1) = xx(i , j); input_vector(2) = yy(i , j); f(i…

close all clear all clc %Opening the file that is stored in the file Thermo = fopen('THERMO.dat','r') %Getting the first line from the opened file fgetl(Thermo); f1 = fgetl(Thermo); A = strsplit(f1, ' '); Global_low_temperature = str2num(A{2}); Global_mid_temperature = str2num(A{3}); Global_high_temperature = str2num(A{4});…

clear all; close all; clc disp(' RANKINE CYCLE SIMULATOR') disp('1-2 Isentropic Expansion in the Turbine') disp('2-3 Constant Pressure Heat Rejection by the Condenser') disp('3-2 Isentropic Compression in the Pump') disp('4-1 Constant Pressure Heat Addition by the Boiler') % Required Inputs P1 = 30 T1 = 400 P2 = 0.5 disp…