Week 4.1 - Solving second order ODEs

close all clear all clc %-----------------------------------user inputs---------------------% Temp1 = "What is Temperature 1 in celsius? "; T1 = input(Temp1); %turbine inlet temperature Pressure1 = "What is pressure 1 in Bar? "; P1 = input(Pressure1); %turbine inlet temperature Pressure2 = "What is pressure 2 in Bar? ";…

clear all close all clc f1 = fopen('THERMO.dat','r'); fgetl(f1); G_temp = str2num(fgetl(f1)); %Global temparture G_low = G_temp(1); G_mid = G_temp(2); G_high= G_temp(3); for i = 1:3 fgetl(f1); end R = 8.3145; %J/Mol*Km (units) for j = 1:53 line = fgetl(f1); A = strsplit(line,' '); % species species =(A{1}); % Local_temparture…

https://drive.google.com/file/d/1J3BnGoUm0o-joDnuvQDoaQltZbuvokzZ/view?usp=drive_linkAim : To simulate the pendulum motion Code : clear all close all clc %% Main Code L = 1; % Length of the pendulum (m) m = 1; % Mass of the ball (kg) b = 0.05; % Damping coefficient g = 9.81; % Gravity (m/s^2) % Initial conditions…

%Everything done properly but i didnt get the result, I dont know why % Define the arm length L1 = 1; % Length of first link L2 = 0.75; % Length of second link % Define the joint angles theta1 = linspace(0,pi/2,50); % Joint angle 1 in degrees theta2 = linspace(0,pi/2,50); % Joint angle 2 in degrees % Initialize the figure…