بررسی کد نوشته شده متلب

با سلام من یک کد متلب نوشته ام که نمودار مقاله را پیدا کنم اما نمی‌دانم اشکال کار کجاست که نمودار مشابه مقاله پیدا نمیشه کسی را نیاز دارم که مسلط به متلب باشد و بتونه اشکال کار را بگه و کد صحیح را به من بده و رسم کنه کد نمودار شماره سه مقاله را ارسال میکنم % Parameters from the article qc = 0.7; % Nonextensive parameter for cold electrons qh = 0.7; % Nonextensive parameter for hot electrons beta = 0.12; % Ratio of temperatures (Tc/Th) delta = 0.24; % Fraction of cold electrons % Derived parameters Sc = qc - 1; Sh = qh - 1; Pc = (3 * qc - 1) / 2; Ph = (3 * qh - 1) / 2; % Define σ and M values sigma_values = [0.13, 0.15, 0.21]; % σ values M_values = [2.09817, *******, *******]; % Corresponding M values % Range of φ phi = linspace(-7, 0, 1000); % φ range from -7 to 0 % Initialize figure figure; hold on; % Initialize arrays to track min and max of V(φ) V_min = inf; V_max = -inf; % Loop through combinations of σ and M for idx = 1:length(sigma_values) sigma = sigma_values(idx); M = M_values(idx); % Calculate R1 and R2 R1 = M + sqrt(3 * sigma); R2 = M - sqrt(3 * sigma); % Compute Sagdeev Potential V(φ) V_phi = zeros(size(phi)); for i = 1:length(phi) cubic_term = (1 / (6 * sqrt(3 * sigma))) * ... (R1^3 - R2^3 + (R2^2 - 2 * phi(i))^(3/2) - (R1^2 - 2 * phi(i))^(3/2)); nonlinear_cold_term = delta * Pc * (1 - (1 + Sc * phi(i))^(Pc / Sc)); nonlinear_hot_term = (1 - delta) / (beta * Ph) * ... (1 - (1 + beta * Sh * phi(i))^(Ph / Sh)); % Total Sagdeev Potential V_phi(i) = cubic_term + nonlinear_cold_term + nonlinear_hot_term; % Track min and max of V(φ) V_min = min(V_min, V_phi(i)); V_max = max(V_max, V_phi(i)); end % Plot the Sagdeev Potential plot(phi, V_phi, 'DisplayName', sprintf('\\sigma = %.2f, M = %.5f', sigma, M)); end % Labels and legend xlabel('\phi (Potential)'); ylabel('V(\phi) (Sagdeev Potential)'); title('Sagdeev Potential V(\phi) for Double Layers'); legend('show'); grid on; % Adjust axis limits xlim([-7, 0]); % φ from -7 to 0 ylim([1.1 * V_min, 0]); % V(φ) from 1.1 * V_min to 0 hold off;