Erasing part of contourf diagram (L-shaped diagram using contourf)
try ind = x_points > 2 & y_points < 2; temp(ind) = nan;

27 jours ago | 0

solve on vector equation
Try this (solution exists not for any s1 vector) function main nt = 1; ni = 2; s1 = [1 2 2]; s1 = s1/norm(s1); fun...

environ un mois ago | 0

How to detect these points in red circle from point cloud. The attached data.txt is the point cloud data. Many thanks.
Here is an idea: load data.txt x = data(:,1); y = data(:,2); z = data(:,3); plot3(x,y,z,'.r') axis equal datacursormode...

environ un mois ago | 0

| accepted

Determine cross sectional area from patch data
I think there is no built-in function for your purpose Try to search something like THIS

environ un mois ago | 0

| accepted

slice plot gives black surfaces
Planes are not black. THe mesh is just too fine Add to your code h = slice(z,y,x,abs(Eg(:,:,:)),zslice,yslice,xslice); set(h,...

environ un mois ago | 0

| accepted

ODE solver - division by zero at time boundaries
As I see it, the problem is at t = 0, y(1,2) = 0 Can you replace 0 with 1e-3?

environ un mois ago | 0

| accepted

Solving a system of ODEs with conditionally defined coefficient
replace f=piecewise(0<=t<1, -2*t^2+3*t, t>=1, 1/t); with if t >= 1 f = 1/t; else f = -2*t^2+3*t; end

environ un mois ago | 0

| accepted

How to implement Forward Euler Method on a system of ODEs
I've tried and reached a success par_c.k = 14.39964548; % electronvolt per AngstrÃ¶m par_c.xj = [0 0]; % position of particle j...

environ un mois ago | 0

| accepted

Solving two second order BVPs and the boundary conditions are related
How to set up initial guess? Any ideas? function main %-----------------------------------------------------------------------...

environ un mois ago | 1

| accepted

Position of elements'nodes intersected by a line
Look HERE Check each triangle/element. If intersection exists then pick/highlight this element/triangle

environ un mois ago | 1

| accepted

Plot dB/dt vs. time by ODE45 (with dB/dt on the y axis and time on the x axis)
What I would like to do is also plot dB/dt vs. time pvar = gradient(varsol(:,3)); dt = gradient(tsol); dB = pvar./dt; plot(t...

environ un mois ago | 1

| accepted

Heat transfer differential equation
Call the function this way: [t,y] = ode45(@f,[0,5],1673)

environ un mois ago | 0

2D information to 1D colorplot
Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

environ un mois ago | 0

Quiver with a system of ODEs.
For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

environ un mois ago | 0

Solving nonLinear equation and Linearized Eqaution using ODE45
YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

environ un mois ago | 0

| accepted

How to plot magnitude and directions between 2 points?
try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

environ un mois ago | 0

| accepted

surface fixed point value Z
I want to know the value of Z from the surfc function Here is it

environ un mois ago | 0

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?
Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

environ 2 mois ago | 1

| accepted

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?
alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

environ 2 mois ago | 0

Euler 3D rotation between two vectors
Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

environ 2 mois ago | 0

| accepted

Simultaneous differential equations - derivative of y^3 wrt t
Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

environ 2 mois ago | 0

Numerically derive a continous, non-symbolic function
Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

environ 2 mois ago | 0

3D Point Cloud - Gaussian Curvature
one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

environ 2 mois ago | 0

Runge Kutta for system of eqs
I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

environ 2 mois ago | 0

The problem of finding intersection point of triangle mesh
To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

environ 2 mois ago | 1

| accepted

How can i solve a system which contains matrices
Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

environ 2 mois ago | 1

| accepted

Using a for-loop to score values in different intervals
Correct form if 0 <= x(i) && x(i) < 1

environ 2 mois ago | 0

| accepted