clear allclose all clc %inputl1 = 1;l2 = 0.5; theta1 = linspace(0,90,40);theta2 = linspace(0,90,40); ct=1for 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); %plotingplot([x0 x1],[y0…

clear allclose all clc b = 0.05;g = 9.81;l = 1;m = 0.1; % initial condition theta_0 = [0;5]; % time pointt_span = linspace(0,10,500); % solve ODE[t, results] = ode45(@(t,theta) ode_func(t, theta, b,g,l,m),t_span, theta_0); plot(t,result(:,1))hold onplot(t,result(:,2))ylabel('plot')xlabel('time')    additional…

   function [f] = stalagmite(input_vector) a=input_vector(1);b=input_vector(2); f1_x=[sin((5.1*pi*a)+0.5)].^6;f1_y=[sin((5.1*pi*b)+0.5)].^6;f2_x=exp(-4*log(2)*(((a-0.0667)^2)/(0.64)));f2_y=exp(-4*log(2)*(((b-0.0667)^2)/(0.64))); f=-(f1_x*f1_y*f2_x*f2_y);end       clear allclose allclcx = linspace(0,0.6,150);y…

CODE :- clear allclose allclc % make directorymkdir 'c:usersadmindesktopnasa file parsing' % reading data file - nasa thermodynami% open a file & 'r' open the file for readingf1=fopen('THERMO.dat','r')gt=fgetl(f1); %opting temperature - by reading header informationglobal_temp=fgetl(f1); %splitingt=strsplit(global_temp);…