close all clear all clc   %inputs l1=1; l2=0.5; theta1=linspace(0,90,10); theta2=linspace(0,90,10);   ct=1; for i=1:length(theta1) THETA1=theta1(i); for j=1:length(theta2) THETA2= theta2(j);   x0=0; y0=0; x1= l1*cosd(THETA1); y1=l1*sind(THETA1); x2= x1+l2*cosd(THETA2); y2= y1+ l2*sind(THETA2);   %plotting…

clear all close all clc   %inputs 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)) hold on plot(t,results(:,2)) xlabel('time') ylabel('plot')…

clear all close all clc   %defining our search space x=linspace(0,0.6,150); y=linspace(0,0.6,150);   %creating a 2 dimensional mesh [xx yy]=meshgrid(x,y);   %evaluating 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,j)=stalagmite(input_vector);…

close all clear clc   %Inputs R= 8.314;   %Reading the data file f1 = fopen('THERMO.dat','r');   %Reading the first line of the file L1 = fgetl(f1);   %Getting the global minimum, maximum and mid temperatures L2 = fgetl(f1); t = strsplit(L2); global_min = str2double(t{2}); global_mid = str2double(t{3});…